La Ciencia 
para Todos 


La colección La Ciencia desde México, del Fondo de Cul- 
tura Económica, llevó, a partir de su nacimiento en 1986, 
un ritmo siempre ascendente que superó las aspiraciones 
de las personas e instituciones que la hicieron posible: 
nunca faltó material, y los científicos mexicanos desarro- 
llaron una notable labor en un campo nuevo para ellos, 
escribir de modo que los temas más complejos e inaccesi- 
bles pudieran ser entendidos por los jóvenes estudiantes 
y los lectores sin formación científica. 

Tras diez años de trabajo fructífero se ha pensado 
ahora dar un paso adelante, extender el enfoque de la 
colección a los creadores de la ciencia que se hace y pien- 
sa en lengua española. 

Del Río Bravo al Cabo de Hornos y, cruzando el 
océano, hasta la Península Ibérica, se encuentra en mar- 
cha un ejército compuesto de un vasto número de inves- 
tigadores, científicos y técnicos, que desempeñan su 
labor en todos los campos de la ciencia moderna, una 
disciplina tan revolucionaria que ha cambiado en corto 
tiempo nuestra forma de pensar y observar todo lo que 
nos rodea. 

Se trata ahora no sólo de extender el campo de acción 
de una colección, sino de pensar una ciencia en nuestro 
idioma que, imaginamos, tendrá siempre en cuenta al 
hombre, sin deshumanizarse. 

Esta nueva colección tiene como fin principal poner el 
pensamiento científico en manos de los jóvenes que, si- | 
guiendo a Rubén Darío, aún hablan en español. A ellos 
tocará, al llegar su turno, crear una ciencia que, sin des- 
deñar a ninguna otra, lleve la impronta de nuestros 
pueblos. 
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INTRODUCCIÓN 


Han pasado tres décadas desde que se inició una nueva línea 
de investigación científica llamada “teoría del caos”. : 
A diferencia de lo que ocurre en otros campos de la física, 
como la mecánica cuántica, las investigaciones sobre las par- 
tículas fundamentales que constituyen la materia o las teorías 
sobre el origen del Universo, se está intentando aplicar esta 
“ciencia del caos” a muchos eventos vinculados directamente 
con la experiencia humana habitual, y explicar así fenómenos 
tan disímiles como las arritmias en el funcionamiento del' 
corazón, o aspectos de la economía como. las fluctuaciones de 


la bolsa de valores, o también la aparición de la vida sobre la 
Tierra, además del comportamiento de los sistemas físicos di- 
námicos con un número elevado de componentes como pue- 
den ser la atmósfera o un líquido en estado turbulento. 

El físico Joseph Ford, proclamó en un artículo del libro The 
New Physics la nueva ciencia del caos como “el comienzo de la 
tercera revolución de la física del presente siglo”, siendo las 
dos anteriores la teoría de la relatividad y la teoría cuántica. 

Parece aún prematuro, sin embargo, darle la categoría de 
tercera revolución de la física, ya que a diferencia de los otros 
dos casos no es correcto hablar de una “teoría del caos”, teoría 
que aún no existe. Se trata de una nueva y muy promisoria 
manera de aplicar las leyes conocidas de la física, con la ayuda 
fundamental de la computadora, a fenómenos muy variados 
que abarcan, además de los tradicionales en física, a los que se 
presentan en las ciencias biológicas y las ciencias sociales, 
siempre y cuando se les pueda encarar como si se tratara de 
sistemas dinámicos complejos. 

En la física se está trabajando con intensidad creciente en 


, 
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los temas de medios desordenados y de sus propiedades espe- 
cíficas, de gran interés científico y tecnológico. 

Pero hay que evitar las confusiones que se pueden generar 
alrededor de este tema, en particular, por las expectativas que 
puede despertar el estudio del caos en los que Mini en 
otros campos del conocimiento. 

En realidad aparece aquí un problema de interpretación: 
entre quienes no están familiarizados con las ciencias físicas O 

matemáticas —y debido en parte a las declaraciones de algu- 
nos científicos— se ha instalado una especie de mitología del 
Caos o Desorden, que asigna un significado trascendente al 
azar (real o aparente) de la naturaleza, y que proclama la 
muerte definitiva del determinismo, cuando todo indica que 
para los sistemas caóticos sigue siendo válido el determinismo, 
si bien se requiere una descripción probabilística de su com- 
portamiento. 

De manera similar a lo que ocurrió con la teoría de las catás- 
trofes desarrollada por René Thom —que fue un primer in- 
tento de estudiar matemáticamente ciertos fenómenos com- 
plejos— hay quienes esperan que el estudio del caos permita 
desentrañar los misterios de las grandes transformaciones so- 
ciales o de la relación entre las redes neuronales y la psi- 
cología, y también ha desencadenado formidables especula- 
ciones sobre el significado del tiempo y del desorden en el 
Universo. 

Por supuesto que es muy beneficiosa la extensión de los des- 
cubrimientos realizados en un campo del conocimiento a 
otras áreas, pero cuando se trata, por ejemplo, del compor- 
tamiento humano, individual o colectivo, que tiene una com- 
plejidad incomparablemente mayor que la de los sistemas fisi- 
cos, esa extensión debe hacerse con mucha prudencia, y en 
general sólo puede tener un carácter de analogía. 

Así, puede resultar fructífero encarar la psicología de un 
grupo familiar aplicando ciertas pautas que tienen analogías 
con las de sistemas dinámicos físicos, pero difícilmente pueda 
trabajarse en este tema aplicando la matemática del caos, y 
buscando dimensiones fractales y atractores extraños. 

Debería evitarse la utilización de un lenguaje que parezca 
atribuir un alcance mágico al “caos”: en los textos científicos, 
este concepto tiene un sentido preciso, que nos remite a fenó- 
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menos complejos, particularmente difíciles de formular mate- 
máticamente, pero que no manifiestan, en principio, relación 
alguna con el Caos primordial concebido por las antiguas 
mitologías. 

Es importante dejar en claro que las leyes fundamentales de 
la física siguen rigiendo y que el hecho de que se utilicen, 
como veremos, las características estadísticas para predecir com- 
portamientos no es un “drama epistemológico” como algunos 
han sugerido. 

La fuerte carga emocional que tiene la palabra caos es en 
parte la causante de las confusiones antes mencionadas, y ha 
contribuido a ello el hecho de que aún no está definitivamen- 
te establecido el nombre de esta nueva disciplina. 

Hace treinta años, cuando comenzó a desarrollarse, se ha- 
blaba de la “ciencia del caos”, que pronto pasó a denominarse 
“caos determinista”, para diferenciarlo del caos producto del 
puro azar. Actualmente tiende a afianzarse la palabra “comple- 

- jidad”, que designa el estudio de los sistemas dinámicos que 
están en algún punto intermedio entre el orden en el que nada 
cambia, como puede ser el de las estructuras cristalinas, y el 
estado de total desorden o caos como puede ser el de un gas 

_ Ideal en equilibrio termodinámico. 

Los fenómenos de “caos determinista” o de “complejidad” 
se refieren a muchos sistemas que existen en la naturaleza 
cuyo comportamiento va cambiando con el transcurrir del 
tiempo (sistemas dinámicos). Dichos fenómenos aparecen 
cuando los sistemas se hacen extremadamente sensibles a sus 
condiciones iniciales de posición, velocidad, etcétera, de modo 
que alteraciones muy pequeñas en sus causas son capaces de 
provocar grandes diferencias en los efectos. Como consecuen- 
cia de ello no es posible predecir con exactitud cómo se com- 
portarán dichos sistemas más allá de cierto tiempo, por lo que 
parecen no seguir ninguna ley, cual si estuviesen regidos por 
el azar. 

Pero los investigadores han encontrado que los sistemas 
dinámicos en estas condiciones presentan pautas de regulari- 
dad colectiva aunque no sea posible distinguir el compor- 
tamiento individual de cada uno de sus componentes. 

Se ha comprobado que hay ciertas características comunes 
que permiten incluir en el estudio de procesos complejos no 
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sólo los sistemas físcos y químicos inertes sino también organis- 
mos vivos, abordados todos mediante herramientas matemáti- 
cas comunes. La herramienta fundamental es, por supuesto, la 
computadora, sin la cual hubiera sido imposible desarrollar 
este nuevo enfoque de los sistemas dinámicos. De manera si- 
milar al impulso que dio a la ciencia la utilización del telesco- 
pio y el microscopio en los siglos XVI y XVII, el uso de esta má- 
quina facilita enormemente la comprobación de las teorías 
mediante la experimentación, gracias al inmenso incremento 
de la capacidad de cálculo y a la posibilidad de hacer simula- 
ciones de los procesos reales y de crear modelos de sistemas 
complejos. No se podrían haber desarrollado conceptos como 
los de atractores extraños, fractales, o complejidad algorítmica 
sin su existencia. 

De hecho, todo este vasto campo de fenómenos no lineales, 
con alta sensibilidad a las condiciones iniciales, no era desco- 
nocido por los grandes matemáticos y físicos del siglo pasado; 
pero en esa época la solución de los correspondientes sistemas 
de ecuaciones requería cálculos numéricos tan engorrosos que 
los hacían impracticables, por lo que se debió esperar a la se- 
gunda mitad de este siglo, con la aparición de las computado- 
ras veloces, para poder encararlos. 

Uno de los resultados más positivos debidos al surgimiento 
de este nuevo campo de investigación es que se han formado 
grupos interdisciplinarios —integrados, por ejemplo, por bió- 
logos, físicos, matemáticos, o sociólogos, economistas y exper- 
tos en computación— para estudiar los problemas inherentes 
a sistemas dinámicos complejos. Estos abarcan desde líquidos 
turbulentos hasta sistemas ecológicos o los modelos económi- 
cos de las sociedades. 

Varios centros de alto nivel han surgido, en particular en 
Estados Unidos, Rusia y Francia; así se pueden mencionar el 
Santa Fe Institute en Nuevo México, y departamentos dedica- 
dos al tema en Los Alamos National Laboratory (Center for 
Non-Linear Studies), el Instituto de Tecnología de Georgia, la 
Universidad de California en Berkeley, y el Centro de Investi- 
gaciones en Saclay, Francia, además del grupo en la Universi- 
dad Libre de Bruselas, Bélgica. 

Un papel preponderante en el estudio de estos fenómenos 
fue el cumplido por la escuela rusa de físicos y matemáticos, 
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como Landau, Kolmogorov, Andronov, Lyapunov y otros, quie- 
nes desarrollaron las técnicas necesarias mucho antes de que 
el caos determinista pasara a ser un tema de moda. Dicha es- 
cuela sigue haciendo grandes contribuciones a esta impor- 
tante área. 

En los últimos años el avance de esta actividad se manifiesta 
por la creciente cantidad de congresos sobre sus diversos as- 
pectos. Las posibilidades de aplicación a las más diversas cien- 
cias se pusieron de relieve en la Primera Conferencia de Caos 
Experimental, realizada en Estados Unidos en octubre de 1991, 
donde se demostró que hay casos en los que se puede aprove- 
char el caos en lugar de evitarlo, y obtener entonces sistemas 
de más flexibilidad que los que por ser ordenados tienen un 
comportamiento “bueno”, o sea predecible. 

Este libro tiene como propósito describir las características 
fundamentales de los sistemas complejos y de los métodos que 
se utilizan para estudiar su comportamiento. 


Deseo, por último, expresar mi agradecimiento a Marcos Sara- 


ceno por sus valiosas sugerencias acerca del desarrollo de este 
tema y de la presentación del texto. 
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I. El universo, ¿funciona como un reloj? 


Es COMÚN a muchos mitos de la creación el concepto de una 
batalla eterna entre Orden y Caos, lo que sugiere que estos dos 
conceptos están profundamente arraigados en la mente hu- 
mana. 

Para la humanidad primitiva, la Naturaleza era el Caos, un 
ente caprichoso, sujeto a los antojos de dioses poderosos e in- 
descifrables, cuyas acciones no se podían prever. 

Así, las cosmologías de muchas culturas imaginaron un esta- 
do inicial de Desorden, o Caos, del que surgieron las cosas y 
los seres. 

En La teogonía, Hesíodo dice: “El Caos fue lo primero, y lue- 
go la Tierra”. Caos es una palabra de origen griego, equivalen- 
te a abismo, o también a una entidad sin forma; cosmos, por 
su parte, designa el orden y, por extensión, el Universo. 

Un mito chino de la creación dice que del caos surgieron los 
dos principios, Yin y Yang, que con sus aspectos femenino y 
masculino crearon luego todo el Universo. 

En la Biblia leemos: “La Tierra estaba informe y vacía, y 
había oscuridad sobre la haz del abismo; y el Espíritu de Dios 
se cernía sobre la haz de las aguas”. 

El caos es entonces la sustancia primordial sin forma, con la 
que el Creador moldea el mundo. Pero los seres humanos ne- 
cesitan descubrir un orden.en la naturaleza, buscar las leyes 
detrás de su complejo comportamiento que les permitan co- 
nocer la duración de los días y las noches, las fases de la Luna, 
la época de las cosechas. Se introduce así la noción de orden 
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ante la necesidad de prever con anticipación, de predecir, in- 
dispensable para la supervivencia. 

Lentamente, a lo largo de muchos milenios, el hombre fue 
descubriendo en la naturaleza comportamientos ordenados, 
regulares, los que aprendió a registrar, predecir y explotar. 

En el pensamiento griego aparecen dos visiones diferentes 
acerca de qué es lo importante para nuestra comprensión del 
Universo, y que se mantienen en esencia en Occidente hasta 
nuestros días: el monismo estático de Parménides y el ser en 
perpetuo movimiento concebido por Heráclito. 

Platón tomó de Parménides el énfasis en el orden, en los 
aspectos inmutables de la realidad, y de Pitágoras el estudio de 
las leyes matemáticas y geométricas que son las expresiones 
de las “Formas”, ideas eternas, de las que el hombre percibe 
sólo las sombras. Además, como Heráclito, tenía presente la 
inestabilidad y el incesante fluir de todo lo que se manifiesta 
en la naturaleza. 

Con este enfoque, el pensamiento platónico busca algo que 
está detrás del proceso de cambio y de la noción del paso del 
tiempo, para poder hacer una descripción inteligible de la rea- 
lidad. El Cosmos es racional, y si en el comienzo existía el Caos, 
fue el Demiurgo, el Supremo Ordenador, quien dispuso la 
transformación del Universo desde el estado de desorden al de 
orden. Este Ser divino es matemático, y el orden que establece 
también lo es. El campo del conocimiento donde este con- 
cepto queda más claramente expresado es el de la geometría, y 
el pensamiento platónico hace una clara distinción entre las 
Formas, que son idealizaciones, son perfectas, y los objetos tal 
como aparecen en la naturaleza, en la que nunca se encon- 
trarán círculos, esferas o planos perfectos. 

De ahí la clara aversión de los pensadores platónicos a uti- 
lizar herramientas para estudiar la geometría —sólo toleraban 
el uso de la regla y el compás—, pues cualquier objeto material 
es imperfecto, pálido reflejo de las Ideas. 

Por acción de la “causa errante”, la materia se resiste a ser 
modelada por las ideas, e impide así que en ella encarnen del 
todo las formas puras de la geometría. El resultado es entonces 
un Cosmos agitado por innumerables pequeñas convulsiones, 
a las que modernamente podríamos denominar movimientos 
de torbellinos. 
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Para Platón, en consecuencia, hay una jerarquía formada 
por tres niveles fundamentales: las Ideas y Formas matemáti- 
cas, que son el modelo perfecto de todas las cosas; el Caos ori- 
ginal, y un estado intermedio, que es nuestro mundo imper- 
fecto, complejo, resultado de la tarea que hizo el Demiurgo 
partiendo del Caos y modelándolo sobre la base de las Ideas. 

En este orden jerárquico, el valor supremo es el de las Ideas 
y Formas matemáticas, que expresan las cualidades divinas de 
simplicidad, armonía, regularidad. 

- Esta filosofía fue incorporada al pensamiento cristiano me- 
dieval, y los científicos a partir del Renacimiento han coincidi- 
do con sus aspectos fundamentales. Galileo, dos mil años des- 
pués de Platón, dice que el libro de la naturaleza está escrito 
en lenguaje matemático, y en su época se sientan las bases para 
el método científico, que consiste en buscar las leyes eternas 
que rigen un fenómeno natural, las que deben ser formulables 
en el lenguaje matemático, e interrogar a la naturaleza para 
descubrir si los hechos corroboran o niegan la teoría propuesta. 

Para los pensadores como Galileo, Kepler, Newton o Einstein, 
todo ocurre como si Dios (o la naturaleza) hubiese elegido el 
orden de las “formas puras”, las que cuanto más simples, se 
consideran tanto más bellas y verdaderas, y el ideal de la cien- 
cia es hallar a través de la razón ese orden y regularidad que 
están detrás del desorden aparente de la naturaleza. 

Como expresara con tanta claridad Albert Einstein: 


Reconocemos en la base de todo trabajo científico de cierta enver- 
gadura una convicción comparable al sentimiento religioso, pues 
acepta un mundo basado en la razón, un mundo inteligible. Esta 
convicción, ligada a un sentimiento profundo de una razón supe- 
rior que se revela en el mundo de la experiencia, expresa para mí 
la idea de Dios. 


La ciencia postula así que detrás de la complejidad del mun- 
do hay leyes matemáticas que muestran una armonía subya- 
cente, en la que no hay lugar para el desorden y lo imprevisi- 
ble; éstos se disiparán al ser iluminados por la luz de la razón. 

Pero hay otra visión que, hasta hace poco, ha sido menos 
difundida que la búsqueda platónica de aquello que es inva- 
riante en la naturaleza. Es la perspectiva de Aristóteles, que 
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pone su énfasis en el cambio, en los procesos observables en el 
mundo, en lugar de las Formas invariantes detrás de ellos que 
no pueden ser observadas. | 

En esta visión tuvieron mucha influencia sus observaciones 
del mundo viviente, para las que estaba muy bien dotado. Para 
los abogados de este enfoque, la naturaleza es semejante a un 
organismo, complejo, cambiante, y aceptan como reales y no 
como aparentes sus aspectos de desorden e-imprevisibilidad, 
tratando de comprender sin eliminar esas características, así 
como nadie intentaría transformar un organismo vivo en una 
máquina tan previsible como un reloj. 

Volveremos sobre este enfoque más adelante, pero detallare- 
mos ahora las características del método científico tal. como 
fuera aplicado a partir de la obra de los grandes pensadores 
como Galileo, Descartes, Huyghens y Leibniz. 

El método científico aplicado por Galileo, tuvo su expresión 
matemática gracias a Descartes, quien enfatizó la necesidad de 
analizar, es decir, dividir lo que se examina en sus componen- 
tes más simples, para luego recomponerlo en una síntesis que 
permita comprender el fenómeno con certeza. Creó además la 
geometría analítica basándose en su introducción de las coor- 
denadas, identificando así al espacio con una inmensa grilla a 
cada punto de la cual se le puede asignar un valor numérico, 
con lo que logró unir la geometría con el álgebra. 


EL MENSAJE DE NEWTON 


Hace tres siglos se publicó la monumental obra de Isaac New- 
ton, Principios matemáticos de filosofía natural, cuyo mensaje ha 
sido decisivo para la cultura de Occidente. 

Según éste: 


El Universo es ordenado y predecible; tiene leyes expresables en 
lenguaje matemático, y podemos descubrirlas. 


Como hemos visto, para los científicos a partir del Renaci- 
miento se reafirma el enfoque platónico, ya que hay un orden 
en el Universo detrás de su aparente complejidad, con leyes 
simples que contienen aspectos inmutables, expresados por las 
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famosas magnitudes que se conservan invariantes en la física: 
energía total, impulso, carga eléctrica. 

Las leyes descubiertas por Newton parecen sencillas, y su 
aplicación al comportamiento de los cuerpos permitió descri- 
bir con precisión el movimiento de los astros en el firmamen- 
to, la caída de los cuerpos, etcétera. 

Newton formuló sus leyes mediante ecuaciones matemáti- 
cas, que relacionan las magnitudes que podemos medir de un 
cuerpo, tales como su posición y su velocidad, con la forma en 
que estas varían con el tiempo. 


LEYES DE NEWTON DE LA DINÁMICA 


1) Inercia: Todo cuerpo material que no tiene aplicada una fuerza 
permanece en reposo o se mueve en línea recta con velocidad 
uniforme. 

2) Fuerza: Cuando se aplica una fuerza a un cuerpo libre su im- 
pulso cambia en el tiempo proporcionalmente a dicha fuerza, 
y la dirección de su movimiento es la de la línea de acción de la 
fuerza. 

3) Acción y reacción: Para dada acción ejercida sobre un cuerpo hay 
siempre una reacción igual y opuesta. 


La primera ley fue la formulación del descubrimiento de 
Galileo, quien notó que se debe atender no a la velocidad de un 
cuerpo sino al cambio de esa velocidad con el transcurso del 
tiempo, poniendo fin así a la creencia aristotélica que había 
bloqueado el avance de la física por muchos siglos. 

Newton aplicó además un concepto que ha sido esencial 
para el método científico: el de aislar idealmente el sistema 
dinámico que se desea examinar del resto del universo del que 
forma parte. Esto permite que su comportamiento se pueda 
comprender, ya que no es necesario considerar todas sus in- 
finitas relaciones con el universo, lo que sólo sería posible para 
un ser infinito. Basta entonces que se consideren sólo aquellas 
características del sistema que son relevantes al fenómeno que 
se desea estudiar. Así, la primera ley pregunta qué se puede 
decir del movimiento de un cuerpo aislado, o'sea al que no se 
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le está aplicando una fuerza. Aristóteles había dicho que de- 
berá mantenerse en reposo. 

Newton, al igual que Galileo antes, establece que, en esa 
situación, el cuerpo puede estar en reposo q puede moverse 
en línea recta con velocidad uniforme. 

Si un cuerpo cae hacia la tierra, esto se debe a que sobre él 
actúa la fuerza de gravedad, y por lo tanto su velocidad no 
puede ser uniforme sino que será cada vez.mayor; por eso es 
más peligroso caer desde una altura grande que de una me- 
nor. 

¿Hay algo que se mantenga constante en este proceso de caí- 
da? Sí, la aceleración, es decir, la velocidad con la que aumenta 
la velocidad del cuerpo al caer. 

La segunda ley formula este concepto, que se expresa ma- 
temáticamente como F = mxXa (la fuerza Faplicada al cuerpo 
es proporcional a la aceleración a de su movimiento, según una 
constante m, la masa o cantidad de materia del cuerpo). 

Para cualquier sistema dinámico del Universo, las leyes de 
movimiento se pueden expresar como F' = mx a, no importa 
cuál sea el origen de la fuerza aplicada. Puede tratarse de la 
fuerza de gravedad, la eléctrica o la magnética, y para todas 
vale la misma asombrosamente simple ecuación. 

La velocidad con que cambia una magnitud se determina 
por la diferencia entre sus valores para dos tiempos sucesivos, y 
de allí el término “diferencial” que aparece en el análisis ma- 
temático: las ecuaciones de Newton involucran velocidades de 
cambio y son, por lo tanto, ecuaciones diferenciales. 

Las ecuaciones algebraicas se distinguen de las diferenciales 
porque no involucran velocidades de cambio. No son siempre 
fáciles de resolver. Pero resolver ecuaciones diferenciales es, 
en general, mucho más dificultoso, y resulta realmente asom- 
broso que tantas ecuaciones importantes por sus aplicaciones 
a la física tengan una solución definida. 

Un principio matemático básico de las ecuaciones diferen- 
ciales es que su solución, es decir, su integración, está determi- 
nada univocamente, y da un único resultado para cada conjun- 
to de datos numéricos que se introduzca en las ecuaciones; si, 
por ejemplo, se quiere conocer qué altura alcanzará un pro- 
yectil, introduciendo en las ecuaciones los datos de la veloci- 
dad inicial, ángulo del cañón, etcétera, se obtiene un resulta- 
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do que define sin ambigúedad esta altura; son por lo tanto 
ecuaciones deterministas: hay un único efecto para cada causa. 

La importancia de estas ecuaciones diferenciales estriba en 
que se les puede aplicar a un sistema dinámico, es decir a todo 
proceso que cambia en el tiempo. 

Muchos físicos y matemáticos se han maravillado ante el 
hecho de que estas ecuaciones parezcan tan efectivas para des- 
cribir la estructura del mundo físico. 

En vista de la complejidad del mundo que nos rodea, es 
realmente notable que haya fenómenos naturales que admitan 
una descripción mediante leyes físicas simples. ¿Cómo es esto 
posible? 

Fue Isaac Newton quien tuvo la visión que abrió el camino 
para que las ciencias naturales tuvieran tanto éxito en estos 
tres últimos siglos. 

Éste se debe en parte a que inicialmente restringió su aten- 
ción al estudio de sistemas naturales simples con sólo unos 
pocos componentes. 


LEYES FÍSICAS Y CONDICIONES INICIALES 


De acuerdo con el enfoque que hizo Newton de la mecánica, 
un sistema material puede ser dividido conceptualmente en: 
1) las “condiciones iniciales”, que especifican su estado físico 
en cierto tiempo inicial (estas condiciones pueden ser, por 
ejemplo, la posición y la velocidad de un proyectil, o de la 
Luna respecto de la Tierra); 2) las “leyes físicas” o haturales, 
que especifican cómo cambia ese estado. 

Las condiciones iniciales son usualmente muy complicadas, 
una complicación que refleja la complejidad del mundo en 
que vivimos. Las leyes naturales, por otra parte, pueden ser y 
son más bien simples, y se expresan mediante las ecuaciones 
diferenciales. Esta división —leyes y condiciones iniciales— se 
mantiene hasta hoy. 

En la práctica sólo se pueden resolver en forma estricta aque- 
llas ecuaciones que representan las leyes físicas simples para 
sistemas con condiciones iniciales simples: el disparo de un 
proyectil, el movimiento de la Tierra alrededor del Sol sin 
tener en cuenta la influencia de los demás planetas. Esto es, 
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ante la infinita complejidad de la naturaleza, que hace que 
cualquiera de sus componentes esté vinculado con los demás 
por una inmensa cantidad de relaciones, se hace una abstrac- 
ción, considerando idealmente el sistema que se desea estu- 
diar como separado del resto, y seleccionando aquellas ca- 
racterísticas del sistema que parecen ser lo suficientemente 
importantes frente a las otras que no producen casi efectos 
sobre el fenómeno que se está examinando: Así, para calcular 
la trayectoria de un proyectil sólo se considera cómo influye 
sobre él la atracción de la gravedad de la Tierra, pues la atrac- 
ción que ejercen la Luna o el Sol es tan pequeña que no se ne- 
cesita incluirla en las condiciones iniciales. Del mismo modo, 
para estudiar el movimiento de la Luna alrededor de la Tierra 
no se tendrá en cuenta la atracción de las estrellas. 

Se tienen entorices condiciones iniciales simples porque se 
han seleccionado aquellas variables que son las que más afec- 
tan el fenómeno en estudio, y además éstas obedecen a leyes 
que se expresan con ecuaciones en las que pequeñas variacio- 
nes en las condiciones iniciales dan soluciones que difieren 
poco entre sí. De esta manera, el hecho de que las condiciones 
iniciales se conozcan en general con cierto margen de error 
afecta relativamente poco el resultado que se puede esperar de 
estas ecuaciones. Por ejemplo, si se dispara un proyectil con 
un fusil, el cálculo de dónde golpeará al blanco partiendo de 
una cierta posición de su cañón o de otra ligeramente distinta 
producirá una diferencia proporcionalmente pequeña en el 
resultado. 

Otro aspecto importante, también estudiado por los funda- 
dores de la mecánica, fue el de la reversibilidad en el tiempo 
de las trayectorias en los sistemas dinámicos: las ecuaciones 
muestran que si se invierte el sentido de las velocidades de 
todos los componentes del sistema, reemplazando y por —v, el 
resultado es matemáticamente equivalente a cambiar el tiem- 
po £ por —t, como si el sistema fuese “hacia atrás” en el tiempo. 
Esto es la forma matemática de expresar que si a partir de un 
cierto instante hay un cambio en el sistema dinámico, otro 
cambio, definido mediante la inversión de las velocidades de 
los componentes, puede restaurar las condiciones originales. 
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TODO PROCESO FÍSICO ESTÁ DETERMINADO Y POR ESO 
ES POSIBLE PREDECIR SU COMPORTAMIENTO 


Ateniéndose a estas reglas de juego, se llega a la conclusión de 
que en un sistema que responde a las leyes de la mecánica clá- 
sica, y que, por lo tanto, es determinista, si se conocen las posi- 
ciones y las velocidades de sus componentes en un instante 
dado, se pueden calcular las posiciones y velocidades en todo 
instante siguiente o anterior. 

Así si en el sistema dinámico formado por dos bolas de billar 
conocemos en un momento dado inicial la posición y veloci- 
dad de cada una, podremos mediante las ecuaciones diferen- 
ciales de Newton predecir los movimientos respectivos desde 
que se les da el impulso inicial hasta que abandonan el paño. 

Más aún, al finalizar el siglo xvin se llegó a la convicción de 
que si se conociera la posición y velocidad de cada uno de los 
planetas que componen el sistema solar en un instante dado, 
se podría calcular la posición de los mismos en el futuro, y tam- 
bién su posición en el pasado, mediante ecuaciones que deter- 
minan sin ambigúedad las trayectorias, Precisamente la aplica- 
ción de las ecuaciones al movimiento de los astros, la mecánica 
celeste, fue lo que significó el mayor triunfo para Newton. 

Las diversas disciplinas científicas que se fueron desarrollan- 
do en los siglos siguientes a partir de la teoría de Newton estu- 
diaron otras magnitudes además de la posición y la velocidad 
de los cuerpos. Pero el procedimiento, consistente en introdu- 
cir los valores de dichas magnitudes para un tiempo inicial en 
una ecuación matemática que, una vez resuelta, determina esos 
valores para cualquier otro tiempo, es el mismo ya sea que se 
trate de la configuración de un átomo, el movimiento de un 
cometa, la temperatura de un gas o el voltaje en un circuito 
eléctrico. 

Al finalizar este proceso de desarrollo científico parecía que 
el comportamiento de todo el Universo podría llegar a expre- 
sarse matemáticamente, ya que está determinado por las leyes 
mencionadas, inmutables, que dictan el movimiento de cada 
partícula en forma exacta y para siempre, consistiendo la tarea 
del científico en aplicar estas leyes a los fenómenos particu- 
lares. 

El esquema newtoniano posibilitó así la construcción de la 
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majestuosa estructura de la mecánica clásica, que dio a la hu- 
manidad la visión de un Universo ordenado y predecible. 


EL DEMONIO DE LAPLACE 


Esta revolución del pensamiento tuvo su más clara expresión 
con Pierre Simon de Laplace, quien en la época de Napoleón 
expresó, en su Teoría analítica de las probabilidades: 


Debemos considerar el estado presente del Universo como el efec- 
to de su estado anterior y como la causa de su estado futuro. Una 
Inteligencia que, por un instante, conociese todas las fuerzas de 
que está animada la naturaleza y la situación respectiva de los se- 
res que la componen, si además fuese lo bastante profunda para 
someter a éstos al análisis [matemático], abrazaría en la misma 
fórmula a los movimientos de los más grandes cuerpos del Univer- 
so y del átomo más ligero: nada sería incierto para ella y el de- 
venir, como el pasado, estaría presente ante sus ojos. 

El espíritu humano ofrece, en la perfección que ha sabido dar a 
la astronomía, un pálido ejemplo de esta Inteligencia. Sus descu- 
brimientos en mecánica y en geometría, junto con el de la gravi- 
tación universal, lo han puesto en condiciones de abrazar en las 
mismas expresiones analíticas a los estados pasados y a los futuros 
sistemas del mundo. 


Se tiene así la visión del Universo como un mecanismo gi- 
gantesco que funciona “como un reloj” (no un reloj electrónico 
digital, sino de los clásicos, formados por piezas móviles, resor- 
tes y engranajes). Dicho mecanismo es, por encima de todo, 
absolutamente determinista y, por lo tanto, predecible: está 
gobernado por leyes eternas que hacen que bajo circuns- 
tancias idénticas resulten siempre cosas idénticas. Y si las cir- 
cunstancias, en vez de repetirse en forma idéntica, cambian 
levemente, el resultado cambiará también en forma propor- 
cionalmente leve. 

Un especialista en estas leyes de la mecánica que conozca 
sus características y su estado en un momento dado puede, en 
principio, establecer qué hará exactamente en cualquier mo- 
mento del pasado y del futuro. Si en la práctica esto se puede 
aplicar sólo para sistemas relativamente simples, relojes, má- 
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quinas, planetas, y hay muchos objetos en el mundo que pare- 
cen tener comportamientos irremediablemente desordenados, 
caóticos, impredecibles, se trata sólo de algo aparente, y a me- 
dida que se vaya perfeccionando el análisis matemático y se 
descubran las leyes físicas ocultas correspondientes, llegará el 
día en que ese aparente caos desaparecerá. 

Desde este punto de vista el futuro está rígidamente deter- 
minado desde el comienzo del Universo. El tiempo deja de 
tener mucho significado físico, ya que es como si esa Inteligen- 
cia propuesta por Laplace que muchos llamaron el “demonio 
de Laplace” (del griego “dáimon”, divinidad secundaria inter- 
media entre los dioses y el hombre), tuviera toda la historia del 


FIGURA Ll. Pierre Simon de Laplace leyendo su Mecánica Celeste. Tomado de W. F. 
Freeman, Mathematics in the Modern World, San Francisco, Morris Kline. 
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Universo grabada en una película cinematográfica, la cual se 
puede contemplar yendo hacia adelante o atrás en el tiempo. 
Como este tiempo es reversible, sólo marca la dirección en 
que se observa un proceso que no puede modificarse. 


TK. Donde aparece el caos en la máquina 


DURANTE los siglos XVIII y XIX se aplicó la mecánica de Newton 
con impresionante éxito. El punto de vista mecanicista se popu- 
larizó y, combinado con el método experimental, dio un gran 
impulso a la física, la química y la biología. Además se transfor- 
mó en la base de nuevas teorías políticas, económicas y sociales. 

La antigua concepción del caos como señor de la naturale- 
za, donde las cosas sucedían por azar, por capricho, sin rela- 
ción alguna entre causa y efecto, cedió paso a la visión de un or- 
den en el mundo tan determinado como el que muestra un 
buen reloj suizo. 

Pero parecería que hay procesos cíclicos en la historia, o 
más bien similares a una espiral circular donde un ciclo no se 
repite exactamente, sino que pasa a un nuevo nivel. 

Algo así ocurre con el tema de orden y caos en nuestra vi- 
sión del Universo. El caos primitivo fue reemplazado por el 
orden newtoniano. Sin embargo, a medida que el conocimien- 
to de la naturaleza se afianzaba, aparecían dificultades para el 
modelo mecanicista. En la segunda mitad del siglo x1x quedó 
claro cuáles eran los límites de la mecánica clásica: su validez 
no alcanzaba para velocidades extremadamente grandes o para 
el mundo extremadamente pequeño. Producto de esta crisis 
surgieron, en la segunda mitad del siglo xx, dos nuevas ramas 
de la física que estudian, respectivamente, la teoría de la rela- 
tividad y la mecánica cuántica. 


La teoría de la relatividad marcó el límite de la validez de las ecua- 
ciones de Newton, las que deben ser corregidas cuando se está 
frente a velocidades cercanas a la de la luz. 
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La mecánica cuántica establece, mediante el principio de incerti- 
dumbre, un límite a la precisión con que se pueden medir simul- 
táneamente variables tales como la posición y la velocidad de una 
partícula atómica. 


Hoy, los científicos de numerosas disciplinas comienzan a 
convencerse de que hay un tercer límite a la posibilidad de co- 
nocimiento de la naturaleza, y que además es válido para el 
mundo de nuestra experiencia cotidiana: en muchas circuns- 
tancias no sólo no se puede predecir el comportamiento de los 
componentes individuales en sistemas dinámicos complejos, 
que involucran las interacciones de gran cantidad de compo- 
nentes, sino que lo mismo puede ocurrir aun en el caso de sis- 
temas simples, formados por unos pocos componentes que es- 
tén sometidos a la acción de dos o más fuerzas. 

Todos sabemos que el mundo en que vivimos es complejo, y 
nadie se sorprende del poco alcance que tienen las prediccio- * 
nes sobre la economía de un país, o el tiempo meteorológico, 
o el comportamiento de cualquier ser humano o de los orga- 
nismos vivos en general. Esto se ha tomado siempre como la 
indicación de que si es que hay leyes para este mundo tan com- 
plicado, tendrían que ser complicadas, y no como las que ri- 
gen para los sistemas dinámicos que ha estudiado la física, que 
exhiben características de orden y predecibilidad. Pero resulta 
ahora que aun sistemas físicos simples, sometidos a leyes sim- 
ples, pueden tener un comportamiento caótico, no previsible. 

De manera que se nos plantea el dilema: ¿nuestro mundo 
complejo está gobernado por leyes simples que iremos descu- 
briendo progresivamente mediante los métodos desarrollados 
por la ciencia, de acuerdo con una visión que hemos llamado 
platónica? ¿O adoptaremos una visión que recuerda a la aris- 
totélica, en cuanto a poner el énfasis en los procesos de cam- 
bio con el transcurrir del tiempo, aceptando que en muchos 
casos no se pueden predecir con exactitud los comportamien- 
tos de los procesos mediante leyes simples que rigen detrás de 
los fenómenos? 

Para responder este dilema es necesario examinar en qué 
condiciones parece ser válido uno u otro enfoque. 
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SENSIBILIDAD A LAS CONDICIONES INICIALES 


Muchos sistemas dinámicos, cualquiera sea su naturaleza (físi- 
cos, químicos, electromecánicos, biológicos), son extremada- 
mente sensibles a los valores de sus condiciones iniciales, tales 
como la posición, la velocidad, etcétera.Esto pone un límite a 
la posibilidad de predecir el estado futuro del sistema, ya que, 
como hemos visto, dicha predicción está basada en la suposi- 
ción de que las causas pequeñas producen efectos también pe- 
queños y que, por lo tanto, un pequeño cambio en los valores 
iniciales de las ecuaciones diferenciales que describen el com- 
portamiento del sistema, producirá un cambio proporcional- 
mente pequeño en la solución de estas ecuaciones que nos da 
a conocer el estado futuro. 

Corresponde aquí distinguir entre ecuaciones diferenciales 
lineales y no lineales. La solución de una ecuación diferencial 
se denomina integración; la clase más importante entre las 
ecuaciones integrablés es la de las ecuaciones lineales. La más 
simple de entre ellas depende de una sola variable, y su solu- 
ción queda representada gráficamente mediante una línea 
recta; de ahí el nombre de ecuación lineal. 

La familia de las ecuaciones lineales tiene la característica 
de que las soluciones que se obtienen al resolverlas para diver- 
sos valores numéricos de las variables, se pueden sumar entre 
sí, dando como resultado también una solución. Un ejemplo 
"simple es el de la ecuación lineal de onda, que describe el mo- 
vimiento sobre una superficie líquida de ondas de pequeña am- 
plitud. La ecuación tiene muchas soluciones diferentes, cada 
una con diferentes amplitudes y longitudes de onda, y éstas se 
pueden sumar produciendo así una nueva solución de la ecua- 
ción. Esto expresa matemáticamente el hecho físico de que, 
como se ve en el agua en un lago, varias ondas diferentes se pue- 
den superponer, y estas superposiciones corresponden tam- 
bién a una solución que es la suma de soluciones de la ecua- 
ción lineal de onda. 

En general, las ecuaciones lineales son mucho más fáciles de 
resolver que las no lineales, y por eso han sido las más estudia- 
das; además, cuando un fenómeno físico requería ser expre- 
sado mediante una ecuación no lineal de difícil solución, el 
procedimiento usual era linealizarla eliminando aquellos tér- 
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minos que influían menos, es decir, se hacia una aproximación 
lineal. 

Pero en la naturaleza, la mayoría de los fenómenos son ex- 
presables mediante ecuaciones no lineales. Las ecuaciones con 
buen comportamiento de la mecánica clásica, tales como las 
que determinan el movimiento de la Luna, y que permiten 
predecir con exactitud el futuro del sistema, son las excepcio- 
nes, no la regla. 

Recién fue posible encarar la resolución de cualquier ecua- 
ción no lineal con la aparición de las computadoras, ya que 
éstas permiten el análisis numérico de todo tipo de ecuación 
—lineal o no lineal—, por complicada que sea. 

Hoy se puede entonces, gracias a las computadoras, encarar 

el estudio de los sistemas dinámicos cuyo comportamiento 
responde a ecuaciones no lineales, y que son precisamente los 
que "presentan sensibilidad a las condiciones iniciales. 

“Como ejemplo de esta sensibilidad consideremos un sistema 
de comportamiento no lineal, e inestable: un cono de metal pa- 
rado sobre su vértice, similar al de la figura 1.1. Por más que 
hagamos vertical a su eje, terminará por tumbarse, y el lado 
sobre el que caiga dependerá de diferencias minúsculas que 
rompen el equilibrio: un leve soplo, una pequeñísima mota de 
polvo. Predecir de qué lado caerá el cono requeriría el cono- 
cimiento preciso de todas las fuerzas a las que está sometido en 
el momento inicial de equilibrio, lo que equivaldría al imposible 
de introducir la totalidad de una inmensa cantidad de paráme- 
tros como condiciones iniciales en las ecuaciones de movimiento. 


FIGURA 11.1. 
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Otro ejemplo de sistema inestable de comportamiento no 
lineal es el de la ladera cubierta de nieve de una montaña, que 
puede estar en una situación tal que la energía producida al 
lanzar un grito provoque como efecto una avalancha de mu- 
chas toneladas de nieve, evidentemente desproporcionado a 
su causa. E E 

Si se trata, en cambio, de sistemas dinámicos que repiten pe- 
riódicamente un comportamiento, en aquellos casos en que 
un pequeño cambio en las condiciones iniciales se reproduzca 
multiplicado en cada uno de los periodos siguientes, de modo 
que haya una situación llamada de realimentación positiva, 
y que se describe mediante ecuaciones no lineales, se puede 
llegar a situaciones tan diferentes que no haya posibilidad de 
predecirlas. 

Consideremos un sistema dinámico, formado por una esfera 
sometida a una fuerza que la haga rodar sobre un riel circular 
que tiene una longitud de 10 metros, como en la figura 11.2. 


Referencia 


FIGURA 11.2. 


Para ubicar la posición de la esfera en cualquier momento 
mediremos la longitud de riel que la separa de un punto fijo 
de referencia que hemos marcado sobre el mismo. Suponga- 
mos que la posición inicial de la esfera es de 1 metro respecto 
de esa referencia, y que medimos esa posición con una cinta 
métrica que nos da un error de +1 milímetro, es decir, que la 
posición real puede ser de 0.999 m, o de 1.001 m o cualquier 
valor comprendido entre esas cifras. 

Supongamos además que cada vez que la esfera pasa por el 
punto de referencia recibe un impulso, y que éste provoca un 
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desplazamiento de la esfera sobre su órbita que incrementa en 
un 10% la posición, de modo que si estaba entre 0.999 y 1.001 m 
en el momento inicial, se multiplica por un factor 1.1 y dicha 
posición pasa a ser entonces de entre 1.099 y 1.101 m al cabo 
del primer periodo, y a su vez vuelve a multiplicarse por 1.1 
para cada uno de los periodos sucesivos por lo que la posición 
será de: entre 1.209 y 1.211 m en el segundo periodo, entre 
1.329 y 1.332 m en el tercero, etcétera. Al cabo de 25 vueltas, la 
posición será prácticamente la misma que al principio: entre 
10.825 y 10.845 m, que como la longitud total del riel es de 10 
metros, ubica a la esfera entre 0.825 y 0.845 m respecto de la 
referencia. Pero obsérvese que si bien la esfera ha vuelto en el 
periodo 25 a estar prácticamente en la posición inicial, la inde- 
finición en su posición ha crecido: es ahora una zona de +10 mm 
de longitud sobre el riel. La imprecisión inicial de +1 mm en 
su ubicación se ha multiplicado por 10. 

En la figura 11.2 se ha ilustrado la situación a las 25 vueltas y 
después de 60 vueltas, donde la esfera puede estar en cual- 
quier posición entre 3.81 y 4.19 metros, y a las 70 vueltas don- 
de estará entre 8.21 y 9.79 metros. Después de 97 vueltas, esta 
indefinición en la posición se habrá amplificado 10 mil veces y, 
por lo tanto, será de 10 metros. De modo que la esfera puede 
estar en cualquier punto del riel de 10 metros de longitud, y 
no podemos saber de antemano cuál es ése punto mediante 
un cálculo. 

Por supuesto, esto parece fácilmente remediable: ya que se 
ha utilizado para medir la posición inicial una cinta métrica 
que da un error de 1 mm, se la puede cambiar por un instru- 
mento mucho más preciso, que dé un error 10 mil veces me- 
nor, de solamente 1 décimo de micrón. 

Ahora sí, con este dato inicial tan preciso se puede ubicar la 
esfera mediante un cálculo que predecirá su posición dentro 
de 97 vueltas con un error no mayor que 1 mm. 

Pero, ¿y si se pretende seguir prediciendo mediante el cálcu- 
lo la ubicación para un mayor número de vueltas? 

Ocurre que aun con esta medición inicial mucho más pre- 
cisa, a las 193 vueltas se está otra vez en la situación primera: 
una indefinición de 10 metros, que no permite predecir dón- 
de está la esfera. Como es obvio no se puede seguir aumentan- 
do de manera indefinida la precisión de la medida, pues rápi- 
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damente aparece la necesidad de medir magnitudes con un 
error menor que el tamaño de un átomo, y aun cuando se pu- 
diese realizar semejante maravilla, unos cientos de vueltas más 
llevarían a la misma frustrante situación. 

Estamos aquí entre la espada y la pared: o utilizamos una 
medición aproximada, lo que no es lo suficientemente exacto 
como para hacer la predicción para cualquier estado futuro, o 
intentamos realizar una medición tan extremadamente precisa 
que es impracticable. 

¿Un sistema dinámico que tiene este comportamiento no es 
raro en la naturaleza. Se trata de sistemas que repiten perió- 
dicamente un cierto estado, y donde éste sufre una perturba- 
ción que se multiplica en cada periodo (es el caso de la reso- 
nancia), lo que es suficiente para que se produzcan después de 
cierto tiempo efectos de una magnitud desproporcionadamen- 
te grande en comparación con la perturbación inicial. En es- 
tos sistemas, regidos por leyes estrictamente deterministas, en 
ciertas condiciones el cálculo del comportamiento puede re- 
sultar crecientemente impreciso, haciendo que todo intento 
de predecir el futuro a largo plazo sea derrotado. 


Los PROCESOS FÍSICOS ESTÁN DETERMINADOS, 
PERO ¿SE PUEDEN PREDECIR SIEMPRE? 


“Se hace necesario ahora distinguir la diferencia entre deter- 


minismo y predecibilidad, dos palabras que desde la época de 
Laplace se consideraban sinónimas. Se ve que, en ciertas con- 
diciones, un sistema dinámico puede ser determinista y, sin 
embargo, tener un comportamiento no predecible. 

Por supuesto que los físicos sabían que el Universo muestra 
una inmensa complejidad, pero se suponía que, en general, se 
puede analizar la misma para descomponerla en sus compo- 
nentes simples, y que de los mismos se pueden seleccionar unas 
pocas variables como condiciones iniciales para introducir en 
las ecuaciones, ya que la inmensa mayoría de las otras condi- 
ciones es de tan pequeña magnitud, que el no considerarlas 
no afecta el resultado. Pero aprendemos ahora que estos casos 
son algo especial, y que mucho más común es la situación en 
que aun cuando se trate de sistemas para los que rige el deter- 
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minismo, hay una limitación de la posibilidad de predecir el 
comportamiento futuro que puede llegar a ser prácticamente 
equivalente a una situación de caos. 

No significa, sin embargo, que no se pueda decir nada sobre 
dichos sistemas, ya que para ciertas condiciones iniciales, el 
comportamiento si es ordenado y, por lo tanto, predecible a 
largo plazo. Además, aun en el estado de comportamiento caó- 
tico, estos sistemas presentan muchas propiedades que se pue- 
den comprender con la ayuda de la teoría de probabilidades, 
siendo esta mezcla de determinismo y probabilidad una forma 
muy fructífera de atacar los problemas característicos de los 
sistemas complejos. 

Se trata entonces de una nueva física de fenómenos no li- 
neales cuyo objeto es el estudio de procesos tales como, por 
ejemplo, las turbulencias y oscilaciones físico-químicas y bioló- 
gicas. Estos fenómenos tienen un aspecto aparentemente azaro- 
so, pero con similitudes insospechadas en su comportamiento, 
lo que conduce a que se encare su estudio mediante métodos 
matemáticos, donde resultan comunes nociones como la bifur- 
cación, los atractores extraños y los exponentes de Lyapunov. 

Este nuevo enfoque de los fenómenos complejos implica un 
cambio de paradigma, ya que el planteo de Laplace con su de- 
monio capaz de predecir cualquier evento queda puesto en 
sus justos términos: no se trata de predecir absolutamente to- 
dos los fenómenos que aparecen en la naturaleza merced a los 
cálculos basados en las leyes deterministas; ni Laplace ni nin- 
gún físico han soñado que esto sea posible. Pero el demonio 
pone de relieve que esta imposibilidad se debe a nuestra im- 
perfección, como seres limitados que somos. 

Sólo un Ser absoluto podría conocer el futuro a cualquier 
plazo, siempre y cuando conociera todas las condiciones ini- 
ciales con infinita precisión, y, en sus cálculos, además mane- 
jara números con una infinita cantidad de cifras, como los 
números irracionales. Ni los seres humanos ni las computado- 
ras pueden hacer esto, ya que necesariamente utilizan una 
cantidad finita de condiciones iniciales y un número limitado 
de cifras en sus cálculos. 

Éste es un enfoque platónico, que postula que existe una 
realidad matemática que sí es perfecta, un ideal que marca la 
meta hacia donde se deben dirigir nuestros esfuerzos. La ac- 
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ción coherente con este enfoque consiste en concentrar los es- 
fuerzos en hallar la solución de las ecuaciones que describen 
al fenómeno, para introducir en ellas los datos empíricos ob- 
tenidos cada vez con mayor precisión, en la confianza de que 
con el transcurso del tiempo se estará así más y más cerca de 
predecir un creciente repertorio de comportamientos. 

Es esta visión de continuo progreso la que se cuestiona ac- 
tualmente para el estudio de los sistemas complejos, ya que la 
sensibilidad a las condiciones iniciales puede ser tal que ningu- 
na observación, no importa cuán precisa sea, nos llevará a de- 
terminar esas condiciones con la exactitud suficiente. Si bien 
en un sistema determinista el estado inicial determina el esta- 
do final, en estos sistemas no lineales el conocimiento. aproxi- 
mado del estado inicial no permite deducir un conocimiento 
aproximado del estado final, por lo que pierde sentido poner 
la prioridad en resolver las ecuaciones para introducir datos 
cada vez más exactos en ellas. 

Esta dificultad para las predicciones en sistemas determi- 
nistas complejos ya se conocía en el siglo pasado: el físico Ja- 
mes C. Maxwell, que tuvo un papel fundamental en la teoría 
del electromagnetismo, expresó en una de sus conferencias en 
1873, que si bien las mismas causas producen los mismos efec- 
tos, cuando hay sensibilidad a las condiciones iniciales, las cau- 
sas semejantes pueden no tener efectos semejantes. Da entre 
otros ejemplos la explosión del algodón pólvora y el efecto de 
pequeñas acciones humanas que pueden desencadenar gran- 
des cambios sociales. 

El gran matemático Henri Poincaré escribió en 1908, en su 
Ciencia y método, que “una causa muy pequeña, que se nos es- 
capa, determina un efecto considerable que no podemos dejar 
de ver y entonces decimos que ese efecto se debe al azar”. Lue- 
go se extiende en este concepto: 


¿Por qué los meteorólogos encuentran tan difícil prever el tiempo 
con alguna certeza? ¿Por qué las lluvias, las tempestades mismas 
nos parecen llegadas al azar, de modo que mucha gente cree natu- 
ral rezar para que caiga la lluvia o haga buen tiempo, cuando en 
realidad encontrarían ridículo pedir con una plegaria un eclipse? 
Vemos que las grandes perturbaciones se producen en general en 
las regiones donde la atmósfera está en equilibrio inestable. Los 
meteorólogos bien ven que ese equilibrio es inestable, que algún 
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ciclón se dará en alguna parte, pero ¿dónde? No están en condi- 
ciones de decirlo; una décima de grado más o menos en un pun- 
to cualquiera y el ciclón estalla aquí y no allá y extiende sus estra- 
gos en comarcas que de otra manera no habrían sido devastadas. 
Si se hubiera conocido esa décima de grado podría habérselo 
sabido de antemano, pero las observaciones no fueron ni bastante 
rigurosas ni bastante precisas y por eso todo parece debido a la 
intervención del azar. 


Para Poincaré esa aparición del azar en la naturaleza tiene 
dos raíces principales: 1) la ya mencionada sensibilidad a las 
condiciones iniciales, que hace que aunque las leyes del sis- 
tema se conozcan, los pequeños errores iniciales tienen una 
enorme influencia sobre el estado final, y 2) la complejidad de 
las causas, por lo que cuando en el proceso de análisis del sis- 
tema se le aísla idealmente, se toma en cuenta apenas una 
parte de todas las innumerables influencias a las cuales está 
sometido. 

Este procedimiento es efectivo sólo para los sistemas integra- 
bles, que, en general, son sistemas lineales. En cambio para los 
no lineales, esta separación para su análisis se hace pagando el 
precio de disminuir el alcance de la predicción. 

Tanto Maxwell como Poincaré pusieron como ejemplo de 
los sistemas dinámicos con sensibilidad a las condiciones ini- 
ciales, el caso de un gas formado por muchas moléculas que 
vuelan a gran velocidad en todas direcciones y tienen muchas 
colisiones entre sí. Si se desea determinar la evolución de un 
sistema de una cantidad de componentes tan inmensa como 
los que forman un gas, se debe tener en cuenta que es imposi- 
ble calcular con las ecuaciones diferenciales de Newton el 
movimiento de cada una de las partículas. 

Piénsese que un centímetro cúbico de aire contiene 27 tri- 
llones (millones de millones de millones) de átomos, y que no 
alcanzaría toda la superficie de la Tierra para escribir nada más 
que las ecuaciones diferenciales correspondientes, sin men- 
cionar el cálculo para determinar con ellas el movimiento de 
cada átomo en ese centímetro cúbico de gas. 

Supongamos que pudiésemos observar una pequeñísima 
porción de ese centímetro cúbico de aire, una porción lo sufi- 
cientemente limitada como para poder registrar el movimien- 
to de las pocas moléculas que estuviesen allí presentes. Pode- 
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mos imaginar que al aislar así idealmente nuestra zona de ob- 
servación, podríamos seguir las trayectorias en línea recta indi- 
viduales, las colisiones entre moléculas y cómo rebotan, y com- 
prenderíamos que éstas siguen las leyes de la mecánica y que 
son perfectamente predecibles, ya que podríamos calcular sus 
movimientos individuales... ¡Inútil esfuerzo, ya que sólo se- 
ría válido por un brevísimo tiempo! Ocurre que repentina- 
mente llegarían desde fuera de nuestra zona de observación 
más moléculas, que no figuraban entre nuestros datos inicia- 
les, y desaparecerían de la zona moléculas que estaban en 
observación, y el resultado sería ahora un movimiento distinto, 
con nuevas colisiones imprevistas, imposible por lo tanto de 
calcular, tan imprevisible como si se tratara de un proceso regi- 
do por el azar. Comprendemos así que si se observa sólo una 
fracción ínfima de un proceso muy complejo, formado por 
innumerables componentes, con una enorme cantidad de dis- 
tintas variables que actúan sobre el sistema, éste puede apare- 
cer como aleatorio, desordenado. 

Esto es válido no solamente para un conjunto de moléculas 
en un gas: hay cierto paralelismo con los fenómenos que estu- 
dian las ciencias sociales, ya que si se desea examinar la eco- 
nomía de una nación o su evolución política, estudiando una 
pequeña parte aislada del conjunto, se va hacia un fracaso, 
pues esa parte estará sujeta continuamente a las influencias in- 
controlables e imprevistas del resto del sistema. 

¿Cómo proceder entonces? Parecería que al estar vedado el 
conocimiento detallado de cada comportamiento individual, 
habría que resignarse a aceptar que los sistemas sociales, bioló- 

- gicos, económicos y físicos con una gran cantidad de compo- 
nentes no se pueden encarar con el método científico, ya que 
tienen un comportamiento que aparentemente depende del : 
puro azar, aunque se pueda pensar que están regidos por leyes 
deterministas. 

La respuesta a esta pregunta surgió en el campo de las cien- 
cias sociales del siglo XIX, cuando se quisieron estudiar las ca- 
racterísticas sanitarias y económicas de las naciones..Se aplicó 
para ello la rama práctica de las probabilidades, que es la esta- 
dística. 

Su desarrollo matemático fue iniciado por A. Quetelet en 
1820, con su libro Mecánica social, inspirado en el modelo de la 
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Mecánica celeste de Laplace, y que tenía como objetivo desarro- 
llar las “ciencias morales y políticas”. Aparece allí la idea del 
hombre promedio, con una longitud de vida, ingreso econó- 
mico, etcétera, que se calculan como promedio de los de toda 
la población. 

La experiencia con el estudio de los gases enseña que, a pesar 
de que a escala microscópica hay un número de componentes 
tan grande que es imposible conocer sus comportamientos in- 
dividuales, éstos muestran en conjunto a escala macroscópica 
propiedades globales como la temperatura, presión, densidad, 
volumen, que se pueden medir, y que se relacionan entre sí 
con leyes bien definidas. Esto es análogo a lo que acontece con 
una población humana de un número grande de individuos, 
para la que se pueden establecer parámetros como la tasa de 
nacimientos, muertes, producto per cápita, etcétera, sin nece- 
sidad de seguir cada historia individual. 

El método estadístico dio tan buenos resultados para las cien- 

cias sociales que luego lo aplicaron los físicos, creando así la 
mecánica estadística, una rama de la física desarrollada alrede- 
dor de 1900 por el austriaco Ludwig Boltzmann y el nortea- 
mericano J. Willard Gibbs. 
Es así que James Clerk Maxwell propuso en 1872 la intro- 
ducción en la física del cálculo de probabilidades para sis- 
temas de gran número de componentes, y escribió: “hemos 
hallado un nuevo tipo de regularidad, la regularidad de los 
promedios”. 


Los JUEGOS DE AZAR Y LAS PREDICCIONES 


La teoría de probabilidades se originó en los juegos de azar: 
Laplace establece en su libro que: 


La probabilidad de un evento es la cantidad de veces en que apa- 
rece, dividido por el total de los que podrían ocurrir suponiendo 


que todos ellos sean igualmente posibles. 


Así por ejemplo, un dado de seis caras nos da una probabili- 
dad de 1/6 para que caiga con la cara del 2 hacia arriba, y la 
misma probabilidad vale para las otras cinco caras. 
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Si se arroja un dado un gran número de veces y se hace la 
estadística de cuántas veces aparece cada uno de sus seis nú- 
meros, se comprobará que cuanto mayor sea el número de 
tiradas tanto más cerca se estará del número 1/6. Pero no se 
sabe antes de arrojar el dado qué cara quedará esta vez arriba y 
ésta es, por supuesto, la razón de la existencia de los juegos de 
azar. 

En teoría, si se conociesen las fuerzas aplicadas al lanzar el 
dado, los rozamientos, resistencia del aire, etcétera, se podría 
predecir cómo caerá, pero es imposible disponer de todos los 
datos necesarios, y, por lo tanto, los jugadores podrán seguir 
disfrutando de la emoción de apostar sin saber qué les depa- 
rará el azar. E 

De modo que cuando se trata de sistemas dinámicos en los 
que no se pueden aplicar las ecuaciones de Newton para pre- 
cisar el comportamiento de cada caso individual, lo más cer- 
cano a una predicción es establecer la probabilidad de que 
ocurra un evento determinado. 

Se debe responder ahora a un nuevo tipo de pregunta. No 
ya: ¿cuál es la velocidad de esta partícula?, sino: ¿qué probabili- 
dad tengo de encontrar partículas cuya velocidad esté entre tal 
y cual valor? 

El método se aplicó en la física inicialmente en el estudio de 
los gases, lo que es lógico por la analogía a la que nos hemos 
referido entre la relación de las moléculas con el gas y la de los 
individuos con una sociedad. Este concepto fue introducido 
en la termodinámica, que surgía como el desarrollo científico 
de los estudios sobre el calor. 

A fines del siglo xIx existían entonces dos enfoques científi- 
cos para la formulación matemática de los fenómenos natura- 
les: el análisis mediante las ecuaciones diferenciales para sis- 
temas con pocos componentes, y el análisis estadístico para los 
sistemas con un número elevado de componentes. 

Pero en la comunidad científica, la distinción jerárquica en- 
tre ambos enfoques era bien clara: la jerarquía máxima corres- 
pondía a la primera de estas visiones, ya que expresaba el en- 
foque que hemos llamado platónico, de que la naturaleza está 
regida por leyes eternas escritas en lenguaje matemático, mien- 
tras que la segunda debe su existencia a la ignorancia que sur- 
ge de nuestras limitaciones humanas, la que se compensa par- 
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cialmente mediante el uso de las probabilidades, pero que 
inevitablemente irá viendo achicar su campo de aplicación a 
medida que avance el conocimiento, gracias a ese progreso 
continuo e inevitable que fue una creencia básica en Occi- 
dente a partir del Iluminismo. 

Como vimos, tal visión cambió radicalmente en el transcur- 
so de este siglo, siendo Henri Poincaré el primero que perci- 
bió las características esenciales del nuevo enfoque. Al igual 
que Maxwell, advirtió que sistemas simples, de pocos compo- 
nentes y que, por lo tanto, están en la primera de las catego- 
rías que hemos descrito, pueden comportarse no obstante caó- 
ticamente y requerir entonces de métodos estadísticos para su 
descripción. 

Poincaré advirtió esto al examinar el caso de un sistema di- 
námico muy simple, formado por sólo tres cuerpos que se 
atraen por acción de la fuerza de la gravedad, cuando uno de 
ellos es muy pequeño con relación a los otros dos. Trató de de- 
terminar cómo se movería y de representar gráficamente su 
evolución como una trayectoria en un espacio matemático lla- 
mado “espacio de las fases” (sobre este espacio volveremos más 
adelante). La estructura de esta trayectoria, llamada homoclí- 
nica, resultó ser tan extraordinariamente complicada y alejada 
de la intuición que Poincaré desistió de intentar dibujarla y 
escribió en sus Nuevos métodos de mecánica celeste: 


Intentemos hacernos una idea de la figura formada por estas dos 
curvas y sus infinitas intersecciones, cada una de las cuales corres- 
ponde a una solución doblemente asintótica. Dichas interseccio- 
nes forman una especie de entramado, de tejido, de red de mallas 
infinitamente finas. Ninguna de esas curvas se ha de cortar a sí 
misma; además, se ha de plegar sobre sí misma de un modo muy 
complejo para cortar infinitas veces todas las mallas de la red. 

Sorprende la complejidad de esta figura que ni siquiera intento 
dibujar. Nada más indicado para poder hacernos una idea de la 
complejidad del problema de tres cuerpos y, en general, de todos 
los problemas de la Dinámica en los que no hay una integral uni- 
forme. 


Esta inesperada complejidad, que asombró tanto al genial 
matemático, es la que hoy se estudia en los sistemas dinámicos 
con estructuras como las homoclínicas, las heteroclínicas y 
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también en los atractores extraños, representación fundamen- 
tal para el examen del caos determinista. 

La obra de Poincaré ha sido básica para el estudio de los sis- 
temas dinámicos y lo consagra como el fundador de los méto- 
dos que hoy permiten abordar las matemáticas del caos de- 
terminista. Sin embargo, hubo un periodo de más de sesenta 
años en los que sus aportes para funciones no lineales fueron 
prácticamente ignorados, y es recién a partir de los trabajos de 
A. Kolmogorov y V. Arnold en Rusia, S. Smale en Estados Uni- 
dos, D. Ruelle y R. Thom en Francia, que se vuelve a trabajar 
seriamente en este campo. 

A pesar de que Poincaré fue un visionario, no avanzó en el 
estudio de estos sistemas, lo que se puede atribuir a dos moti- 
vos principales. El primero es que en su época los cálculos nu- 
méricos se hacían manualmente, lo que imposibilitaba el tra- 
tamiento numérico de los sistemas de ecuaciones no lineales, 
que debieron esperar a la aparición de las computadoras. El 
otro motivo es su actitud filosófica, que compartía con la ma- 
yoría de los matemáticos y físicos de la época, y que se refleja 
claramente en un escrito de Poincaré donde comenta la exis- 
tencia de ciertas funciones matemáticas que no responden a 
las exigencias clásicas de ser continuas y derivables: “La lógica 
a veces engendra monstruos. Desde hace medio siglo se han 
visto surgir una multitud de funciones raras que parecen esfor- 
zarse en asemejarse lo menos posible a las honestas funciones 


_ que sirven para cualquier cosa. No más continuidad, no más 


derivadas, etcétera.” 

¡Cuál sería el asombro de Poincaré ante el atractor que 
Lorenz descubrió al modelizar fenómenos meteorológicos! No 
por el aspecto matemático, ya que, como vimos, comprendía 
perfectamente la imposibilidad de hacer previsiones a largo 
plazo, sino porque ese atractor extraño, que permite graficar 
las características esenciales del sistema dinámico de la atmósfe- 
ra, se engendró utilizando los monstruos matemáticos “raros” 
de los que tanto desconfiaba. El había llegado a vislumbrar el 
universo fascinante de las matemáticas de lo imprevisible, pero 
retrocedió ante esa falta de regularidad, de continuidad, que 
le producía tanto desagrado. 

Una consecuencia de esta visión fue que los científicos im- 
buidos de ella han tendido a ignorar ciertos campos de la in- 
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vestigación, como el que alberga a los “monstruos” que mencio- 
naba Poincaré (por ejemplo, las curvas sin tangentes y todos 
los algoritmos que hoy llamaríamos “caóticos”). 

Para poder considerar de real valor científico al caos deter- 
minista, hay que renunciar a la creencia en una jerarquía don- 
de la cima está ocupada por las formas perfectas (círculo, 
esfera, etcétera). Pues, ¿por qué aceptar que un círculo es su- 
perior a una figura compleja, como lo es un fractal o cual- 
quier elemento de la naturaleza? Como dijera el matemático 
B. Mandelbrot, “las montañas no son conos, ni las nubes son 
esferas, ni tampoco los ríos ni los rayos van en línea recta; en 
definitiva, las formas que nos rodean a lo largo de nuestra 
vida, siempre han estado muy alejadas de la simplicidad de lo 
euclidiano”. 

Este cambio de valoración tiene un origen cultural y filosófi- 
co, y si se pasa esto por alto, se arriesga reducir el estudio mo- 
derno del caos al de un mero conjunto de nuevas técnicas. 

Un aspecto importante de este nuevo enfoque es el indis- 
pensable uso de la computadora, que permite hacer los cálcu- 
los y además traducir en imágenes el comportamiento de los 
sistemas dinámicos. 

También en este aspecto hay una clara distinción entre las 
actitudes de los matemáticos que adhieren a uno u otro en- 
foque. Así como la utilización de herramientas que no fueran 
la regla y el compás era rechazada por los geómetras griegos, 
hoy continúa la resistencia de muchos matemáticos “puros” al 
uso de las computadoras para las demostraciones matemáticas, 
basados en el hecho cierto de que el número de decimales de 
todo resultado numérico que puede manejar una computado- 
ra es limitado. 

Otro aspecto nuevo es la introducción deliberada de lo alea- 
torio en muchos de los cálculos, por ejemplo, mediante series 
de números generados al azar, con el objeto de modelar fenó- 
menos naturales que tengan componentes imprevisibles. 

El uso de las probabilidades había sido iniciado, como vi- 
mos, en el pasado por Quetelet en las ciencias sociales y por 
Maxwell en la teoría cinética de los gases, pero se lo considera- 
ba como un índice de ignorancia, que señalaba los límites de 
las zonas donde están los temas que la verdadera ciencia había 
logrado esclarecer; y en esas zonas lo aleatorio, las probabili- 
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dades sólo tenían cabida mientras no se pudiera aplicar la 
verdadera ciencia. 

Pero en el primer cuarto de este siglo, la mecánica cuántica 
mostró que el concepto de probabilidad y los métodos estadís- 
ticos son esenciales para la formulación de las leyes físicas, al 
menos en la escala atómica, con lo cual se comenzó a apreciar 
su significado conceptual. 

Los matemáticos y físicos tenían presente. que hay cierta cuo- 
ta de desorden en la naturaleza, pero no creían que su estudio 
fuera digno de constituir una verdadera ciencia. 

Ahora, en cambio, el caos está plenamente reconocido y se 
acepta que además de los fenómenos que exhiben orden y re- 
gularidad, hay un vasto universo de fenómenos desordenados, 
irregulares, que no son reducibles a las formas puras, y que 
pueden también aparecer insospechadamente en sistemas muy 
simples, como el formado por tres cuerpos. 

Los fenómenos caóticos pueden ser estudiados científica- 
mente, pues aunque no hay posibilidad alguna de prever el 
comportamiento detallado de sus componentes individuales, 
sí se puede hacer una predicción cualitativa de la evolución 
del sistema en su conjunto, y también se pueden buscar las 
condiciones para que un sistema dinámico que está en estado 
de orden pase al comportamiento caótico y viceversa. 

Como se puede ver, algunas de las ideas más interesantes de 
las ciencias actuales están referidas a estos temas, los que han 
obligado a precisar los verdaderos alcances del determinismo, 
de nuestra capacidad de predecir el futuro y del verdadero sig- 
nificado de las leyes naturales. Sólo el tiempo nos dará una 
idea exacta de la verdadera trascendencia de estos cambios 
sobre el método y la filosofía de la ciencia. 

Pasaremos ahora a examinar las características de los sis- 
temas dinámicos que pueden presentar un comportamiento 
caótico. 


EN LA NATURALEZA HAY MUCHOS PROCESOS CÍCLICOS 
Con respecto a los procesos dinámicos que se observan en la na- 
turaleza, cuando se les quiere estudiar se debe subrayar la uni- 


versalidad de los que son periódicos o de oscilación, procesos 
ubicuos en la física, la astronomía y la biología, ya sea en la 
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estructura de las nubes en la atmósfera, los remolinos en un 
río torrentoso, el sonido melodioso de un violín o el flujo de 
energía eléctrica en un circuito electrónico, y van desde el co- 
nocido movimiento de los péndulos y de las órbitas de los pla- 
netas a los complejos ritmos biológicos: fluctuaciones en los 
ciclos de reproducción de una población animal, respiración, 
ritmo cardiaco, alternancia de vigilia y sueño, procesos neu- 
rofisiológicos, etcétera. En todos ellos se halla esta característi- 
ca de periodicidad en el comportamiento. 

Como son legión los dispositivos y fenómenos que presen- 
tan este comportamiento oscilatorio, es natural que el oscila- 
dor suscite un creciente interés científico. Pero ahora vamos 
mucho mas allá, pues la conclusión que se extrae de la visión 
anterior es que las propiedades de los osciladores pueden apli- 
carse al conjunto de los procesos dinámicos. 

En efecto, se puede describir matemáticamente cualquier 
proceso que fluctúa en el tiempo como la suma algebraica de 


FIGURA IL. 
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un conjunto de oscilaciones periódicas, por un método llama- 
do análisis de Fourier. 

Como se puede ver en la figura IL.3, si se grafica la variación 
en el tiempo del proceso O que queremos analizar, se obtiene 
una curva irregular, que no se repite en el tiempo. Pero Fourier 
demostró que se le puede considerar como la resultante de su- 
mar varias curvas periódicas regulares, como las 4, By C, y, por 
lo tanto, toda teoría de procesos dinámicos puede requerir la 
aplicación del concepto de los osciladores, que combinados 
dan cuenta de cualquier variación en el tiempo. 

Es por este motivo que, como iniciación al estudio de los sis- 
temas dinámicos examinaremos en detalle el comportamiento 
de un oscilador elemental, tal como el péndulo simple, y luego 
veremos qué ocurre cuando se combinan dos o más de estos 
osciladores elementales para formar un sistema dinámico. 


III. Los péndulos simples 
pueden ser muy complejos 


ANALIZAREMOS aquí en qué condiciones será caótico el movi- 


miento de un sistema tan poco complicado como puede serlo 
el formado por un par de péndulos. 

El péndulo es un sistema dinámico que repite su compor- 
tamiento a intervalos regulares que llamamos periodos. 

Ha sido estudiado por lo menos desde la época de Galileo, 
de modo que se podría suponer que ya no tendría por qué 
proporcionarnos ninguna sorpresa. Por otra parte, desde hace 
siglos el péndulo constituye el paradigma de la previsibilidad y 
de la regularidad. Antes de la utilización de los osciladores de 
cuarzo, los relojes eran regulados por péndulos. 


SISTEMAS DINÁMICOS CONSERVATIVOS O DISIPATIVOS 


Se debe ante todo distinguir entre sistemas dinámicos conserva- 
tivos, que son aquellos en los que la energía del sistema se man- 
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tiene constante porque hay ausencia de frotamiento o fricción 
interna, y aquellos en los que debido al frotamiento hay una 
disminución continua de energía, llamados sistemas disipativos. 

A los primeros se los suele denominar también “sistemas ha- 
miltonianos”, ya que resulta muy fructífero describir su com- 
portamiento con una función matemática desarrollada por 
W. Hamilton que utiliza como variables las posiciones y, en 
lugar de las velocidades, los impulsos (producto de velocidad 
por masa). 

Ejemplos de sistemas hamiltonianos son el sistema solar y el 
plasma en un acelerador de partículas, y de los disipativos la 
atmósfera terrestre, los océanos, los organismos vivientes y to- 
das las máquinas. 

Examinemos las características de los sistemas hamiltonia- 
nos, comenzando con el caso más sencillo, o sea el del péndu- 
lo ideal, que nadie hallará en ningún laboratorio, pero que 
permite fijar las ideas esenciales para estudiar los osciladores. 

Este péndulo ideal se mueve en un espacio de dos dimen- 
siones, o sea, sobre un plano vertical, el hilo del que cuelga es 
rígido y no tiene peso, no fricciona en el pivote y está en el 
vacío, por lo tanto no hay resistencia del aire a su movimiento. 

Se trata entonces de un sistema hamiltoniano simple, y es 
posible describir matemáticamente su movimiento mediante 
la ecuación de Newton: una ecuación diferencial que vincula la 
derivada segunda del ángulo que forma el hilo con la vertical, 
respecto del tiempo (o sea la variación en el tiempo de la 
velocidad con que cambia el ángulo), con otros parámetros 
como la longitud del hilo y la aceleración de la gravedad. 

Si se quiere representar gráficamente su movimiento, nece- 
sitaremos calcular para cada momento su ángulo respecto de 
la vertical y su velocidad, que llamaremos 0 (theta) y v respecti- 
vamente. Deseamos saber cómo varían ambos a medida que 
transcurre el tiempo y para ello nos será de suma utilidad el 
uso del espacio de las fases. 


EL ESPACIO DE LAS FASES Y LOS PÉNDULOS ACOPLADOS 


Este es un espacio matemático abstracto, que no se debe con- 
fundir con aquel en el que realmente se. mueven los compo- 
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nentes del sistema, pero que contiene en sus formas geométri- 
cas una información concreta: las variables que describen el 
movimiento del sistema dinámico. 

Se utilizará aquí un método que es exactamente el opuesto 
al que utilizó René Descartes cuando concibió el uso de coor- 
denadas. 

Descartes descubrió cómo transformar la geometría en nú- 
meros al imaginar el espacio como una inmensa grilla, el sis- 
tema de coordenadas cartesianas, de forma que la posición de 
cualquier punto en el espacio esté definida mediante números 
que miden su distancia a estas coordenadas de referencia. Por 
lo tanto, toda forma geométrica se puede expresar en números 
mediante una ecuación matemática que la refiera a las coorde- 
nadas cartesianas. 

Al introducir ahora el espacio de las fases hacemos al revés: 
transformamos números en formas geométricas pretendiendo 
que estos números sean coordenadas en dicho espacio imagl- 
nario. 

La ventaja de proceder así es que basta con observar esa for- 
ma geométrica para saber cómo varía el comportamiento del 
sistema a medida que transcurre el tiempo. Para ello se utili- 
zan como ejes de coordenadas las variables dinámicas del sis- 
tema, o sea aquellas magnitudes que van cambiando con el 
tiempo, como pueden ser la velocidad del péndulo y su ángulo 
respecto de la vertical. Para un momento inicial, el péndulo, 
entonces, queda representado por un punto en el gráfico del 
espacio de las fases que nos indica su velocidad y ángulo. Un 
tiempo después tendrá otra posición y velocidad, y a esto le 
corresponde un punto distinto, y así podremos hacer la histo- 
ria del péndulo mediante puntos sucesivos que trazan una tra- 
yectoria. De esta forma, la dinámica cambia de estudiar largas 
listas de números a visualizar cómo va evolucionando con el 
tiempo una figura geométrica. 

Si el resultado después de un tiempo es un único punto fijo, 
el sistema está estático, no evoluciona. Si hay una curva y ésta 
es cerrada, esto indica que el sistema repite periódicamente su 
comportamiento; si es una curva abierta, se deberá examinar 
sus características para ver si hay o no regularidades subyacentes. 

Fue el gran matemático Henri Poincaré quien tuvo la bri- 
llante idea de proponer este método, por el que la dinámica se 
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puede hacer visual. Los investigadores intentan descubrir las 
características generales de un sistema, que se aprecian mejor 
si se observan las formas que aparecen en el espacio de las fases. 

En general, esta manera de representar un sistema dinámico 
tiene entre sus ventajas el que las coordenadas del espacio de 
las fases pueden representar cualquier característica del sis- 
tema dinámico que varíe en el tiempo, como, por ejemplo, las 
señales eléctricas del corazón, la población de abejas en una 
colmena, el valor del dólar, etcétera. 

Un principio general de las trayectorias en el espacio de las 
fases es que ninguna puede tocar a otra como consecuencia 
del carácter determinista de esta descripción: si se intersecta- 
ran dos trayectorias en un punto, a partir del mismo habría 
dos curvas diferentes, que corresponderían a dos soluciones 
distintas de las ecuaciones diferenciales del sistema, o sea dos 
comportamientos distintos al mismo tiempo. 

Así, si estuviésemos representando como sistema dinámico a 
un automóvil, si dos trayectorias se tocasen en un punto indi- 
carían que tenemos a un mismo vehículo moviéndose simul- 
táneamente a, por ejemplo, 30 y 90 km por hora, lo que ob- 
viamente, se opone al carácter no ambiguo, determinista, de 
estos fenómenos. 

Para representar el comportamiento del péndulo ideal en el 
espacio de las fases, comencemos por dibujar dos ejes de refe- 
rencia, o coordenadas, uno horizontal O para los ángulos, y 
otro vertical v para las velocidades. Marquemos ahora el cami- 
no o trayectoria del péndulo en este mapa como en la figura 
11.10 de la página siguiente. Supongamos que a partir de un 
instante inicial en que le damos un impulso, medimos cada 
décimo de segundo el ángulo y velocidad, y marcamos el pun- 
to correspondiente en la figura 111.10. Como podemos ver en 
la figura Ml. 1 a, en el tiempo inicial 7'= 0, el péndulo forma un 
ángulo —0 respecto de la vertical, situación que se representa 
con el punto T'=0 en la figura MI.10. Al liberarlo se moverá 
hacia la derecha con una velocidad creciente, que será la máxi- 
ma cuando pase por la vertical en el instante 7'= 1, momento 
en el que O = 0, y esto se representa en la figura 111.15 con el 
punto T'= 1. Al llegar al máximo desplazamiento a la derecha, 
para 0 = +6, su velocidad se habrá reducido a cero, y esto co- 
rresponde al punto T= 2 en la figura ML. 1b. 


45 


FIGURA HI 1. 


Continuando este procedimiento se obtendrá un conjunto 
de puntos, que trazan la trayectoria dinámica del oscilador, tra- 
yectoria que representa el movimiento entero del péndulo en 
un ciclo. 

Dado que el camino se repite ciclo tras ciclo, el mapa para 
este péndulo simple es una única trayectoria cerrada, también 
llamada órbita, por analog con el movimiento de los pla- 
netas. 

Si le damos al péndulo un mayor impulso inicial, el ángulo 
máximo será mayor. De hecho, en un único gráfico podemos 
representar el movimiento del mismo péndulo con diferentes 
velocidades iniciales, y para cada una de ellas se tiene una órbi- 
ta distinta. Esta es una de las características de los sistemas con- 
servativos o hamiltonianos. 

Se obtiene así una familia de curvas, que puedan cubrir el 
plano 0,v (figura 111.2). Para el péndulo ideal, y para ángulos 
pequeños de alejamiento de la vertical, estas curvas son círcu- 
los concéntricos, que corresponden a la solución más simplifi- 
cada de la ecuación diferencial, llamada en este caso “ecuación 
del oscilador armónico simple”, y que es una ecuación lineal. 

Como puede verse, hay un punto central A, en el cruce de 
los ejes, es decir, para 0 y y de valor cero. Éste representa al 


Velocidad 
Tiempo 


Vía 
S 


FIGURA 11.2. 


Posición O 
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péndulo cuando tiene velocidad cero y descansa sobre la verti- 
cal, o sea que está en reposo. 

Si ahora apartamos mucho de la vertical al péndulo, la re- 
lación entre la fuerza que lo mueve y el ángulo se describe me- 
diante una ecuación diferencial mucho más complicada, y que 
es no lineal. Su solución se grafica en la figura MI.3. La familia 
de posibles trayectorias se asemeja ahora a un ojo. El punto 
central A sigue representando al péndulo inmóvil (velocidad y 
ángulo cero). Las elipses concéntricas B corresponden a ciclos 
del péndulo cada vez más alejados de la vertical, como un co- 
lumpio al que se le da un impulso cada vez mayor, hasta que 
forma 90 grados con la vertical, y luego comienza a subir más 
allá del pivote. 


Velocidad 


Posición 


Figura M1.3. 


¿Qué pasa cuando recibe un envión tal que sobrepasa los 
180 grados? 

Como sabemos, no oscilará, sino que girará en círculo en 
uno u otro sentido, como una hélice. 

Esto se representa en el espacio de las fases con las familias 
de curvas indicadas con Cy D. Como puede verse, las C tienen 
velocidad positiva, es decir, representan giros completos en 
una dirección, por ejemplo en el mismo sentido que las agujas 
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del reloj, y las D representan el movimiento en el sentido con- 
trario, lo que se indica con una velocidad de signo negativo. 

Este método de representar el comportamiento del sistema 
dinámico en el espacio de las fases nos permite apreciar las ca- 
racterísticas esenciales de un vistazo: si el gráfico es un punto, 
el sistema está quieto, si hay círculos concéntricos, está oscilan- 
do con poca amplitud; si hay elipses, las oscilaciones son más 
amplias; 'si está en las zonas Co D, está rotando en lugar de 
oscilar. 

Las dos curvas S que forman el límite entre las regiones B de 
oscilación y C, D de rotación se llaman “separatrices” y corres- 
ponden a las posiciones del péndulo en que está suspendido 
exactamente sobre la vertical en su posición más alta (figu- 
ra 111.3), o sea que tiene la máxima inestabilidad. Si está en esa 
posición caerá rotando ya sea en un sentido o el opuesto, y ad- 
quirirá una velocidad cada vez mayor, hasta pasar por el punto 
más bajo de la vertical y ascender luego con velocidad decre- 
ciente hasta volver al punto más alto sobre la vertical. Sobre la 
curva separatriz, la energía total de que dispone el péndulo es 
exactamente igual a la que necesita para salir de ese punto 
más alto sobre la vertical y volver a alcanzarlo. 

¿Qué pasa con un sistema hamiltoniano más complicado, 
formado por dos péndulos ideales acoplados (A y B), es decir 
que se influencian mutuamente en su movimiento? 

Cada uno de los péndulos tiene su propio periodo de oscila- 
ción, pero éste es ahora afectado por el del otro. Si ignoramos 
al péndulo A, entonces el movimiento del péndulo B trazará 
una curva cerrada en el espacio de las fases bidimensional. Si 
ignoramos al péndulo B, el movimiento de A trazará otra curva 
cerrada en un espacio bidimensional diferente. Pero si los dos 
péndulos interactúan, ya no son independientes y para re- 
presentar las trayectorias los dos espacios planos deben combi- 
narse en uno cuya dimensión debe crecer por lo tanto de dos 
a tres. 


Supongamos que el periodo de A es, por ejemplo, nueve ve- 


ces mayor que el de B. Si A fuese independiente de B, se le po- 
dría representar en la figura 1114 como una curva cerrada 
sobre un plano horizontal, y a Bsobre un plano perpendicular 
al horizontal. Pero cuando sus movimientos se acoplan, esto se 
representa en el espacio de las fases combinando ambas curvas 
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cerradas: a medida que A FIGURA JILA. 
se va desplegando horizon- 
talmente, la curva B la va 
desviando del plano hori- 
zontal, en un movimiento 
comparable al de arrollar 
una soga alrededor de un 
neumático. El resultado de 
un círculo enrollándose al- 
rededor de otro es la crea- 
ción de una figura en forma 
de anillo solenoide, sobre la 
superficie de lo que los ma- 
temáticos llaman toro (fi- 
gura 11.4). Aquí podemos 
considerar que el ciclo A es 
el eje del toro, y el Bes, por 
lo tanto, un ciclo perpen- 
dicular a dicho eje. 

Ahora veamos el toro de tres dimensiones con más detalle. 

Si los periodos o frecuencias de dos péndulos acoplados 
están en una relación simple, por ejemplo si uno tiene un pe- 
riodo nueve veces mayor que el otro, la relación es 9/1, y la 
línea que da vueltas sobre la superficie del toro va trazando un 
solenoide que pasa siempre sobre los mismos puntos del toro 
no importa cuántas revoluciones realicen estos péndulos com- 
binados, lo que demuestra que el sistema dinámico es exacta- 
mente periódico. Así, en la figura 111.4, el punto inicial para 
T=0 y el final al cumplirse un periodo, T'= 1, coinciden. 

Pero ¿qué ocurre cuando los periodos de los osciladores 
acoplados son inconmensurables, o sea que el cociente de esos 
periodos es un número irracional? Un número irracional no 
puede ser escrito como una razón y su expresión decimal con- 
tiene un número infinito de términos sin una pauta repetitiva. 

Si el sistema acoplado tiene una relación irracional entre 
periodos, la curva en el espacio de las fases se irá enrollando 
alrededor del toro pasando en cada vuelta por puntos dife- 
rentes, como se ve en la figura IL5, donde el punto inicial 
T=0 y el que corresponde a un periodo, T'= 1, no coinciden, 
de modo que la curva con el transcurso del tiempo pasará por 


V; 


Periodo A= 1 por segundo 
Periodo B = 9 por segundo 


Sistema periódico 


50 


todos los puntos de su su- FIGURA UL5. 
perficie hasta cubrir total- 
mente el toro, y nunca se 
repetirá a sí mismo. Un sis- 
tema con estas caracterís- 
ticas se llama casi perió- 
dico. 

Los matemáticos son Ca- 
pacés de trabajar con toros 
de cualquier número de di- 
mensiones. Lo que equiva- 
le a que es perfectamente 
posible juntar más de dos 
osciladores y representar su 
movimiento combinado en 
la superficie de un toro mul- 
tidimensional. 


Periodo A/B= número irracional 


Sistema casi periódico 


SECCIONES DE POINCARÉ Y ESPACIOS MULTIDIMENSIONALES 


Si bien no hay ninguna restricción en principio al número de 
dimensiones del espacio de las fases, evidentemente son mu- 
cho más fáciles de visualizar las formas con solamente tres o 
dos dimensiones (volúmenes o superficies). 

Hay un método, concebido por H. Poincaré, para visualizar 
las propiedades esenciales de trayectorias complicadas de espa- 
cios de tres o más dimensiones, que consiste en rebajar en uno 
ese número de dimensiones. 

Consideremos el ejemplo de las trayectorias que rodean el 
toro de un sistema biperiódico y que forman la superficie de 
un solenoide tridimensional. 

Cortemos ahora el toro con un plano transversal (figura IM1.6) 
y marquemos sobre él los puntos donde la trayectoria intersec- 
ta el plano. Como la trayectoria es periódica, se marcará un 
punto por vuelta, de modo que después de un tiempo suficien- 
temente largo se tendrá un verdadero mapa sobre el plano, lo 
que resulta muy ventajoso en la práctica, por la simplifica- 
ción que significa el rebajar en uno el número de dimensio- 
nes, y además porque se ha pasado de la descripción continua 
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en el tiempo de la trayecto- FIGURA IL. 
ria en el espacio de las fa- 
ses, a tomar solamente los 
datos de movimiento cada 
vez que la trayectoria atra- 
viesa esa sección, y eso im- 
plica que es necesario ma- 
nejar un número de datos 
mucho menor. 

¿Qué nos puede decir la 
sección de Poincaré para 
este sistema biperiódico? 

Si la relación de frecuen- 
cias es un número racional, 
la curva está fija sobre el 
toro, cada vuelta se super- 
pone a la anterior, de mo- 
do que en la sección apare- 
cen solamente los correspondientes puntos aislados cuyo 
número y posición dependen de la relación de periodos. 

Si, en cambio, el sistema es casi periódico, con una relación 
de frecuencias que es un número irracional, la curva pasará 
por un punto distinto en cada vuelta, cubriendo toda la super- 
ficie del toro con el tiempo, de modo que en la sección de 
Poincaré se tiene una curva cerrada. 

El análisis de sistemas de osciladores más complicados re- 
quiere la introducción de un espacio de las fases de más dimen- 
siones, ya que como hemos visto, el número de dimensiones 
depende de la cantidad de variables independientes del sis- 
tema, que pueden ser la velocidad o el impulso, la posición, o 
algunas otras características dinámicas que definan el compor- 
tamiento de cada uno de los componentes. 

Para poder encarar estos conceptos de espacio multidimen- 
sional, cruciales para el estudio de los sistemas complejos, se 
requiere hacer una generalización de la geometría de coorde- 
nadas. ) 

Decimos que vivimos en un espacio tridimensional, ya que 
nuestro movimiento en el espacio tiene tres grados de liber- 
tad: podemos hacer tres tipos de movimiento que tengan di- 
recciones perpendiculares entre sí (izquierda-derecha, delan- 
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te-detrás, arriba-abajo), y cualquier punto del espacio puede ser 
alcanzado combinando esos tres tipos posibles de movimiento, 
de modo que su posición puede ser indicada mediante tres 
números que llamaremos coordenadas x, y, z del punto, y que 
dan las distancias del punto a una referencia en esas tres direc- 
ciones perpendiculares. 

También podemos referirnos a un espacio abstracto de cua- 
tro dimensiones, con cuatro coordenadas w, x, y, z, a uno de 
cinco dimensiones, con coordenadas vu, 14, Xx, y, 2, y así sucesiva- 
mente, siempre y cuando tengamos en cuenta que ya no nos 
estamos refiriendo al espacio físico tridimensional en el que 
vivimos y nos movemos, sino a un espacio matemático. Este 
resulta de enorme utilidad cuando se trata de comprender las 
matemáticas de muchas variables. 

En todo problema, sea de física, biología, economía, cual- 
quier magnitud significativa puede ser considerada, y visualiza- 
da, como una dimensión del problema. Un economista puede 


trabajar con un “espacio” multidimensional con variables ín- . 


dice de costo de vida, costo de la vivienda, valor del dólar, pre- 
cio del petróleo, trimestres de la última década, etcétera. 

Un físico puede estudiar un sistema dinámico formado por 
tres cuerpos, por ejemplo, donde a cada uno le corresponden 
las tres coordenadas de posición x, y, z más las tres coordena- 
das de velocidad o de impulso, o sea un total de 18 coorde- 
nadas, de modo que está tratando con un espacio de 18 di- 
mensiones. 

Nada nos impide continuar aumentando el número de com- 
ponentes del sistema, de modo que podemos decir que, en ge- 
neral, un sistema dinámico con » variables independientes 
—n grados de libertad— se puede representar en un espacio 
de n dimensiones. 

En general, esto es válido para los sistemas ordenados y esta- 
bles que, aunque estén formados por un gran número de com- 
ponentes, y que por lo tanto se los debería representar en un 
espacio de las fases de gran número de dimensiones, en la 
práctica se mueven en un muy pequeño subespacio de este vas- 
to espacio multidimensional que representa un estado físico 
en todos sus menores detalles, y cuyo volumen llamaremos H. 
Así un cuerpo sólido, por ejemplo, una roca, si bien está for- 
mado por una enorme cantidad de moléculas, como las mis- 
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mas están rígidamente li- 
gadas entre sí se mueven 
todas al unísono, y basta 
con representar el movi- 
miento de un único pun- 
to: el centro de masa de 
la roca. 

La evolución del siste- 
ma con el tiempo se repre- 
senta con la trayectoria 
de H en el espacio de las 
fases (figura 1.7), don- 
de se marcaron los pun- 
tos que corresponden a 
los estados para los tiem- 

FIGURA HI. 7. pos. T=L Tel Pad, 

etcétera. 

E Un sistema dinámico de comportamiento ordenado puede, 
40 si se alteran sus características, pasar a tener un comporta- 
1] miento caótico, y también puede ocurrir el proceso contrario. 
¿Es posible representar esos cambios en el espacio de las fa- 
lh ses? La respuesta es afirmativa: como veremos más adelante, el 
' estudio de la transición de un sistema dinámico del orden al 
| 
1 


! caos es, en cierto sentido, el análisis de cómo un movimiento 
ps que puede ser muy simple, limitado y repetitivo se quiebra en 
cierto punto crítico, desarrollando un nuevo comportamiento 
que corresponde a un desplazamiento de la trayectoria del sis- 
“tema a zonas mucho más vastas del espacio de las fases. 

Ahora podremos visualizar la diferencia entre procesos pre- 
decibles o procesos caóticos. Representemos (figura 111.8) la 
evolución de un sistema dinámico de comportamiento pre- 
decible como una trayectoria indicada por 1 en el espacio de 
las fases. Si variamos ligeramente las condiciones iniciales, te- 
nemos la trayectoria 2, que se mantiene cercana a la 1, indi- 
cando que el comportamiento es prácticamente el mismo (fi- 
gura 111.84). En cambio, en un proceso caótico, las trayectorias 
1 y 2 que inicialmente estaban próximas se alejan cada vez más 
con el tiempo, señalando la divergencia creciente de los com- 
portamientos (figura 111.80). 
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FIGURA IL8. 


CONSERVACIÓN DE LOS VOLÚMENES 
EN EL ESPACIO DE LAS FASES 


Como hemos visto, la energía total de un sistema hamiltoniano 
es invariable. ¿Cómo se expresa esto en el espacio de las fases? 

Consideremos la representación en el espacio de las fases de 
la evolución del sistema dinámico. Vimos que hay una región 
de volumen HA para representar ese sistema, y que se va despla- 
zando en el espacio de las fases, trazando así una trayectoria 
que expresa su evolución en el tiempo. 

En el siglo pasado, el matemático J. Liouville demostró que 
para todo sistema hamiltoniano el volumen de esta región A se 
mantiene constante en el tiempo, tal como si se tratase del movi- 
miento de una gota de un líquido que no es posible comprimir. 

Esta propiedad de los sistemas hamiltonianos es mucho más 
restrictiva que la simple conservación de la energía o que la 
reversibilidad en el tiempo de las ecuaciones de movimiento. 

La conservación del volumen de la región en el espacio de 
las fases puede estar asegurada de dos maneras diferentes: 1) la 
región considerada se va desplazando alrededor de la trayecto- 
ria, rotando y deformándose en forma periódica, y entonces 
las trayectorias vecinas se “enrolian” sin apartarse nunca dema- 
siado entre sí; 2) el volumen A se estira con el tiempo en una 
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dirección, contrayéndose en la dirección perpendicular. Mien- 
tras que en el primer caso dos trayectorias inicialmente próxi- 
mas se mantienen vecinas, en el segundo tienden a alejarse. 

Desde el punto de vista dinámico, la diferencia es conside- 
rable. En efecto, las trayectorias son estables en el primer caso, 
inestables en el segundo, pues aquí un débil apartamiento int- 
cial se puede amplificar en el transcurso del tiempo. 


REVERSIBILIDAD EN EL TIEMPO 


Otra propiedad importante de las ecuaciones que describen el 
comportamiento de un sistema hamiltoniano es que el cambio 
de signo del tiempo, es decir, el reemplazo de +1 por —t no tie- 
ne incidencia alguna; las ecuaciones son idénticas. O sea que si 
se filma el movimiento de un péndulo ideal no se puede dis- 
cernir en qué sentido se pasa la película. 

Se dice, en general, que los sistemas conservativos tienen 
una mecánica reversible. 

Como veremos más adelante, la irreversibilidad en el tiempo 
es característica de los sistemas en los que la energía, en lugar 
de conservarse, se disipa. 


DONDE APARECE EL CAOS EN LOS PÉNDULOS 


El estudio del caos determinista comenzó en los años sesenta 
con los trabajos pioneros de E. Lorenz, D. Ruelle y F. Takens, 
quienes utilizando las computadoras, que en ese entonces co- 
menzaban a mostrar su fantástica potencialidad, demostraron 
que aun sistemas dinámicos simples, formados por unos po- 
cos osciladores, podían comportarse de manera imprevisible, 
caótica. 

Es fácil y divertido comprobar esto con un sistema de dos 
péndulos acoplados—sencillo de construir y hasta adquirible 
como juguete. En la figura ULO9, el sistema consta de un pén- 
dulo liviano formado por dos pequeñas esferas unidas por un eje 
alrededor del cual pueden rotar. Este eje cuelga sobre un pén- 
dulo más pesado, que es el que impone la oscilación básica. 

Ambos péndulos tienen adosados imanes permanentes, de 
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FIGURA IM.9. 


modo que sus movimientos se acoplan. En la base del sistema 
hay un pequeño electroimán alimentado por un circuito eléc- 
trico oscilador, que mantiene las oscilaciones del péndulo prin- 
cipal para que éste no se amortigúe (péndulo entretenido). 

Una vez que se da el empujón inicial, el péndulo más pesa- 
do oscila con la regularidad de un reloj, mientras que cada vez 
que una de las esferas del otro gira cerca del péndulo grande 
recibe un envión por la atracción entre los respectivos imanes. 
Pronto se asistirá a un espectáculo sorprendente: el péndulo 
liviano hace una extraña danza errática, oscilando por momen- 
tos en forma rítmica, para saltar imprevisiblemente a movi- 
mientos que son caóticos. 

¿Cómo se representará este comportamiento en el espacio 
de las fases? 
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SISTEMAS DISIPATIVOS Y ATRACTORES 


Tratemos el caso de los sistemas dinámicos disipativos, que 
abundan en este mundo, y cuya energía va disminuyendo por 
fricciones y otros efectos. 

Sabemos que si le damos impulso a un péndulo real, oscilará 
o rotará pero pronto, a diferencia del péndulo ideal que 
hemos estudiado antes, su movimiento se irá amortiguando 
cada vez más, hasta que finalmente quedará quieto, si es que 
no recibe nuevamente energía. 

Este comportamiento queda representado como una trayec- 
toria en espiral en el espacio de las fases, que desemboca en el 
punto A para la posición final de reposo (figura M1.10). . 


Impulso 


Posición 


FIGURA 111.10. 


No importa cuánto impulso inicial recibe el péndulo, en to- 
dos los casos la pérdida de energía finalmente lo inmoviliza, es 
decir, la trayectoria desemboca inevitablemente en el punto de 
reposo A, como si éste atrajera las curvas en el espacio de las fa- 
ses. De allí el nombre de “atractor” que recibe el punto fijo A. 

Por supuesto que un reloj a péndulo sería de muy poca utili- 
dad si esta trayectoria desde el primer impulso al reposo se 
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hiciera en unos pocos minutos, y por ello se han inventado 
diversos mecanismos que reponen la energía que va perdien- 
do al funcionar: pesas que caen resortes elásticos, electro- 
imanes que cambian cíclicamente su polaridad. El resultado es 
que el péndulo se mueve a un ritmo regular a pesar de la fric- 
ción y de la resistencia del aire. 

Por lo tanto, en un péndulo que recibe energía, las curvas 
corresponden en el espacio de las fases a la figura 111.11. 


FIGURA 111.11. 


De hecho, si al péndulo se le da un empuje adicional, o se le 
frena momentáneamente, volverá eventualmente a su ritmo 
original, que corresponde a la trayectoria C. Esta curva consti- 
tuye evidentemente un nuevo tipo de atractor, ya que en lugar 
de ser atraído el sistema a un punto fijo, es llevado a una 
trayectoria que forma una curva cerrada. 

Obsérvese que hay una diferencia importante con el péndu- 
lo ideal que se mueve sin fricciones ni pérdida de energía: en 
éste, la más pequeña perturbación causada al agregarle un 
impulso mayor o al frenarlo, hace que cambie la órbita del 
péndulo para contraerse o expandirse un poco, o sea salta de 
una curva a otra concéntrica más grande o más pequeña. En 
contraste, la trayectoria de un péndulo ayudado mecánica- 
mente tiene estabilidad, resiste pequeñas perturbaciones (figu- 
ra 111.11), y así cuando se le da un impulso mayor lo va disipan- 
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do de modo gradual, y si se lo frena va recibiendo energía de 
su fuente, por lo que en ambos -asos regresa finalmente a esa 
única curva cerrada, que por lo tanto es también un atractor, 
al que llamaremos “atractor de ciclo límite”. 

Hay dos clases básicas de ciclos límites: el que acabamos de 
presentar para el péndulo, es estable; los puntos de trayecto- 
rias próximas se mueven hacia él. Hay también sistemas di- 
námicos con ciclos límite inestables, donde los puntos de tra- 
yectorias cercanas se alejan de esa curva, que actúa como un 
repulsar. 

La importancia de lo que hemos analizado para el péndulo 
es que éste permite comprender las características esenciales 
del comportamiento de los sistemas que actúan cíclicamente y 
que son tan frecuentes en la naturaleza: un oscilador eléctrico, 
las mareas, las vibraciones del aire en un tubo de órgano, los 
impulsos eléctricos que hacen latir el corazón, la cantidad de 
individuos en una población de animales... 

Recordemos los dos tipos de atractores descritos hasta aho- 
ra: 1) el punto atractor, que corresponde a un estado estaciona- 
rio del sistema, nada ocurre al transcurrir el tiempo; 2) el 
atractor de ciclo límite, que indica un comportamiento periódi- 
co, lo que implica, además, que, si bien el sistema es disipativo 
y, por lo tanto, va perdiendo su energía, ésta se va reponiendo 
por la entrega de energía de alguna fuente exterior. 


CONTRACCIÓN DE LOS VOLÚMENES 
EN EL ESPACIO DE LAS FASES 


¿Cómo queda representado en el espacio de las fases el hecho 
que estamos tratando aquí con sistemas disipativos? Como vi- 
mos, el atractor de ciclo límite nos indica que se está ante un 
sistema disipativo estable, que repone la energía que va per- 
diendo. Si se tiene, por ejemplo, un péndulo entretenido y se 
le perturba frenándolo temporariamente, al disminuir el án- 
gulo y la velocidad se obtiene la trayectoria B (figura 111.12), 
pero el sistema recibe energía de su fuente para compensar 
esta disminución y la trayectoria termina por confundirse con 
el ciclo límite A. 

Lo mismo ocurre si se da al péndulo un impulso adicional 
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que lo saque de su órbita; con el tiempo, la disipación de esta 
energía en exceso hará converger su trayectoria sobre el ciclo 
límite A. 

Por lo tanto, el atractor de ciclo límite está dentro de una 
zona C del espacio de las fases tal que toda trayectoria que se 
inicie en un punto cualquiera de dicha región termina por ser 
atraída inexorablemente por el atractor. Esta zona se denomi- 
na “cuenca de atracción” (figura 111.12). 


atracción 


FIGURA 107.12. 


Supongamos que en esta cuenca C haya una región R que 
representa un conjunto de valores iniciales de posiciones e im- 
pulsos para el sistema. A medida que las trayectorias son atraí- 
das por el ciclo límite, se van acercando entre sí y acaban por 
converger sobre el mismo, quedando finalmente una única 
trayectoria: hay una contracción de la región R, la que va dis- 
minuyendo hasta desaparecer en la curva del atractor. 

Como una curva es una línea de una sola dimensión, no 
puede pasar por todos los puntos que conforman un volumen; 
siempre quedan infinitos puntos no cubiertos por la curva y es 
evidente, entonces, que en un espacio de las fases de tres di- 
mensiones un atractor debe tener una dimensión inferior a 3, 
y esto se puede generalizar: 
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La dimensión d de un atractor en un espacio de n dimensiones es 
menor que n: 
d<n 


Ésta es la característica principal que distingue a los sistemas 
disipativos de los conservativos. En un sistema disipativo todas 
las condiciones iniciales convergen hacia regiones del espacio 
de las fases que tienen una dimensión menor que la del espa- 
cio original. 


SISTEMAS DISIPATIVOS FORMADOS 
POR DOS PÉNDULOS ACOPLADOS 


Pasemos ahora a estudiar el caso de un sistema disipativo for- 
mado también por dos péndulos acoplados. Aquí, de modo 
similar a lo visto en el estudio de sistemas conservativos forma- 
dos por la combinación de dos péndulos ideales, se requerirá 
un espacio de las fases de tres dimensiones, en el que las tra- 
yectorias se desarrollan sobre la superficie de un toro, de modo 
que tenemos para el caso de un sistema disipativo un atractor 
más complicado que el de punto fijo o el de ciclo límite: se tra- 
ta de una curva en forma de solenoide que pasa por todos los 
puntos de la superficie del toro, representando así el compor- 
tamiento combinado de dos péndulos acoplados de periodos 
que son inconmensurables, o sea cuyo cociente es un número 
irracional. Por lo tanto, lo llamaremos “atractor casi perló- 
dico”. 

El comportamiento de tal sistema es predecible, es decir, 
que conociendo las velocidades y posiciones de sus compo- 
nentes en un momento dado, se lo puede determinar para 
cualquier otro instante. Y aun si tal conocimiento tiene un 
cierto margen de error, o incertidumbre, éste se mantiene en 
el mismo orden de magnitud para toda determinación futura. 
Tal propiedad se traduce en el espacio de las fases por el 
hecho de que dos trayectorias vecinas siguen siendo vecinas a 
medida que pasa el tiempo. 
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Los ATRACTORES EXTRAÑOS 
Y LAS DIMENSIONES FRACTALES 


Volvamos ahora a un sistema que presenta un comportamien- 
to caótico, como el de los dos péndulos acoplados por un 
imán que habíamos descrito. 

Su representación en el espacio de las fases también reque- 
rirá que éste tenga tres dimensiones, en el que la trayectoria 
será una curva continua, similar al caso del atractor casi perió- 
dico. Pero lo que lo distingue de éste es que si uno examina 
dos trayectorias vecinas en el espacio de las fases, se ve que di- 
vergen rápidamente, alejándose cada vez más. Se trata aquí de 


un “atractor extraño”, cuya característica esencial es la amplifi- 


cación de los apartamientos, por mínimos que éstos sean, en- 
tre trayectorias en el espacio de las fases. Esta característica se 
denomina sensibilidad a las condiciones iniciales. He aquí la clave 
para comprender por qué el determinismo que rige a un sis- 
tema dinámico no implica necesariamente la predictibilidad. 

Si aparece esta sensibilidad, el sistema es impredecible des- 
pués de un cierto tiempo, no importa cuáles sean las otras ca- 
racterísticas del espacio de las fases, ni tampoco su número de 
dimensiones. Sólo podríamos predecir su evolución con exac- 
titud si conociésemos con precisión infinita absolutamente 
todos los factores que actúan sobre el sistema, y sabemos que 
esto es imposible, ya que necesariamente conocemos las condi- 
ciones iniciales sólo en forma aproximada. 

En cambio, las situaciones correspondientes a los otros atrac- 
tores pertenecen a estados estables y predecibles, porque no 
tienen dicha sensibilidad, sino todo lo contrario: hay una in- 
sensibilidad a las condiciones iniciales. 

Esto se evidencia porque, de acuerdo con la figura anterior, 
hay en el espacio de las fases una cuenca de atracción dentro 
de la cual una trayectoria que parta de cualquier condición ini- 
cial termina inexorablemente por fundirse en el atractor; no 
importa cómo se inicie el proceso, el sistema acaba después de 
un tiempo por seguir el comportamiento estable y previsible 
descrito por dicho atractor. 

La sorpresa que se experimenta ante la aparición del caos 
en un sistema tan simple como el formado por dos osciladores 
acoplados es natural, pues a nadie le resultaba novedoso o ex- 
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traño que sistemas tan complicados como la atmósfera terres- 
tre o un torrente turbulento de agua tengan un comporta- 
miento difícil de pronosticar. Sin embargo, el estudio de estos 
sistemas simples ha arrojado luz sobre el de los formados por 
un gran número de componentes. 

El atractor que caracteriza el comportamiento caótico se 
transforma de manera espectacular y aparentemente contraria 
a la intuición, ya que debe reflejar en su geometría esa sl- 
tuación. 

Su estructura debe exhibir dos tendencias opuestas: hacien- 
do honor a su nombre, las trayectorias vecinas deben conver- 
ger hacia el atractor y contrariamente, para reflejar el estado 
de sensibilidad a las condiciones iniciales, las trayectorias de- 
ben diverger distanciándose cada vez más. La velocidad con que 
divergen estas trayectorias sucesivas se mide con un coeficiente 
llamado “exponente de Lyapunov”. 

Además, la otra condición esencial es que, como hemos vis- 
to antes, las curvas formadas por las trayectorias no se pueden 
cruzar tocándose en un punto (condición de determinismo). 
Resulta entonces un atractor muy difícil de visualizar, cuya for- 
ma no puede estar sobre una superficie, no importando que 
sea esta plana o curva como el solenoide que rodea el toro en 
el caso simple que examinamos. 

La razón de esta imposibilidad es que la divergencia entre 
trayectorias adyacentes a la que nos obliga el estado caótico no ' 
se puede lograr sobre la superficie del toro, porque tras una 
cierta cantidad de vueltas éstas, al no ser paralelas, deberían 
cruzarse para poder alojar todas las trayectorias. No alcanza 
entonces el espacio que hay en la superficie y existe una única 
solución alternativa, que es la que aparece en el atractor: las 
trayectorias se despegan de la superficie saltando hacia afuera 
para ocupar la región que la rodea. Por lo tanto las dimen- 
siones de este atractor deben ser mayores que las dos que co- 
rresponden a una superficie, o sea: 


2<d 
Pero, y aquí surge la dificultad, estamos tratando un sistema 


simple, en un espacio de tres dimensiones, y hemos visto que, 
como las condiciones del sistema dinámico (velocidades, posi- 
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ciones, etcétera) están necesariamente limitadas a cierto rango 
de valores, el atractor debe ocupar sólo una zona restringida 
del espacio de las fases, en lugar de llenarlo por completo, y 
entonces su dimensión d debe ser menor que 3: 


2<d<3 


¡O sea que debemos imaginar una figura geométrica de di- 
mensión mayor que 2 y menor que 3, es decir que esté en una 
situación intermedia entre una superficie y un volumen! Ob- 
viamente, una situación así no es normal en la geometría de 
Euclides, y éste puede haber sido uno de los motivos de que se 
haya descubierto hace tan poco la posibilidad del caos en sis- 
temas simples. 

Pero recientemente se han desarrollado matemáticas en las 
que sí existen formas irregulares o fragmentadas que se pue- 
den caracterizar con “dimensiones” que, a diferencia de las eu- 
clidianas no son números enteros, y que han sido llamadas 
fractales por el matemático Benoit Mandelbrot, que ha sido uno 
de los principales impulsores para el estudio de estas extrañas 
geometrías. 

Este atractor es entonces un fractal, está localizado en una 
región que incluye la superficie del toro pero, como es una cur- 
va, no ocupa todos los puntos del volumen de dicha región, 
quedando siempre infinitos puntos por los que no pasa. 

Otros sistemas dinámicos, como los formados por oscilado- 
res eléctricos, fluidos tur- 
bulentos, reactivos quími- 
cos, exhiben atractores con 
estas características y que 
han sido bautizados “atrac- 
tores extraños” por Ruelle 
en 1971. 

Un atractor extraño y su 
sección de Poincaré pue- 
den tener el aspecto ge- 
neral de la figura 111.13, 
que podríamos comparar 
con el de un ovillo de hilo 
después de que un cacho- FIGURA HIT. 13. 
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rro lo haya usado para jugar durante unas horas: si bien sigue 
siendo un hilo único (por lo tanto, de trayectoria determina- 
da), será imposible al seguir las vueltas prever si un centímetro 
más adelante va a replegarse, ir hacia adentro del ovillo, o ha- 
cia afuera, etcétera. Como tiene sensibilidad a las condiciones 
iniciales, la más ínfima alteración de éstas se representará con 
otro ovillo enmarañado cuyas vueltas no tienen nada que ver 
con el primero, aunque el volumen que ocupa sea práctica- 
mente el mismo. Para obtener la forma de un atractor extraño 
vamos a utilizar el siguiente procedimiento: consideremos el 
flujo en tres dimensiones de las trayectorias en el espacio de 
las fases (figura 111.14). 


Dirección 
de contracción 


Dirección 
de estiramiento 


iniciales Qu 


Condiciones 


FIGURA 111.14. 


Este flujo está sometido a dos influencias opuestas, ya que se 
debe contraer porque hay un atractor, y se debe estirar por la 
sensibilidad a las condiciones iniciales. 

Para analizar mejor el proceso, separemos los dos efectos, 
de manera que la contracción por la condición de ser un atrac- 
tor se ejerza en la dirección vertical, mientras que el de esti- 
ramiento se ejerza en la dirección horizontal. 

Secciones transversales de Poincaré sucesivas muestran, por 
lo tanto, un rectángulo que se va deformando, contrayéndose 
verticalmente y dilatándose horizontalmente. Este proceso de 
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contracción y estiramiento debe continuar a medida que el 
flujo de trayectorias recorre el atractor. Si se mide el estira- 
miento en cada vuelta sucesiva y se compara con el anterior, 
cuando crece exponencialmente con un factor (exponente de 
Lyapunov) positivo, hay realmente estiramiento de la trayec- 
toria. 

Pero hemos visto que el volumen que ocupa el flujo en el es- 
pacio de las fases debe mantenerse constante para sistemas ha- 
miltonianos, porque la energía se conserva, o debe disminuir 
para sistemas disipativos, que pierden energía por fricción. Por 
lo tanto, el área ocupada por las sucesivas deformaciones de 
dicho rectángulo no puede aumentar, sólo puede mantenerse 
prácticamente constante si hay poca disipación de energía o 
disminuir, según el caso. 

Además, aparece aquí una limitación adicional, que es la res- 
ponsable de la aparición de la forma fractal: como las variables 
que describen el sistema (impulsos, posiciones, etcétera) no 
pueden tener cualquier valor imaginable, sino que están nece- 
sariamente restringidas, el flujo debe mantenerse confinado a 
una región del espacio de las fases y, por lo tanto, hay un límite 
para el estiramiento del rectángulo. 

En consecuencia, la única solución que permite la geome- 
tría para cumplir todas estas condiciones asegurando al mismo 
tiempo que el flujo se mantenga en una zona restringida, es 
que el rectángulo se repliegue sobre sí mismo. 

Al mantenerse de manera simultánea las tres operaciones, 
contracción, estiramiento y replegado, el rectángulo se trans- 
forma progresivamente en una herradura que, a su vez, se 
aplanará, estirará, plegará, dando nacimiento a una estructura 
de doble horquilla, y así sucesivamente. El área de la figura 
resultante en cada deformación va disminuyendo, o puede 
mantenerse prácticamente constante según el caso, pero no 
puede crecer. 

En cada estiramiento aumenta la distancia en xentre puntos 
que antes eran contiguos, que se separan entonces exponen- 
cialmente, como corresponde a la condición de sensibilidad a 
las condiciones iniciales medida por el coeficiente de Lyapu- 
nov. El atractor extraño se fabricó así de manera similar a la que 
usa el panadero para la masa del pan, y no es sorprendente 
que su estructura sea por lo tanto parecida a la del hojaldre. 
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FIGURA IH1.15. 


Si se examina la sección de Poincaré del atractor de la figura 
anterior, se ve una estructura de bandas que se repite sin fin 
(figura 111.15). 

Cada banda está compuesta por subbandas y éstas, a su vez, 
por otras de estructura similar; sin importar qué escala utilice- 
mos para examinar la microestructura, su aspecto será seme- 
jante. Esta propiedad se llama de “autosemejanza” y es una ca- 
racterística de las formas fractales, las que como hemos visto, 
tienen dimensiones que no son números enteros. 

En resumen, para que el comportamiento de un sistema di- 
námico pueda transformarse en caótico basta con que tenga 
un mínimo de tres grados de libertad y, por lo tanto, un espa- 
cio de las fases con un mínimo de tres dimensiones. Si en el 
correspondiente atractor encontramos que la sección de Poin- 
caré presenta estructuras que son autosemejantes, sabremos 
que el atractor es extraño, lo que implica sensibilidad a las 
condiciones iniciales y, por lo tanto, imposibilidad de predecir 
el comportamiento del sistema más allá de un cierto tiempo. 

Estos atractores extraños, curvas que no tienen fin y que se 
alojan en el espacio de las fases, son, en general, figuras geo- 
métricas de una rara belleza. No es posible calcularlos con exac- 
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titud con base en ecuaciones matemáticas, ya que ninguna los 
puede describir en forma precisa, y la única vía para construir- 
los y visualizar su aspecto en el espacio de las fases es mediante 
la computadora, tal como fueron descubiertos. 

Es importante señalar que los atractores extraños son figuras 
que ocupan sólo una zona en el espacio de las fases, y esto per- 
mite diagnosticar de un vistazo si un sistema que aparece 
como caótico tiene regularidades subyacentes o si su compor- 
tamiento es puramente aleatorio. 

Si se puede construir un atractor con base en los datos, esto 
nos indica que hay algún mecanismo no lineal operando sobre 
el sistema y que, por lo tanto, se trata del caos determinista. Si, 
en cambio, el comportamiento representado da una disper- 
sión de puntos por todo el espacio de las fases, esto indica que 
se trata de fenómenos que ocurren al azar, lo que técnicamen- 
te se llama ruido aleatorio. 


LAS DIMENSIONES FRACTALES 


¿Qué es un fractal exactamente, y cómo se hace uno? 


Un fractal es una forma geométrica que consiste en un motivo que: 


se repite a sí mismo en cualquier escala a la que se le observe. 


Esta forma puede ser muy irregular, o muy interrumpida o frac- 
cionada, y de aquí el origen de su nombre, que B. Mandelbrot 
derivó del latín “fractus” (interrumpido o irregular). Su carac- 
terística básica es el concepto de la autosemejanza, y por ello la 
geometría de los fractales es una herramienta imprescindible 
para el estudio de todos aquellos fenómenos que exhiban una 
misma estructura no importa cuál sea el aumento con que se 
les examine. 

Esta propiedad aparece con sorprendente frecuencia en la 
naturaleza: pensemos en un torrente turbulento con remoli- 
nos, que si se examinan en detalle contienen remolinos más 
pequeños, o en muchas plantas, helechos, árboles, que se ra- 
mifican por procesos sucesivos hasta llegar a las hojas, o en el 
sistema de circulación sanguínea con sus ramificaciones cada 
vez más pequeñas hasta llegar a los capilares. 
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Los fractales son, al mismo tiempo, muy complejos y particu- 
larmente simples. Son complejos en virtud de su detalle infini- 
to y sus propiedades matemáticas únicas (no hay dos fractales 
iguales); sin embargo, son simples porque pueden ser genera- 
dos por la aplicación sucesiva de una simple iteración, y la in- 
troducción de elementos aleatorios. ; 

Los matemáticos han concebido gran variedad de fractales, 
pero el más antiguo es probablemente , el del matemático G. 
Cantor, en 1883. Este quería sorprender a sus colegas con dos 
características aparentemente contradictorias para un conjun- 
to de números comprendidos entre el 0 y el 1: a) que el con- 
junto tuviera medida cero, o sea que si se le representase con 
puntos sobre una recta, en cualquier porción de la misma 
los puntos que no pertenecen al conjunto fuesen muchísimos 
más que los que son parte de él, y, simultáneamente, 5) que la 
cantidad de miembros del conjunto fuera tan innumerable 
como el conjunto de todos los números reales también com- 
prendidos entre 0 y 1. 

Muchos matemáticos, incluyendo al mismo Cantor al princi- 
pio, no creían que semejante monstruo pudiese existir, pero 
finalmente él lo descubrió. Su construcción parece asombrosa- 
mente sencilla. Comencemos con un segmento de recta y mar- 
quemos los puntos extremos con 0 y 1. Borremos el tercio in- 
terno, pero manteniendo sus puntos extremos 1/3 y 2/3. 
Quedan dos segmentos con un total de cuatro puntos extremos, 
y a cada uno de éstos le borramos nuevamente el tercio inte- 
rior, con lo que quedan cuatro segmentos con dos puntos 
extremos cada uno, o sea 
puntos que correspon- y 4 
den a los 8 números: 0, 


1/9, 2/9, 3/9, 6/9, 7/9 4 1/3 2/3 1 

8/9 y 1. Si continuamos == A 

con el mismo procedi- o 08 A A 10 

miento ad infinitum, lle- “—= mu a 

garemos eventualmente a 

que cada segmento esté mm ma as 

formado por un único 

punto. Ba 28 BY E BB 8 se 29 
En la figura 1.16 se 

muestran las primeras cin- FiGURA 111.16. 
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co Operaciones de borrado, pero no hay forma de dibujar el 
resultado final. El conjunto de Cantor tiene agujeros en cual- 
quier escala en que se le examine, y está formado únicamente 
por infinitos puntos aislados, ninguno de los cuales tiene por 
vecino a Otro punto del conjunto. Si se quiere medir la longi- 


tud que queda, se comprueba que es cero por estar hecha de. 


una suma de agujeros, ya que los segmentos que se fueron qui- 


tando suman un total de longitud 1 después de las infinitas 


operaciones de borrado. 

- Los matemáticos dicen entonces que el conjunto de Cantor 
tiene medida o dimensión cero. Pero actualmente hay un nue- 
vo concepto matemático que tiene más significado para medir 
las dimensiones de los fractales, y de acuerdo con éste se de- 
muestra que la dimensión del conjunto de Cantor es la llama- 
da “dimensión fractal D” 


_log2 
log 3 


D = 0.6309 


que resulta ser no un número entero sino una fracción, y que 
indica una forma geométrica intermedia entre puntos (dimen- 
sión euclideana 0), y una curva (dimensión euclidiana 1). 

Hay que aclarar que no todos los fractales tienen una D que 
es una fracción; éste puede ser entero, pero su forma es siem- 
pre interrumpida o irregular. 

Tal como ocurre con las dimensiones del espacio de las fa- 
ses, se debe tener presente que una dimensión fractal no tiene 
el mismo significado que el de las dimensiones de nuestro es- 
pacio euclidiano, sino que es la expresión numérica que nos 
permite medir el grado de irregularidad o de fragmentación 
de una de estas figuras. Así, una dimensión D entre 1 y 2 signi- 
fica que se trata de ciertas curvas planas muy irregulares que 
casi llenan un plano, y las superficies que parecen hojaldres, 
con numerosos pliegues que llenan parte de un volumen tie- 
nen un D entre 2 y 3. Hay fractales con D= 1 y con D= 2, que 
no se parecen en nada a una línea o a un plano, pero siempre 
son irregulares o interrumpidos. 

Volvamos ahora a examinar el atractor extraño correspon- 
diente a la transformación de la herradura. 
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FIGURA 11.17. 


Esta transformación fue desarrollada por el matemático S. 
Smale, para aplicarla al estudio de osciladores electrónicos 
acoplados que se utilizaban en los equipos de radar. Al hacer 
una sección s-s* (véase la figura 11.17) se observa el conjunto 
de Cantor. 

Las propiedades esenciales de este atractor aparecen en un 
gran número de secciones de Poincaré para sistemas caóticos. 
El conjunto de Cantor es un ejemplo de fractal interrumpido. 

Un fractal, ejemplo clásico de una curva irregular, infinita- 
mente fracturada, es la curva de Koch, propuesta por el mate- 
mático sueco H. von Koch en 1904. Para construirla, se parte 
de un segmento de recta (figura 111.18) de longitud 1 y en su 
tercio medio se construye un triángulo equilátero. La longitud 
de esta línea es ahora 4/3. Si se repite la operación se obtiene 
la figura con una longitud (4/3)?, o 16/9, y reiterando infini- 
tas veces se llega a una forma fractal de longitud infinita y cu- 
yos extremos están, sin embargo, separados por la misma dis- 
tancia que el segmento generador inicial de longitud 1. Su 
dimensión fractal es D= 1.26... 

Otra variante es el copo de nieve de Koch que se construye 
con el mismo método a partir del triángulo equilátero (figu- 
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[e] 
mb 


Longitud 1 


1/3 1/3 


1/3 1/3 


_ 


0 Longitud 4/3 


o 
— 


Longitud (4/3)? 


o 
sb 


Longitud (4/3)" 


FIGURA 111.18. 


ra 111.19), y en el que, si cada lado del triángulo mide 1/3, 
después de n iteraciones se tiene un perímetro de longitud 
total (1/3)”, que se hace infinito cuando n llega a infinito, de 
modo que si uno quisiera recorrer idealmente con un lápiz 
todas las vueltas de la curva, no llegaría jamás al final, a pesar 
de que encierra una figura hexagonal de área perfectamente 
limitada. 

Estas extrañas figuras parecen productos de la imaginación 
de algunos matemáticos desconectados de la vida cotidiana. 

Sin embargo, como tantas veces se ha visto en la historia del 
conocimiento, el pensamiento creativo, que no busca una utili- 
dad práctica inmediata, resulta ser muy fructífero. Así, esta 
nueva geometría concebida a principios de siglo resulta al día 
de hoy indispensable para la dinámica que estamos estudian- 
do, y que tiene tantas aplicaciones. 

Ahora podemos entender por qué los fractales y los atrac- 
tores extraños se hallan tan íntimamente conectados. Como ha- 
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FIGURA 111.19. 


bíamos visto, un atractor extraño está recorrido por el punto 
que representa al sistema dinámico, que va describiendo una 
curva de infinita complejidad, que se estira y al mismo tiempo 
se pliega y repliega ad infinitum, y cuya sección de Poincaré está 
formada por puntos agrupados con características de autose- 
mejanza. 

Un atractor extraño es, por lo tanto, una curva fractal. 


IV. Caos determinista en los cielos 


PODRÍAMOS IMAGINAR la sorpresa de Newton si leyese el artícu- 
lo de la revista Science del 3 de julio de 1992, donde se confir- 
ma que el comportamiento del sistema solar entero evidencia 
signos de caos. Gracias a las computadoras modernas se ha 
podido encarar el viejo problema de cuán estable es este sis- 
tema, ya planteado en la época del gran matemático y físico 
Henri Poincaré. 
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Una característica del sistema planetario, como los demás 
sistemas que estudia la mecánica celeste, es que para los tiem- 
pos considerados, que pueden ser de hasta miles de millones 
de años, no hay prácticamente disipación de energía por roza- 
miento, ya que los cuerpos celestes se mueven en un vacío casi 
perfecto, y las pérdidas de energía por otros efectos, como el 
de radiación, tampoco son importantes, de'modo que se tie- 
nen aquí cuerpos que forman sistemas dinámicos conservati- 
vos O hamiltonianos. 

Para el sistema solar, los resultados son sorprendentes: G. 
Sussman y J. Wisdom describen en el mencionado artículo el 
cálculo por integración numérica de la evolución de todo el sis- 
tema para los próximos cien millones de años, lo que requirió 
un mes de tiempo de una computadora especial para cada co- 
rrida del programa. 

Resultó que el comportamiento de los nueve planetas a par- 
tir de los próximos cuatro millones de años revela que el sis- 
tema planetario está en estado caótico. 

Para nuestra tranquilidad, esto no significa que el caos en el 
sistema solar sea de tales características que se vaya a aniquilar 
dentro de poco tiempo, con planetas chocando entre sí, o hu- 
yendo hacia otras galaxias, sino que sus órbitas son impredeci- 
bles cuando se calculan para tiempos del orden de los cien mi- 
llones de años y, por lo tanto, sólo se puede anticipar que se 
moverán en el espacio dentro de zonas determinadas. Estas no 
se superponen, de modo que no presagian colisiones entre pla- 
netas, al menos hasta donde se ha calculado, que corresponde 
a los próximos cien millones de revoluciones de nuestro plane- 
ta alrededor del Sol. 

También se encontró que el subsistema formado por Júpiter 
y sus satélites es caótico, lo mismo que el de la órbita del plane- 
ta Plutón. 

Los investigadores estimaron el error con que se conoce la po- 
sición inicial de un planeta en su órbita en sólo 0.00000001%. 
Cabría esperar que si se calculan dos órbitas posibles que difie- 
ran inicialmente en su posición en ese porcentaje, las órbitas 
mantendrían a medida que pasa el tiempo dicha diferencia, 
que daría una distancia entre las mismas del orden de 10 m 
hasta 1 km. 

Pero se encontraron con que la distancia entre esas dos órbi- 
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tas alternativas se multiplicaba en el cálculo por 3 cada 4 millo- 
nes de años y, por lo tanto, pasaba a ser 9 veces mayor a los 8 
millones de años, 27 a los 12 millones de años, y así sucesiva- 
mente; a los 100 millones de años, entonces, podrían diferir 
en la posición del planeta en un 100% para la región caótica, 
también llamada “de resonancia”. : 

En particular, para la Tierra, un error de medición de la ubi- 
cación inicial en su órbita de tan sólo 15 metros, hace imposi- 
ble predecir en qué posición de su órbita estará dentro de cien 
millones de años. 

En realidad el cálculo muestra que estas regiones caóticas, 
llamadas de resonancia, están restringidas a una porción del 
espacio, de modo que en el caso del sistema solar, no eviden- 
cien futuras colisiones entre planetas. 

Pero la importancia de este descubrimiento es algo más que 
una curiosidad astronómica de interés sólo para especialistas. 
Entre obras consecuencias, tiene una eminentemente práctica: 
por las rutas del cielo, no solamente circulan “vehículos pesa- 
dos” como son los planetas y sus satélites, sino que ellas son 
cruzadas también por infinidad de fragmentos tales como as- 
teroides y cometas. Mediante la observación con telescopios y 
el cálculo basado en las ecuaciones de la mecánica celeste, es 
posible prever dónde y cuándo se cruzarán los caminos de es- 
tos objetos. Así se sabe que en julio de 1994 un fragmento, de 
tamaño considerable, chocaría con Júpiter, cayendo sobre la 
cara que en ese momento no será visible desde la Tierra. Tam- 
bién se sabe que el asteroide 1989 AC, que aparentemente está 
en resonancia con Júpiter, cruzará la órbita de la Tierra en el 
año 2004, pasando a una distancia de nuestro planeta de 0.011 
unidades astronómicas, o sea 1650000 km. 

Afortunadamente las ecuaciones permiten definir con pre- 
cisión la ubicación del planeta para dentro de unas pocas re- 
voluciones alrededor del Sol, de modo que la humanidad es- 
pera este evento con total tranquilidad. Pero como hemos 
visto, la capacidad de predecir si habrá o no colisiones se va 
perdiendo en una especie de niebla a medida que avanzamos 
en el futuro. 

¡Qué lejos queda el concepto de las leyes de la mecánica 
celeste que permiten predecir para siempre el movimiento de 
los astros! 
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Pero, ¿no es esto contradictorio con el hecho de que cono- 
cemos el movimiento ordenado de los planetas desde hace mi- 
lenios, y que ya los babilonios podían predecir con exactitud 
un eclipse con años de anticipación? ¿Acaso la revolución new- 
toniana no consistió en el descubrimiento de las leyes inmu- 
tables que gobiernan todos los sistemas dinámicos que apare- 
cen en la naturaleza? 

Ahora comprendemos que ese orden es tal para un tiempo 
de observaciones de unas pocas decenas de miles de vueltas de 
la Tierra alrededor del Sol, pero que se desvanece para escalas 
de tiempo miles de veces mayores. 


Los DESCUBRIMIENTOS DE POINCARÉ 


La existencia de este problema era también sospechada por 
los científicos del siglo pasado, pero fue Henri Poincaré (figu- 
ra 1V.1) el que lo evidenció en su verdadera magnitud. 


FIGURA 1V.1. Henri Poincaré. Tomado de La Recherche, 232, mayo 1 991, 
vol. 22, p. 566. 
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En 1887, el rey Óscar HI de Suecia instituyó un premio de 
2500 coronas para quien produjese una respuesta a una pre- 
gunta fundamental para la astronomía: ¿es estable el sistema 
solar», definiéndose la estabilidad como la situación en que pe- 
queños cambios en los movimientos planetarios sólo traen 
pequenas alteraciones en el comportamiento del sistema. 

En aquella época podía caber, por ejemplo, la sospecha de 
que la Tierra terminaría por precipitarse sobre el Sol. Al tratar 
de resolver este caso, Poincaré abrió el rumbo para tratar los 
problemas de estabilidad en sistemas dinámicos complejos, y si 
bien no pudo resolver el problema para los diez cuerpos que 
forman el sistema solar, recibió de todos modos el premio por 
sus importantes aportes, entre ellos la creación de la. topo- 
logía. 


Por supuesto que las leyes de Newton siguen siendo válidas, pero 
sus soluciones exactas requerían una inteligencia como la de La- 


place que introdujera en las mismas datos con precisión infinita. 


Los astrónomos sólo pueden conocer en forma aproximada 
las condiciones iniciales de velocidades y posiciones de los cuer- 
pos celestes, pero esta precisión limitada a una cierta cantidad 
de cifras decimales no ha sido obstáculo para muchos cálculos, 
ya que se ha trabajado normalmente con ecuaciones en las 
que pequeñas variaciones en las condiciones iniciales produ- 
cen efectos proporcionalmente pequeños, pero, ¿qué ocurre 
cuando un sistema obedece a situaciones de alta sensibilidad a 
las condiciones iniciales? 


TRES O MÁS CUERPOS PUEDEN GENERAR CAOS 


Para un sistema hamiltoniano de sólo dos cuerpos, como la 
Tierra y la Luna o la Tierra y el Sol, las ecuaciones de Newton 
se pueden resolver en forma exacta, este problema se llama in- 
tegrable, y su solución corresponde a una órbita elíptica. Pero 
si se agrega un tercer cuerpo (por ejemplo, al introducir el 
efecto que produce Júpiter, con su masa que es la milésima par- 
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te de la del Sol, sobre el sistema Tierra-Sol), hay que utilizar 
un método aproximado, llamado de perturbaciones. 

En dicho método, esa perturbación del orden de un milési- 
mo que produce la atracción gravitatoria de Júpiter sobre el 
sistema Tierra-Sol, se suma a la solución del caso de los dos 
cuerpos, Tierra y Sol, logrando así una mejor aproximación a 
la realidad. A este resultado se le vuelve a sumar el efecto de 
perturbación de Júpiter, pero elevado al cuadrado, o sea un 
factor que es un millonésimo, y así reiteradamente, en una se- 
rie de aproximaciones donde cada una deberá ser de menor 
magnitud que la anterior. 

Se espera que esta serie formada por la suma de términos 
con valores decrecientes converja, es decir, que si para una 
suma de, por ejemplo, 1000 términos se tiene como resultado 
un cierto número, y al agregar el término 1001 la suma crece 
en 1%, el sumar el término 1002 produzca un nuevo aumento 
de la suma menor que ese 1%, por ejemplo 0.9%, y que cada 
vez que se agregue otro término, la nueva suma aumente en 
forma decreciente, tendiendo a una cifra que sea prácticamen- 
te constante para un número suficientemente grande de tér- 
minos. 

Esto permitiría afirmar que el problema de tres cuerpos está 
resuelto: la ubicación de cada uno de ellos en su órbita estará 
dada por esos números calculados con el método de perturba- 
ciones. Pero al aplicarlo al caso del sistema Sol-Tierra-Luna, 
Poincaré descubrió, para su asombro, que algunas órbitas se 
comportaban caóticamente, es decir que su posición era impo- 
sible de predecir mediante ese cálculo, lo que implica que no 
hay aquí un comportamiento lineal de las ecuaciones. 

La suma de aproximaciones divergía —su resultado era un 
número cada vez más grande— con lo que el efecto de estas 
perturbaciones infinitésimas se iba amplificando y que en cier- 
tas situaciones límite podrían llegar a sacar a un planeta por 
completo de su órbita. 

Pero Poincaré no pudo avanzar más y llegar a resolver el caso 
real del sistema solar completo, por las dificultades que signifi- 
caba intentar el cálculo para diez cuerpos. Como hemos visto, 
este problema tuvo que esperar un siglo, para ser encarado con 
las modernas herramientas de la computadora y el cálculo por 
integración numérica. El gran matemático demostró, sin em- 
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bargo, la posibilidad de que un sistema dinámico totalmente 
determinista, como el de los planetas en órbita, llegue a un 
estado de caos donde no se pueda predecir su comportamien- 
to futuro. Basta para ello que responda a una situación no li- 
neal, en la que una pequeñísima fluctuación se vaya amplifi- 
cando al ser reiterada un gran número de veces. - 

Además, Poincaré probó que el caos puede aparecer ya con 
sistemas relativamente simples, como puede-ser uno formado 
por sólo tres cuerpos, para el que las estructuras que los des- 
criben en el espacio de las fases forman la complicada geo- 
metría a la que hacíamos referencia en el capítulo 2. 


Los FENÓMENOS DE RESONANCIA 


Una contribución importante de Poincaré fue demostrar que 
en este sistema las inestabilidades son debidas al fenómeno de 
resonancia. 

La resonancia aparece cuando hay una relación numérica 
simple entre periodos de rotación de dos cuerpos del sistema 
solar. Por ejemplo, Plutón tiene un periodo de rotación en su 
Órbita alrededor del Sol de 248.54 años y Neptuno de 164.79 
años; la relación entre periodos es entonces de 3 a 2, lo que se 
indica como resonancia órbita-órbita 3:2. Puede haber tam- 
bién resonancias entre el periodo orbital de un objeto y el de 
rotación del mismo alrededor de su eje (espín), como, por ejem- 
plo, la Luna, con una resonancia espín-órbita 1:1 que es la 
responsable de que aquélla muestre siempre la misma cara a 
la Tierra. 

Un efecto que puede resultar de la resonancia es la indefini- 
ción en la posición de un planeta en su órbita. Tal es el caso de 
Plutón, con su resonancia órbita-órbita 3:2 con Neptuno. Se 
ha hecho el cálculo de su posición mediante una computado- 
ra, pasando el programa dos veces, cada una con la posición 
inicial muy ligeramente diferente. Las dos órbitas así obteni- 
das ubican a Plutón en lados opuestos respecto del Sol des- - 
pués de cuatrocientos millones de años. 

Otro efecto posible de la resonancia es producir un brusco 
aumento de la excentricidad de la órbita de un asteroide que 
lo arroje sobre un planeta. Se ha comprobado que dicho efec- 
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to es el responsable de la existencia de zonas vacías o “lagunas 
de Kirkwood” en el cinturón de asteroides que existe entre 
Marte y Júpiter. J. Wisdom, utilizando un método de integra- 
ción numérica para esta dinámica caótica, demostró que el 
efecto de resonancia produce cambios violentos en la excentri- 
cidad de las órbitas, los que acabaron por fragmentar a los as- 
teroides y lanzarlos sobre Marte, y también sobre la Tierra: un 
ejemplo dramático de la ubicuidad del caos en el sistema solar. 

Las resonancias son importantes en los sistemas hamilto- 
nianos y a menudo implican situaciones caóticas. 

Para comprender esto, podemos visualizar el movimiento 
del sistema en el espacio de las fases (véase la figura 1V.2). 


DA 


Zona caótica 


Sección de Poincaré 


FIGURA IV.2. 


Aquí como para el péndulo ideal, podemos considerar valo- 
res pequeños de energía, a los que, en forma similar a las oscila- 
ciones de poca amplitud del péndulo, corresponden los movi- 
mientos llamados “de libración”, que son pequeñas oscilaciones 
alrededor de la posición de equilibrio en reposo; y así como 
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para los valores más grandes de energía, el péndulo efectúa mo- 
vimientos de rotación, para un sistema kepleriano (un cuerpo 
único en órbita alrededor del Sol), el movimiento se denomina 
“de circulación”. Entre ambos movimientos, de libración y de 
circulación, está la separatriz, que para el caso de una pertur- 
bación (resonancia), puede dar origen a una zona caótica. 

Analicemos la sección de Poincaré correspondiente, donde 
el punto equivale a una órbita periódica. Si.se combinan dos 
péndulos, se tiene un sistema dinámico de dos cuerpos. Cada 
círculo concéntrico que rodea al punto central en la sección 
de Poincaré corresponde a un movimiento casi periódico, que 
es una combinación de los dos movimientos circulares, cada 
uno con su propio periodo. En el espacio de las fases se tiene 
un toro, sobre cuya superficie se va arrollando la trayectoria 
casi periódica. Para que haya un movimiento casi periódico es 
indispensable que la relación entre los periodos de rotación 
sea un número irracional. En este caso, al permanecer la curva 
sobre el toro, su intersección con la sección de Poincaré está 
contenida en una curva cerrada (curva invariante). Si se va 
aumentando el periodo menor, se tienen curvas concéntricas 
de radio creciente. 

Pero en aquellas zonas en las que la relación de periodos es 
un número racional, como, por ejemplo, 3:2, 5:3, etcétera 
(condición de resonancia), puede aparecer el caos. Esto se evi- 
dencia en la sección de Poincaré porque aparecen regiones de 
inestabilidad, donde los puntos sucesivos que marca la trayec- 
toria caótica al interceptar dicha sección, llenan toda la región 
inestable de manera aleatoria. 

Dentro de esta zona caótica de resonancia aparecen islas de 
estabilidad, en cada una de las cuales hay una estructura pa- 
reada a la que se encuentra en el centro de la sección de Poin- 
caré. Alrededor de cada isla existe una zona caótica en la que 
el movimiento es inestable: una trayectoria puede girar en las 
curvas dentro de las islas (libración), y luego en las curvas de 
la zona central (circulación). La zona caótica corresponde a las 
separatrices en el caso del péndulo, y en esta zona el movimien- 
to es extremadamente inestable, y muy sensible a las condicio- 
nes iniciales. 

El origen del movimiento caótico de Mercurio, Venus, Tie- 
rra y Marte es la existencia de resonancias entre los periodos 
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de precesión (o sea de rotación en el espacio de la órbita pla- 
netaria) de estos planetas. Los cálculos muestran que nuestro 
planeta está lejos de la zona central estable de movimientos 
casi periódicos, y más bien en la zona caótica cerca de una isla 
(figura anterior), motivo por el cual no se puede predecir su 
posición a los cien millones de años. 


¿COMPORTAMIENTO CAÓTICO EN LAS GALAXIAS? 


Se ha descrito hasta aquí el método para estudiar la evolución 
de un sistema hamiltoniano formado por tres o más cuerpos, 
pero ¿qué ocurriría si se le aplicara a un sistema de una canti- 
dad inmensa de componentes como lo es una galaxia, que 
también es un sistema hamiltoniano, ya que para las escalas de 
tiempo involucradas, de cientos de millones de años, la pérdi- 
da de energía por radiación o fricción es despreciable? 

En 1962, el astrónomo M. Hénon encaró este problema tra- 
tando de calcular el movimiento de estrellas individuales alre- 
dedor del centro de una galaxia. Una de las diferencias con el 
sistema solar es que, en éste, el centro de atracción está concen- 
trado en el Sol, una esfera alrededor de la que giran los plane- 
tas en órbitas planas, mientras que una galaxia se puede mode- 
lar como una zona de atracción con forma de disco alrededor 
de la cual giran las estrellas en órbitas de tres dimensiones. 

Hénon computó con el método de Poincaré la intersección 
de las trayectorias sucesivas de una estrella con un plano, deter- 
minando qué cambios se producían para diferentes energías 
del sistema. Como era de esperarse para un sistema hamilto- 
niano, la sección de Poincaré del toro tridimensional genera- 
do por las trayectorias mostraba curvas concéntricas, de áreas 
encerradas proporcionales a la energía (véase la figura 1V.3). 

Pero llegó un momento, para las energías más altas, en el que 
las curvas se quebraban, desapareciendo para ser reemplazadas 
por puntos distribuidos aparentemente al azar en zonas dentro 
de las cuales aparecían otras curvas como islas en un mar bo- 
rrascoso. En los estados del sistema correspondientes se hace 
entonces imposible establecer con precisión las órbitas estelares. 

En resumen: vemos así como aun la mecánica celeste, que 
era considerada el mejor ejemplo de la capacidad predictiva 
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FIGURA IV.3. 


de las ciencias físicas, muestra los límites impuestos por la exis- 
tencia de los fenómenos no lineales, los que aparecen tanto en 
el sistema solar como en las galaxias. 

Estos hallazgos no sólo son importantes para la mecánica 
celeste. Hay otros sistemas dinámicos de gran interés, como los 
plasmas formados por gases ionizados, o sea tan calientes que 
muchos de sus átomos pierden sus electrones. 

El principal interés de estos sistemas se debe a que se está 
investigando cómo utilizarlos para hacer reactores de fusión 
nuclear que permitan suministrar energía segura y barata. * 
Estudiándolos como sistemas hamiltonianos de una inmensa 
cantidad de componentes, sometidos ya no a la atracción gra- 
vitatoria sino a campos eléctricos y magnéticos, se han hallado 
secciones de Poincaré con islas de regularidad, similares a las 
figuras de Hénon, lo que confirma que pueden pasar a un 
estado de caos para ciertos valores críticos. 


V. Caos, entropía y la muerte del universo 


Como consecuencia de la Revolución Industrial apareció la ne- 
cesidad de comprender el fenómeno de la generación y apro- 
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vechamiento del calor en las máquinas de vapor. Esto involu- 
craba estudiar sistemas de gases y de líquidos de los que sólo se 
pueden medir aquellas propiedades globales como la tempe- 
ratura, presión, volumen, viscosidad, que no requieren el co- 
nocimiento de las posiciones y velocidades de cada uno de sus 
átomos. 


LEYES PARA EL CALOR 


El primer gran logro de la termodinámica fue la ley de conser- 
vación de la energía, que establece que el calor es una más de las 
formas en que se presenta la energía, y que la energía total 
involucrada en un proceso puede cambiar de características, 
pasando por ejemplo de calórica a cinética, eléctrica o química, 
pero jamás se pierde. En una máquina de vapor se transforma 
la energía calórica en energía cinética (de movimiento) de un 
pistón, y si se frotan dos barras, se obtiene calor por fricción. 
Los conceptos de ca- 
lor y temperatura se 
hicieron más precisos 
con el desarrollo de la 
teoría cinética de los 
gases: la energía caló- 
rica de un cuerpo es la 
suma de la energía de 
movimiento, Oo cinética, 
dei totaJ de las molécu- 
las que lo componen, 
mientras que la tempe- 
ratura es proporcional a la energía cinética de la molécula pro- 
medio. Nuestra imagen de un gas es hoy (figura V.1) la de un 
número extremadamente elevado de moléculas que se com- 
portan como diminutas bolas de billar, elásticas, que se mueven 
en línea recta a grandes velocidades sin preferencia por ninguna 
dirección en particular en el espacio, hasta colisionar con otras 
y rebotar vigorosamente. Las moléculas que forman el aire en 
condiciones normales de presión y temperatura, por ejemplo, 
se mueven en promedio a 1800 kilómetros por hora y realizan 
cinco mil millones de colisiones cada segundo. Este agitado 
movimiento no cesa nunca y aumenta al crecer la temperatura. 


FIGURA V.1. 
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Hay que resaltar aquí el notable hecho de que la palabra gas 
fue inventada por el médico holandés Van Helmont en el siglo 
xvII. Hasta ese entonces se hablaba de “vapores” y “espíritus”, 
pero él, con rara visión, consideró que los mismos estaban for- 
mados por partículas invisibles disparadas en todas direccio- 
nes, y que eran la imagen del caos, palabra de la que derivó el 
nombre de “gas”. ' 

Si el gas está contenido en un recipiente, sus moléculas ha- 
rán impacto en su movimiento con las paredes del mismo, ejer- 
ciendo globalmente una presión que será tanto mayor cuanto 
más energía cinética tengan sus moléculas, o sea cuanto más 
caliente esté el gas. 


La segunda ley de la termodinámica establece que, si bien la cantidad 
total de energía es constante en un sistema aislado, que es aquel 
que no recibe energía ni materia desde su entorno, la energía útil 
capaz de ser utilizada para efectuar un trabajo dismuniye, pues en 


todo proceso una fracción de la energía se transforma inevitable- 
mente, por fricción, rozamiento, etcétera, en calor, la que no puede 
ya aprovecharse para su conversión en alguna forma de energía. 


La termodinámica nos está dando, entonces, dos noticias: 
una “buena” y otra “mala”. La “buena” es que la energía es ina- 
gotable, que siempre habrá energía en el Universo; la “mala” 
es que la energía viene en dos tipos, de los que sólo uno de 
ellos nos es útil, y que, además, esta energía útil está disminu- 
yendo y algún día desaparecerá. 

¿Qué ejemplo tenemos de energía útil capaz de producir 
trabajo en un sistema aislado? Un caso simple es el de un re- 
cipiente cerrado, dividido en dos compartimentos mediante 
un tabique móvil que actúe de pistón (véase la figura V.2). 

Supongamos que inicialmente llenamos un compartimento 
con un gas caliente, y el otro con el mismo gas a menor tem- 
peratura. La energía del gas es entonces utilizable, pues el gas 
caliente tiene sus moléculas con mayor energía de movimiento 
que el gas frío y, por lo tanto, ejerce mayor presión sobre el 
pistón, empujándolo. 

A este estado organizado, donde hay una diferencia apro- 
vechable entre las dos regiones del sistema, lo consideramos 
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Gas trío Gas caliente 


FIGURA V.2. 


como ordenado, en contraposición con el caso en que el reci- 
piente está lleno uniformemente con gas a la misma tempera- 
tura, el que, por lo tanto, ejerce de ambos lados del pistón igual 
presión y no lo puede mover. 


LA FLECHA DEL TIEMPO 


Si bien las moléculas del gas tienen individualmente movimien- 
tos reversibles —una película cinematográfica de sus trayecto- 
rias proyectada del principio al fin o del fin al principio sería 
igualmente válida, pues las leyes de Newton serían las mismas— 
en conjunto el gas sufre un proceso irreversible: si inicialmente 
había una diferencia de temperatura entre ambos comparti- 
mentos, ésta terminará por desaparecer, porque las moléculas 
del gas más caliente irán perdiendo parte de su energía al em- 
pujar el pistón y, por lo tanto, su temperatura disminuirá. 

Al mismo tiempo, las moléculas del gas frío colisionan con 
el pistón que avanza hacia ellas, con lo que reciben más ener- 
gía cinética, de modo que aumenta su temperatura. Después 
de un tiempo, ambos compartimentos contendrán un gas ti- 
bio, que ejerce igual presión de ambos lados del pistón y que, 
en consecuencia, no puede producir trabajo. 

Nuestra experiencia nos dice que nunca se ha visto que es- 
pontáneamente aparezca una diferencia de temperatura, y 
vuelva por sí sólo a calentarse el gas en un compartimento y en- 
friarse en el otro. 
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No es que haya una prohibición que surja de las leyes de la 
física, pues para que esto ocurriera bastaría con que la mayoría 
de las moléculas, que individualmente se están moviendo en 
cualquiera de las direcciones posibles sin preferencia por nin- 
guna, de manera espontánea y por una simple coincidencia se 
echaran a mover todas en una misma dirección, hacia el pistón 
las de un compartimento y alejándose del pistón las del otro. 

Se ha calculado la probabilidad de que algo así pase y el re- 
sultado es que se requeriría esperar toda la edad del Universo, 
y mucho más, para que se produjera tal fenómeno. Por lo tan- 
to, esta vuelta al estado inicial que equivale a una reversibili- 
dad en el tiempo, en la práctica es tan improbable que a 
mos considerar que jamás ocurrirá. 

De este modo se introdujo en la física la noción de procesos 
no reversibles en el tiempo, una asimetría que se ha denomi- 
nado “flecha del tiempo”, pues los procesos termodinámicos 
para sistemas aislados sólo ocurren en una dirección: aquella 
en que con el tiempo crece la energía no utilizable. 

Más aún, si se considera todo el Universo como un sistema 
termodinámico aislado, ya que por definición no hay un exte- 
rior al mismo desde el que pueda llegar materia o energía, se 
llega a la conclusión de que finalizará con una “muerte térmi- 
ca”, cuando toda la energía que contiene se degrade a energía 
calórica no utilizable, hasta alcanzar un estado final de equi- 
librio. 

Este estado al que se llegaría inexorablemente se identifica 
con el desorden o caos, porque se visualiza como el fin del 
Universo tal como lo conocemos, de la organización que man- 
tiene toda la actividad ordenada que percibimos, y en el cual, 
desde las galaxias y sistemas solares hasta los seres que las habi- 
tan, desaparecerían para transformarse en un caos homogé- 
neo de átomos moviéndose ciegamente para el resto de la eter- 
nidad. 

Así, la formulación más general de la segunda ley de la ter- 
modinámica expresa que todo sistema físico aislado, incluyen- 
do el Universo, aumenta espontáneamente su desorden. 
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LA ENTROPÍA 


En 1865, el físico Clausius introdujo el concepto de entropía 
para plantear en una función matemática precisa esta tenden- 
cia de la evolución de los sistemas termodinámicos. 

La función entropía aumenta en un sistema aislado de la 
misma manera que el desorden, y se la considera como una 
medida de ese desorden. 

Clausius reformuló, además, las dos leyes de la termodiná- 
mica del siguiente modo: 


La energía del Universo es constante. 
La entropía del Universo tiende hacia un máximo. 


Estos conceptos fueron analizados mediante la mecánica es- 
tadística por el físico Ludwig Boltzmann, quien demostró que 
el estado final de un sistema aislado, cuando no hay un cambio 
con el tiempo en las propiedades macroscópicas, tales como la 
densidad, presión, temperatura, etcétera, es el de equilibrio 
termodinámico, hallándose, en ese caso, el conjunto de sus 
moléculas en un estado de máxima entropía. 

En 1875 planteó su definición de entropía, como propor- 
cional al número de posibles configuraciones distintas del sis- 
tema que sean compatibles con las propiedades macroscópi- 
cas. Así, para un gas aislado en un recipiente y en equilibrio, 
hay una cierta cantidad de configuraciones, entendiéndose 
por configuración a cada una de las posibles maneras en que 
se pueden distribuir las moléculas en el recipiente, y que pro- 
“«duzcan como resultado idénticas propiedades macroscópicas 
de presión, temperatura, densidad, etcétera. Todas estas con- 
figuraciones tienen la misma probabilidad de presentarse, de 
modo que basta con calcular su número total para obtener la 
entropía del sistema. 

Cuanto mayor sea este número, mayor es la entropía que se 
puede dar, entonces, como una cantidad numéfica para un sis- 
tema en equilibrio, lo que constituye una gran mejora respec- 
to de utilizar simplemente los conceptos de “orden” y *desor- 
den”, que son vagos. : 
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Se podría considerar que el desorden o caos que reina en 
un gas en equilibrio se debe al gran número de moléculas. La 
existencia de caos determinista en sistemas con pocos grados 
de libertad demuestra que lo anterior no es obvio en absoluto. 
Se lo asocia más bien con la enorme diferencia entre los volú- 
menes ocupados por estados ordenados y caóticos en el espa- 
cio de las fases. 

Recordemos que éste es un espacio matemático, cuyo núme- 
ro de dimensiones depende del número de variables indepen- 
dientes o grados de libertad del sistema dinámico a representar. 

Si el sistema está formado por n partículas que no estén liga- 
das entre sí, habrá 3n coordenadas de posición y 3n coorde- 
nadas de impulso, lo que aun para unos pocos centímetros 
cúbicos de cualquier gas implica un número inmenso de di- 
mensiones del espacio de las fases. Así, por ejemplo, para un 
centímetro cúbico de aire, que contiene unos 2.7 x 101? molé- 
culas, la cantidad de dimensiones del espacio de las fases es de 
alrededor de 6 x 2.7 x 101% o 160000 000 000000 000000. 

Obviamente carece de sentido intentar visualizar tal “espa- 
cio”. De todos modos se le puede esquematizar como si fuera 
tridimensional. Lo que interesa aquí es que si en un instante 
de tiempo dado, esa porción de aire tiene, por ejemplo, las si- 


guientes propiedades macroscópicas observables: una tempe-. 


ratura de 20 grados centígrados, una presión normal de 1013 
milibares, un volumen de 1 cmn?; quedará representada dentro 
de una cierta región que designaremos como HH, en el espacio de 
las fases (véase el capítulo III) y no necesitaremos preocupar- 
nos por conocer el infinito detalle de todas las posiciones e im- 
pulsos individuales de las moléculas componentes. 

Estos quedan representados en cada uno de los puntos den- 
tro de la región. Á cada punto le corresponde un instante en 
el tiempo con una determinada configuración, una distribu- 
ción en el espacio de las moléculas individuales del gas, cada 
una con su velocidad. 

Un punto diferente simboliza otra distribución de posicio- 
nes y velocidades, pero si ambos puntos están contenidos en la 
misma región las propiedades macroscópicas son las mismas. 

Se puede dividir así el espacio de las fases en un número de 
regiones, a cada una de las cuales le corresponde un conjunto 
de características macroscópicas observables diferentes. 
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FIGURA V.3. 


El tamaño de una región, que está dado por la cantidad de 
puntos que la componen, nos indica la cantidad de posibles 
distribuciones diferentes de las moléculas. 

Y cómo evoluciona el gas con el tiempo se representa global- 
mente con la trayectoria que, saliendo del punto que indica la 
distribución de las posiciones y velocidades de las moléculas 
para el instante inicial, se desplaza por el espacio de las fases, y 
va emigrando de esa región inicial a otras a medida que cam- 
bian las características macroscópicas. En la figura V.3 se indi- 
can esas regiones por sus propiedades macroscópicas de pre- 
sión Py temperatura 7. 

En un proceso termodinámico los tamaños de las distintas re- 
giones pueden ser muy diferentes. Para un gas ideal, que esté 
comprimido en un cierto volumen dentro de una caja aislada y 
pase a expandirse hasta ocupar todo el volumen, se podrían 
estimar los tamaños relativos de las dos regiones correspon- 
dientes: la inicial y la final de equilibrio, designadas con Po, T, 
y P, T, en la figura V.3. Al llegar a esta última región, el gas es- 
tará en equilibrio térmico, el que se caracteriza porque sus mo- 
léculas están distribuidas de una manera uniforme en el vo- 
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lumen de la caja y se mueven en todas direcciones con un ran- 
go de velocidades conocido como “distribución de Maxwell”. 

Para poder tener una idea de los factores de magnitud en 
Juego veamos un ejemplo muy simplificado: supongamos que 
hay un total de seis moléculas dentro de una caja dividida en 
dos partes. Hay varias posibilidades para su distribución en la 
caja: dos en una mitad y cuatro en la otra, cinco en una mitad 
y una en la segunda, etcétera. Las moléculas son idénticas, en 
su movimiento no distinguen diferencias entre las dos mitades, 
y se puede calcular fácilmente cuántas configuraciones posi- 
bles hay para cada distribución. 

Resultan veinte combinaciones diferentes que dan una distri- 
bución uniforme de tres moléculas en cada mitad, contra seis 
donde las seis están en una de las dos mitades. Pero a medida 
que aumenta la cantidad de moléculas, la diferencia aumenta 
muy rápidamente: para 10 moléculas son 252, contra 10; para 
20 moléculas son 1.3 x 10%%, contra 20; para 100 moléculas son 
10% configuraciones que dan distribuciones uniformes contra 
sólo 100 en las que las 100 están concentradas en una mitad. 

Téngase ahora en cuenta la cantidad de moléculas que hay 
en un litro de cualquier gas, y se comprenderá que, en conse- 
cuencia, es inmensa la cantidad de posibles configuraciones 
con una distribución uniforme en la caja, en comparación con 
aquéllas donde se hallan todas en una mitad. 

Como a cada distribución posible le corresponde un punto 
en el espacio de las fases, se deduce que la región correspon- 
diente al gas en equilibrio térmico es, de lejos, la mayor en el 
espacio de las fases. 

Supongamos ahora que partimos de la situación en que el 
aire está todo acumulado en una porción de la caja. Inmedia- 
tamente comenzará a expandirse ocupando cada vez más espa- 
cio hasta llenarla; después de un tiempo llegará al equilibrio 
térmico, con una distribución uniforme de temperatura y pre- 
sión en toda la caja. 


LA ENTROPÍA Y EL ESPACIO DE LAS FASES 
¿Cómo representaremos este proceso en el espacio de las fases? 


El punto parte de una región muy pequeña —la región que 
representa la colección de los posibles estados iniciales para 
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los que el gas está acumulado en un rincón de la caja—. Cuan- 
do el gas comienza a expandirse, la trayectoria del punto móvil 
entrará en una región del espacio de las fases de mayor tama- 
ño, y luego a regiones más y más grandes mientras el gas siga 
expandiéndose, hasta que finalmente entra en la región de 
volumen más grande del espacio de las fases —la que corres- 
ponde al equilibrio térmico— que, por lo que hemos visto, ocu- 
pa prácticamente todo el espacio de las fases (figura V.3). 

Por lo tanto, las oportunidades de que la trayectoria, después 
de haber abandonado ese pequeñísimo volumen de espacio de 
las fases para entrar al vasto dominio de la región de equilibrio 
térmico, vuelva a la región inicial son prácticamente nulas. 

No es que haya una prohibición dictada por las leyes de la 
naturaleza: como hemos visto, todos los estados son posibles, 
pero la probabilidad de que la trayectoria encuentre nuevamen- 
te esa región es mucho menor que la de encontrar por casuali- 
dad una aguja en un pajar del tamaño del planeta. 

Es fácil ver que no alcanzaría toda la edad del Universo para 
que espontáneamente volviese a aparecer justamente esa con- 
figuración con el gas acumulado en una porción de la caja. Lo 
que no deja de ser tranquilizador, pues la física nos asegura 
así que, si bien no es imposible, la probabilidad de que, de 
pronto, todo el aire de la habitación en que estamos se retire a 
un rincón, dejándonos sin poder respirar, es tan escasa que no 
ocurre ni siquiera una vez en muchos cientos de miles de mi- 
llones de años, y esto en la práctica, para nosotros equivale a 
una imposibilidad. 

Por todo lo anterior deducimos que, una vez que un gas ha lle- 
gado al estado de equilibrio térmico, no sale espontáneamente 
de él, pues aunque todos los otros estados son también posibles, 
éste es mucho más probable por un factor inmensamente mayor. 

Además, si consideramos que la entropía del sistema es una 
medida del volumen de la correspondiente región del espacio 
de las fases, llegamos a la conclusión de que si la trayectoria 
que representa al gas parte inicialmente de la región de volu- 
men muy pequeño, o sea con una entropía mínima, a medida 
que transcurra el tiempo pasará por regiones del espacio de 
las fases con volúmenes o entropías crecientes, hasta llegar al 
máximo de entropía en el equilibrio térmico. 

Se llega así a la formulación de la segunda ley de la termo- 
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dinámica utilizando el concepto de la representación del sis- 
tema en el espacio de las fases. La entropía crece porque los es- 
tados fuera de equilibrio son muchísimo menos probables que 
los de equilibrio (ocupan volúmenes mucho menores en el 
espacio de las fases). Por lo tanto, una vez que el sistema en su 
evolución a lo largo del tiempo llega en el espacio de las fases a 
esa vasta región, es muy poco probable, aunque no imposible, 
que vuelva a la que corresponde a estados fuera de equilibrio. 


VI. El comportamiento de los sistemas 
con gran número de componentes 


EXAMINEMOS qué ocurre con sistemas formados por una in- 
mensa cantidad de moléculas componentes como son los flui- 
dos (líquidos, vapores, etcétera), cuando pasan bruscamente 
de un comportamiento ordenado a 
una confusión caótica. 

No se trata aquí de buscar algún 
fenómeno exótico, que sólo se pueda 
observar tras las paredes de sofistica- 
dos laboratorios de física; en realidad, 
como veremos, esto es totalmente 
normal en nuestra vida cotidiana. 

Observemos el humo que asciende 
por el aire desde un cigarrillo puesto 
en el cenicero o una vela que acaba 
de apagarse (figura VI.1). Asciende 
verticalmente durante varios centíme- 
tros, en una columna ordenada, casi 
rectilínea (flujo laminar), hasta que 
abruptamente aparecen torbellinos 
cada vez más complejos y termina en 
una nube totalmente desordenada que 
acaba por perderse en el aire. 

Los millones de partículas micros- 
cópicas de ceniza caliente que forman FIGURA VL1. 
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el humo, se mueven todos con la misma velocidad durante el 
ascenso en flujo laminar, como si se tratase de una autopista 
con columnas de autos manejados por conductores respetuo- 
sos de las leyes de tránsito. Tras la fase de remolinos, desembo- 
can en la nube desordenada, y en ella las partículas de ceniza 
están animadas de movimientos en todas direcciones, de modo 
tal que si se quisiera medir la velocidad del humo en un punto 
de esa nube, se comprobaría que va variando de instante en 
instante en forma totalmente aleatoria, tan imprevisible como 
si los datos de las mediciones se obtuvieran jugando a la ruleta, 


TURBULENCIA 


Esta condición se denomina turbulenta; en física, un flujo se 
llama turbulento si la velocidad del fluido parece variar al azar 
tanto en el tiempo como en el espacio. 

La turbulencia aparece con mucha frecuencia en los fluidos 
que encontramos en la naturaleza: en corrientes de aire, torren- 
tes de agua, los procesos atmosféricos, corrientes oceánicas, 
procesos en la atmósfera de planetas tales como Júpiter, etcé- 
tera. Como es sabido constituye además uno de los grandes 
problemas de la tecnología moderna (en la industria aeronáu- 


tica, en las cañerías de transporte de petróleo o de agua), y a : 


pesar de los esfuerzos realizados por muchos científicos, se 
está aún lejos de dominar sus principios fundamentales. 

No obstante, algunas de las vías promisorias pasan por estu- 
diar qué ocurre en un fluido cuando se produce la transición 
del orden laminar al caos de la turbulencia. Hasta 1970, se 
aceptaba la teoría del físico Lev Landau, quien entendía que 
en un sistema formado por tantas partículas, cuando se sale del 
flujo laminar y se inician los primeros remolinos, esto es equiva- 
lente al comienzo de las oscilaciones en un péndulo, es decir, 
aparece un movimiento periódico. Los remolinos iniciales se 
dividen muy pronto en remolinos más pequeños, lo que impli- 
ca que aparecen inestabilidades en el fluido que lo hacen osci- 
lar con otro periodo diferente adicional, y esto, a su vez, produce 
remolinos menores con nuevas oscilaciones. Así, la turbulen- 
cia sería inicialmente una superposición de tres o cuatro movi- 
mientos periódicos. Una vez que se desarrolla totalmente, el 
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número de oscilaciones diferentes sería infinito, como si se tra- 
tase de infinitos péndulos que se perturban mutuamente, y eso 
explicaría el comportamiento caótico de la velocidad del fluido. 
Pero en 1970, los matemáticos D. Ruelle y F. Takens propu- 
sieron una interpretación distinta a la de Landau. Coincidían 
en la primera etapa del paso de flujo laminar a los primeros ' 
remolinos. Durante la etapa de flujo laminar, todas las partícu- 
las se mueven con prácticamente la misma velocidad en una 
misma dirección, de modo que la representación del proceso 
en el espacio de las fases es muy simple: un punto, que atrae 
las trayectorias de aquellas que se desvían por pequeñas per- 
turbaciones (véase la fi- 
gura VI.2a). Aquí el atrac- 
tor puntual representa la A 
velocidad constante del 
fluido. 
Pero en realidad la ve- Atractor 
locidad va aumentando, e 
debido a la fuerza ascen- E 
dente del humo caliente, A 
y cuando supera un cier- 
to valor, se produce un 
cambio abrupto: ante la 
más pequeña perturba- 
ción, como puede ser una 
muy leve corriente de B 
aire, las hileras o láminas 
de humo se desvían for- Atractor 
mando rulos que giran 
sobre sí mismos. Aparece 
aquí un movimiento rít- 
mico, con un cierto pe- 
riodo, que se puede re- 
presentar en el espacio 
de las fases (véase la figu- 
ra VI.2b) de la misma ma- 
nera que para el caso de 
un péndulo, con un ci- 
clo límite. 
Este lazo es relativamen- FIGURAS VI.2.a y b. 
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te estable, de modo que constituye un atractor. Pero al seguir 
moviéndose el remolino pronto sufre el efecto de alguna otra 
corriente de aire, que simultáneamente lo hace oscilar en otra di- 
rección, con lo que tenemos un caso similar al del sistema for- 
mado por dos péndulos de distinto periodo del capítulo III. 
Vimos que la correspondiente representación en el espacio de 
las fases es el solenoide en la superficie de un toro tridimen- 
sional, de modo que nos encontramos con tres variables in- 
dependientes o grados de libertad. ¿Qué ocurre de aquí en 
adelante? 

Es en este punto donde se separan ambas teorías. Para Lan- 
dau, una nueva perturbación introduce otra oscilación, que 
nos llevaría a una representación en un espacio de las fases de 
cuatro dimensiones, y así sucesivamente hasta llegar al caos 
con un espacio de las fases de un elevado número de dimen- 
siones, que corresponden al elevado número de variables inde- 
pendientes del fluido turbulento. Para Ruelle-Takens, el caos 
se produce abruptamente y mucho antes: cuando se perturba 
al remolino que tiene dos frecuencias de oscilación simul- 
táneas. 

Lo anterior se visualiza no con una nueva figura en un espa- 
cio de las fases de cuatro dimensiones, sino alterando de raíz 


la forma del atractor que estaba sobre el toro (véase la figura . 


111.13 de la página 71). Es como si el solenoide que formaba 
esa trayectoria hubiera estallado, produciendo una figura tan 
extravagante que Ruelle la llamó “atractor extraño”. Este nom- 
bre fue desde entonces el que adoptaron los científicos para 
denominar a los atractores que indican un comportamiento 
que no se puede predecir a largo plazo. 

Una sección de Poincaré del atractor extraño muestra una 
estructura fractal, de modo que el atractor que era bidimen- 
sional saltó a una dimensión mayor que 2 pero menor que 3. 

Mediciones realizadas en laboratorios han confirmado la 
existencia de este atractor extraño, lo que señala una vía para 
formular las leyes de la turbulencia. 

Otro aspecto importante de este hallazgo es que evidencia 
que también un sistema dinámico formado por un gran núme- 
ro de elementos puede llegar al comportamiento caótico a 
partir de sólo tres grados de libertad. Hasta entonces se opina- 
ba que, como parece indicar el sentido común, sólo puede ha- 
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ber caos cuando el sistema tiene una gran cantidad de varia- 
bles independientes. 


CONVECCIÓN TÉRMICA 


Un notable ejemplo de un fluido dando origen a un fenóme- 
no de ordenamiento espontáneo se puede encontrar en la con- 
vección térmica, o sea el transporte de calor que efectúa un 
fluido caliente al desplazarse hacia una zona más fría. 

Los primeros estudios los realizó H. Benard en 1900. El físi- 
co francés descubrió que si se calienta una capa delgada de un 
fluido (aceite), por debajo, éste puede organizarse espontá- 
neamente en células de convección de un tamaño característi- 
co (véase la figura VI.3), similares a un panal de abejas. En 
cada una de las células de Benard sube fluido caliente por el 
centro y baja fluido frío por los bordes. 

Hoy sabemos que este fenómeno de organización espon- 
tánea está muy difundido en la naturaleza: la superficie del Sol 
está cubierta de células de convección cada una de las cuales 
tiene un tamaño del orden de los mil kilómetros, y en general 
las mismas células se pueden observar en todos los fluidos en 
los que la convección térmica produce movimientos, cuando 
forman capas mucho más anchas que altas. 

Esto vale también para la circulación del aire en la at- 
mósfera y de las corrientes 
oceánicas, las que determi- 
nan en buena parte el clima 
a corto y mediano plazo. 

Supongamos que pone- 
mos agua o aceite en una 
sartén bien grande y calen- 
tamos uniformemente, de 
modo que el líquido en con- 
tacto con el metal tenga la 
misma temperatura en to- 
das sus partes; al calentarse 
conducirá el calor hacia 
arriba, donde se disipará 
en la superficie del líquido. 
Tenemos así una conduc- 
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FIGURA VÍA. 


ción térmica del calor. Pero si vamos aumentando la tempe- 
ratura del calentamiento por encima de un cierto valor apa- 
rece súbitamente otro fenómeno, el de convección térmica 
con la formación de las células. 

Lo que ocurre es que las capas inferiores de fluido, más ca- 
lientes, se van dilatando, por consiguiente tienen menor densi- 
dad que las capas más frías que están encima, y al ser más livia- 
nas ascienden, siendo reemplazadas por esos volúmenes más 
fríos (figura VI 4). Éstos se calientan al ir acercándose al fon- 
do, y simultáneamente los que ascendieron se enfrían, de 
modo que se produce un movimiento circular. 

Es notable la complejidad de estos movimientos: una sección 
del fluido muestra la forma en que están coordinados los senti- 
dos de circulación de las células adyacentes, alternadamente 
en el sentido de las agujas del reloj, o en el sentido contrario. 

La correlación de los movimientos debidos al proceso de 
convección térmica es tanto más notable si se tiene en cuenta 
que el tamaño de estas células de Benard puede llegar a algunos 
milímetros, mientras que el rango de acción de las fuerzas entre 
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moléculas es del orden de un diez millonésimo de milímetro 
(1077 mm). En cada célula hay alrededor de 10% moléculas. 

Debe destacarse que el experimento es perfectamente re- 
producible, o sea que repitiendo las mismas condiciones de ca- 
lentamiento, se llega, a partir del mismo umbral de temperatu- 
ra, a la misma formación de células con similar geometría y 
velocidad de rotación. Pero lo que no se puede controlar es el 
sentido de rotación: se sabe de antemano que con condiciones 
de la experiencia bien controladas, en determinada zona de la 
capa de fluido aparecerá una célula, pero es imposible prede- 
cir si rotará en un sentido o en el opuesto. 

Repetidas pruebas han demostrado que la probabilidad de 
que la célula rote en un sentido o el otro es la misma para am- 
bos casos, es decir que en un punto cualquiera del fluido, la 
velocidad tanto puede estar dirigida en un sentido como en el 
contrario. Esta propiedad se 
puede representar en una 
gráfica como la de la figura 
VI.5, donde la velocidad V 
está indicada en el eje verti- 
cal, con V* para velocida- 
des de rotación en un senti- 
do y V” para el opuesto; en 
el eje horizontal, la varia- 
ble R representa las condi- 
ciones de viscosidad, densi- 
dad, espesor, y la diferencia 
de temperatura entre la su- 
perficie inferior y la supe- 
rior de la capa. Se comienza 

Eigura Ma. a calentar y R pasa de cero a 

valores crecientes, hasta que, 

a partir del valor crítico R, aparece una bifurcación de la cur- 

va: dos ramas iguales V* y V”, siendo imposible determinar de 

antemano por cuál de ellas transitará el comportamiento del 

sistema dinámico, que evidentemente responde a condiciones 
no lineales, con sensibilidad a las condiciones iniciales. 

Si se sigue aumentando la temperatura de calentamiento, 
súbitamente se borran las células, y comienza un movimiento 
aleatorio que es la iniciación de la turbulencia. 
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En 1987, los investigadores M. Dubois, P. Atten y P. Bergé 
hicieron mediciones en el fluido cuando aparecen las células 
de Benard y representaron las correspondientes oscilaciones 
térmicas en el espacio de las fases. Obtuvieron así un atractor 
extraño cuya sección de Poincaré tiene una dimensión entre 
2.8 y 3.1, según las condiciones de calentamiento, lo que indi- 
ca que el comportamiento del sistema se puede describir con 
un mínimo de cuatro grados de libertad. 


LA CONVECCIÓN EN LA ATMÓSFERA 
Y LAS PREDICCIONES METEOROLÓGICAS 


Como sabemos, cotidianamente los meteorólogos publican sus 
predicciones para el clima de los próximos días, con base en 
los datos de mediciones en la atmósfera, en la superficie te- 
rrestre y los océanos, y a observaciones satelitales. Poderosas 
computadoras procesan esos datos en centros de Estados Uni- 
dos y Europa, mediante modelos de la atmósfera como un sis- 
tema dinámico con más de un millón de variables indepen- 
dientes. Hay planes para aumentar cada vez más la potencia de 
esos centros de procesamiento de datos y también la cantidad 
y frecuencia de las mediciones en todo el planeta. Sin embar- 
go, ningún meteorólogo serio cree que en algún futuro cer- 
cano se pueda afirmar en un pronóstico meteorológico algo 
así como “el viernes 26, o sea dentro de 12 días, en el Gran 
Buenos Aires las temperaturas máxima y mínima serán de 19 y 
14 grados, con un 67% de humedad y un 53% de probabilidad 
de lluvias”. 

Actualmente, con todos los instrumentos modernos, el pro- 
nóstico empieza a ser incierto a partir de cuatro o cinco días. 
Ha habido casos (como en Inglaterra, 1987) en que no se pudo 
prever la aparición de un desastroso huracán 24 horas antes de 
que se presentase. 

Se estima que en el futuro se podrá pronosticar hasta con 14 
o 15 días de anticipación, y que más allá sólo habrá la cre- 
ciente niebla de lo incierto. Pero nos estamos refiriendo a la 
predicción exacta para puntos individuales en el planeta; si, en 
cambio, consideramos el clima para toda una región, éste sí es 
predecible, con bastante más precisión que “en junio en la 
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provincia de Buenos Aires la temperatura oscilará entre 18 y 5 
grados”. 

Al ser la atmósfera un fluido, se le estudia como un sistema 
dinámico formado por una inmensa cantidad de compo- 
nentes, sujeto a turbulencias y convecciones. El Sol calienta la 
superficie de nuestro planeta, que a su vez eleva la temperatu- 
ra del aire: en forma similar a lo que vimos cuando exami- 
namos las células de Benard, al calentarse desde abajo la capa 
de aire que forma la atmósfera, aparecen convecciones. 

La distancia entre el ecuador y los polos es de diez mil kiló- 
metros, mientras que el espesor de la atmósfera es de unos diez 
kilómetros hasta la troposfera. En esta capa esférica de fluido 
en rotación aparecen células de convección distribuidas a lo 
largo de seis anillos que rodean el planeta, tres en el hemisfe- 
rio Norte, y otros tres en el Sur, como se indica en la figu- 
ra VL6. Las flechas gruesas indican el sentido general de la ro- 
tación, y las finas la circulación adicional producida por la 
rotación del planeta (efecto Coriolis). 

Hay además corrientes de aire de muy alta velocidad que cir- 
culan en la estratosfera, y que influyen decisivamente sobre el 
clima a largo plazo. Son células inestables, móviles, con límites 
y tamaños cambian- 
tes. Sus característi- 
cas dinámicas pueden 
variar en magnitudes 
que se amplifican al 
doble aproximada- 
mente cada par de 
días, y a su vez actúan 
sobre las corrientes 
de la estratosfera que 
modifican el clima 
desde las regiones 
templadas hasta las 
polares. 

El comportamien- 
to de la atmósfera 
constituye entonces 
un proceso no lineal, 
con realimentación, Figura VL6. 
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altamente sensible a las condiciones iniciales, de modo que 
cumple las condiciones para que sea un sistema caótico. Es 
decir que no importa cuán complejos sean los modelos di- 
námicos, o cuán precisos sean los datos de las mediciones en 
tierra, agua y aire, las leyes de la física imponen un límite más 
allá del cual es imposible hacer predicciones meteorológicas. 

Pero si bien no se puede predecir a más de cuatro o cinco 
días para un punto geográfico específico, ¿podría hacerse una 
predicción a más largo plazo, que diera la tendencia global 
para toda una región que tenga dimensiones relativamente 
grandes? 

El meteorólogo Edward Lorenz se planteó este problema a 
comienzos de la década de los sesenta. En el Massachusetts 
Institute of Technology, donde investigaba, se estaba desarro- 
llando el estudio de dinámicas no lineales, y decidió hacer un 
modelo atmosférico simple, de tres ecuaciones no lineales con 
sólo tres variables independientes en lugar de la inmensa can- 
tidad que puede haber para ese sistema. 

Representó la evolución del clima en un espacio de las fases 
de tres dimensiones y comprobó, para su sorpresa, que las co- 
rrespondientes trayectorias daban origen a un atractor de una 
forma muy curiosa, con dos lazos semejantes a alas de mari- 
posa (véase la figura VI.7). Había descubierto el primer atrac- 
tor extraño. En este espacio de las fases, cada una de las alas 
del atractor representa un posible estado de la atmósfera, por 
ejemplo tiempo lluvioso en el ala izquierda, y tiempo seco y 
estable en el de la derecha. Si las condiciones iniciales son las 
que marca el punto 1 de la izquierda, la evolución seguirá la 
trayectoria que se mantiene en la misma ala: el tiempo será llu- 
vioso. Pero una pequeña perturbación, que cambie las condi- 
ciones iniciales llevando a la atmósfera a la situación represen- 
tada por el punto 2, nos lleva a la trayectoria 2, que evoluciona 
hacia el ala derecha, y el tiempo será entonces seco y estable. 

El atractor de Lorenz tiene estructura fractal, con dimen- 
sión D= 2.06. 

Lorenz publicó este hallazgo en 1963, siendo el primer ejem- 
plo de un modelo calculado con sólo tres variables indepen- 
dientes donde aparecía un comportamiento impredecible a 
largo plazo. 
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FIGURA VI.7. 


EL EFECTO MARIPOSA 


Lorenz acuñó su famosa expresión “efecto mariposa” como 
ejemplo de esta sensibilidad extrema a las condiciones ini- 
ciales: el batir de alas de una mariposa en el Amazonas hoy, 
podría producir una pequeñísima alteración en el estado de la 
atmósfera, que si se amplificase duplicándose cada par de días, 
iría haciendo diverger cada vez más lo que es respecto de lo 
que hubiera sido sin la mariposa, de modo que varias semanas 
más tarde aparecería un ciclón en el Caribe, que, de no haber 
existido el insecto en cuestión, nunca hubiera surgido. 

El modelo de Lorenz con apenas tres variables describe sólo 
cualitativamente la forma caótica en que evoluciona el clima, 
pero lo hace de manera muy similar a los modelos de muchas 
variables independientes. 

En la actualidad, los modelos de predicción meteorológica 
tienen alrededor de un millón de grados de libertad, y esto per- 
mite hacer pronósticos generales para regiones globales con 
una anticipación de hasta un mes. 


104 


VII. Los sistemas lejos del equilibrio 


EL SURGIMIENTO de la termodinámica fue un desafío para los 
físicos del siglo xtx, cuya escuela se basaba en los conceptos de 


Laplace, Lagrange, Hamilton y sus discípulos. Aquí aparecían. 


para el estudio del calor ecuaciones tan universales como las 
de Newton, y que se referían al comportamiento global de la 
materia. Como hemos visto, para aplicar las leyes de Newton a 
los sistemas dinámicos, es preciso definir las posiciones y ve- 
locidades de cada uno de sus elementos componentes. Pero 
un gas, un líquido o un sólido tienen una cantidad inmensa de 
componentes (del orden de 10% moléculas en un centímetro 
cúbico en un gas), de modo que se ignoraba por qué vía se po- 
dría establecer su comportamiento global partiendo del cálcu- 
lo de los movimientos de sus moléculas. 

Con la termodinámica surge un nuevo enfoque, donde las 
ecuaciones del calor, que también son universales, utilizan co- 
mo parámetros propiedades colectivas, macroscópicas, tales 
como la presión y la temperatura, y no parecen requerir del 
conocimiento detallado de lo que ocurre con cada una de las 
moléculas participantes. 

La termodinámica aparece entonces como una herramienta 
potencialmente valiosa para estudiar los procesos de cambio 
global de los sistemas dinámicos formados por una cantidad 
muy grande de componentes y, efectivamente, logró grandes 
avances durante el siglo pasado, en especial cuando, como he- 
mos visto en el capítulo anterior, aplicó los conceptos de la 
mecánica estadística para vincular la presión y la temperatura 
con los efectos promedio del movimiento de un número muy 
grande de moléculas. 

Se llegó así a interpretar satisfactoriamente mediante la ley 
de entropía la evolución hacia el equilibrio termodinámico de 
las estructuras aisladas que aparecen en la física y la química, 
cuyo ejemplo más conocido es el de los gases ideales. 

La segunda ley de la termodinámica establece que, en general, 
las estructuras aisladas terminan por descomponerse, quedan- 
do reducidas a un movimiento desorganizado de sus elemen- 
tos: las nubes de humo se disipan, las zonas calientes y las 
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zonas frías en un objeto se desvanecen hasta llegar a una tem- 
peratura uniforme. 

La cantidad de entropía en el Universo —lo aleatorio, o des- 
ordenado— sólo puede aumentar, hasta llegar al máximo, dice 
la segunda ley. Pero, curiosamente, en la misma época en que 
se enunciaba esta ley aparecía la teoría de la selección natural 
para explicar la evolución de las especies. Mediante dicha 
teoría, Carlos Darwin intentaba dar cuentá del hecho de que 
los organismos vivos constituyen estructuras cada vez más orga- 
nizadas: a partir de las bacterias se llega hoy a los mamíferos y 
al hombre. ¿Contradice este proceso de organización crecien- 
te la segunda ley de la termodinámica? En realidad, la exis- 
tencia de seres vivos no es un desafío a esta ley, que se aplica a 
sistemas aislados. 


Los SISTEMAS ABIERTOS PUEDEN ORGANIZARSE 


Un sistema vivo es abierto: un hombre absorbe energía y mate- 
ria de fuentes externas (el calor del Sol, el aire, carnes, verdu- 
ras, fuentes que a su vez están estructuradas y, por lo tanto, son 
de baja entropía) y expele sus productos de desecho, que son de 
alta entropía por ser el resultado de la descomposición de ma- 
teria organizada, a otros sistemas abiertos de su medio am- 
biente. Mientras un organismo esté vivo, se mantiene lejos del 
equilibrio termodinámico al que tienden los sistemas aislados. 

Sabemos muy bien que si alguien se aísla por completo del 
ambiente exterior, en muy poco tiempo la inexorable segunda 
ley se cumple llevándolo a la muerte. 

Durante bastante tiempo, muchos científicos pensaron que 
las leyes fundamentales de la física sólo permitían deducir 
que los sistemas deben llegar al equilibrio termodinámico y que, 
por lo tanto, el proceso de evolución biológica, con una cre- 
ciente complejidad de los organismos vivientes, era una rara 
excepción. 

Pero ni Boltzmann ni Darwin habían podido considerar la 
existencia de los fenómenos de formación espontánea de es- 
tructuras en la materia, una propiedad que se ha comenzado a 
estudiar en estas últimas décadas y que aparece con sorpren- 
dente frecuencia en la naturaleza. 

Así, sistemas tan simples como una capa de líquido o una 
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mezcla de productos químicos pueden, en ciertas condiciones, 
exhibir fenómenos de comportamiento coherente espontáneo. 
Una condición esencial para que esto ocurra es que se trate de 
sistemas abiertos, mantenidos lejos del equilibrio termodiná- 
mico por fuentes de energía. 

Es realmente extraordinario que enormes conjuntos de par- 
tículas, sujetas sólo a las ciegas fuerzas de la naturaleza, sean, 
sin embargo, capaces de organizarse en configuraciones de 
actividad coherente. Uno de los grupos que han promovido el 
estudio de este tipo de procesos —desarrollando así la rama de 
la termodinámica de sistemas lejos del equilibrio—, ha sido el 
de Ilya Prigogine y sus colaboradores, razón por la cual Pri- 
gogine recibió el premio Nobel de química en 1977. 

Estos procesos ayudan además a comprender los mecanis- 
mos que llevan a las oscilaciones en ciertas reacciones quími- 
cas, conocimiento que tiene una gran trascendencia para los 
procesos catalíticos industriales, sin mencionar la importancia 
de este tipo de reacciones para los procesos bioquímicos en los 
seres vivos. 


REACCIONES QUÍMICAS MUY EXTRAÑAS 


Hemos visto que los fenómenos de convección térmica pue- 
den llevar a un líquido o gas que parte inicialmente de un es- 
tado homogéneo a evolucionar y estructurarse con el tiempo, 
dando origen a formas regulares. 

Una condición indispensable para que aparezcan células 
como las de Benard es que el sistema sea abierto, es decir con 
una fuente externa de energía, y que esté alejado del equilibrio. 

Esta formación de células por convección es muy llamativa, 
pero parece un humilde fenómeno si se compara con los 
espectaculares efectos de las oscilaciones químicas, que se han 
comenzado a estudiar metódicamente a partir de 1980. 

En 1958, el bioquímico ruso Boris Belousov hizo una mezcla 
de ciertos productos químicos que usualmente forman un lí 
quido incoloro hasta que reaccionan, tornándose entonces de 
color amarillo pálido: Belousov había mezclado”los ingredien- 
tes sin preocuparse por las proporciones utilizadas, y observó 
sorprendido que la solución cambiaba periódicamente su co- 
lor, pasando a intervalos regulares de incolora a amarillo páli- 
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do para volver a hacerse incolora, lo que significaba que la 
reacción retrocedía y volvía a avanzar como si no pudiese de- 
cidir qué sentido tomar. 

El pobre Belousov intentó publicar su descubrimiento en las 
revistas científicas, pero éstas se rehusaron. Los árbitros que 
evaluaron el trabajo consideraron que la única explicación 
posible para ese fenómenó era una mezcla ineficiente de los 
reactivos, pues las leyes de la termodinámica se oponen a la 
existencia de tales oscilaciones. Belousov murió sin haber lo- 
grado que se reconocieran sus investigaciones. 

El escepticismo de los químicos de la época no debe sor- 
prendernos demasiado. En efecto, en una reacción química 
típica entre dos reactivos A y B, sus moléculas se están movien- 
do al azar, y si una A y otra B colisionan pueden combinarse 
para formar una molécula C y otra D, que son los llamados 
productos de reacción. Esto se simboliza como: 


A+B>C+D 


Los reactivos A y B van desapareciendo progresivamente a 
medida que aumenta la proporción de los productos C y D. 
Sin embargo, en un sistema aislado se observa que los reactivos 
A y B nunca se agotan del todo y que, con el tiempo, los cua- 
tro, A, B, € y D, coexisten manteniendo fija la proporción de 
cada uno en la solución. Esta proporción ya no varía y se dice 


entonces que el sistema está en equilibrio químico, que corre-' 


sponde al equilibrio al que llegan los sistemas según la segun- 
da ley de la termodinámica. 

Pero también se puede transformar este sistema en abierto, 
por ejemplo agregando continuamente más reactivos A y B, o 
retirando parte del Co D que se producen. Aquí también se 
puede comprobar que mediante la combinación adecuada de 
los flujos de entrada y salida, las proporciones van variando pro- 
gresivamente con el tiempo hasta hacerse fijas, aunque con va- 
lores distintos a los del sistema aislado, resultado que, una vez 
más, es esperable de acuerdo con las leyes de la termodinámica. 

Pero lo que afirmaba Belousov implicaba que estas combina- 
ciones en lugar de crecer progresivamente con el tiempo hasta 
alcanzar un equilibrio estable, podían retroceder hacia el esta- 
do inicial, lo que equivale a contradecir el segundo principio 
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de la termodinámica, y además podían hacerlo repetidamente, 
oscilando en un sentido y el otro. 

Los fenómenos periódicos o de oscilación son ubicuos en la 
física, la astronomía y la biología, pero los químicos creían que 
las reacciones eran inmunes a esa clase de comportamiento. 

No se espera usualmente que las concentraciones de los pro- 
ductos intermedios de la reacción alcancen un cierto nivel, 


luego caigan a uno inferior, luego suban y caigan repetidamen- 


te hasta que en algún punto resulten productos estables resis- 
tentes a cambios ulteriores. De todos modos, en algunos labo- 
ratorios de Moscú se siguió utilizando la “receta” de Belousov 
como una curiosidad de la química, y al comenzar la década 
de los sesenta, Anatol Zhabotinskii retomó el tema para su 
tesis de doctorado. Hizo investigaciones sistemáticas cambian- 
do los reactivos para obtener colores que contrastaran mejor. 

En la reacción ahora llamada BZ (Belousov-Zhabotinskii), la 
oscilación se manifiesta por el cambio regular entre el rojo y el 
azul. Se logra disolviendo en agua a una determinada tempe- 
ratura ciertas proporciones de ácido sulfúrico, ácido malónico, 
bromato de potasio y sales de cerio y hierro. Las oscilaciones en- 
tre el rojo y el azul resultantes tienen un periodo de casi un 
minuto y pueden durar varias horas. 

Además pueden producirse oscilaciones en el espacio: si se 
vierte algo de esa solución en un disco formando una capa de 
poco espesor, aparecen unas figuras hermosas, con círculos con- 
céntricos o espirales azules sobre un fondo rojo, que giran en 
un sentido o el opuesto y que van cambiando con el tiempo. La 
figura VIL] muestra una sucesión de fotografías de estas ondas. 

Actualmente la reacción BZ y otras similares se estudian en 
lo que se denomina reactor de flujo continuo, introduciendo 
continuamente los reactivos y retirando el exceso de los pro- 
ductos para mantener el volumen constante. Se puede con- 
trolar así rigurosamente la reacción y mantenerla como un sis- 
tema químico abierto, lejos del equilibrio. 


REALIMENTACIÓN Y CATÁLISIS 


Al día de hoy se han estudiado en detalle todas las etapas de 
este tipo de reacciones, planteando las correspondientes ecua- 
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ciones, y con ellas se ha logrado simular por computadora las 
oscilaciones. Se sabe ahora que las condiciones para que éstas 
aparezcan son que, además de tratarse de un sistema alejado 
del equilibrio, debe haber una realimentación, es decir que 
alguno de los productos que aparecen en una etapa del proce- 
so sea capaz de influir sobre su propia velocidad de formación. 
En la química esto se denomina “autocatálisis”. 

Un catalizador influye sobre la velocidad con la que reaccio- 
nan químicamente las sustancias presentes y se mantiene sin 
cambios durante este proceso. En la autocatálisis, si la sustan- 
cia que se va produciendo actúa sobre su propia velocidad de 
producción incrementándola, lo que aumenta su concentra- 
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ción y, a su vez, produce más cantidad de esa sustancia, se tie- 
ne una realimentación positiva, que responde a una ecuación 
no lineal. 

Variando el flujo de entrada de reactivos en el reactor de flu- 
jo continuo, el sistema químico puede pasar a uno cualquiera 
de dos estados estables diferentes, cada uno de los cuales tiene 
su propio periodo de oscilación, es decir que: aparecen bifurca- 
ciones para la posible evolución del sistema, las que recuerdan 
a las encontradas en la dinámica de flujo de fluidos, 

Para algunos rangos del flujo en el reactor aparece un com- 
portamiento más complejo, con una mezcla de varias frecuen- 
cias de oscilación de diversas amplitudes. La explicación co- 
rriente es que, a diferencia de una reacción química “normal”, 
en la que los reactivos y productos se mantienen uniforme- 
mente distribuidos en la solución, aquí una pequeñísima he- 
terogeneidad o diferencia en un punto se incrementa por el 
efecto de autocatálisis, de modo que, en una región, una sus- 
tancia química determinada puede dominar la reacción mien- 
tras que, en la vecina, su concentración se agota. Así se activan 
complejas oscilaciones en el sistema de reacción. 

Un oscilador químico, por ejemplo, puede ser inicialmente 
de un color rojo uniforme y, cuando sigue la reacción, apare- 
cen puntos blancos que se transforman en anillos concéntricos 
azules, que se aniquilan al colisionar entre sí. Se tiene enton- 
ces la situación de sensibilidad a las condiciones iniciales que, 
como hemos visto, caracteriza a los sistemas dinámicos no li- 
neales y, efectivamente, representando la evolución de estos 
sistemas químicos en el espacio de las fases aparecen atrac- 
tores extraños. 

Hay que tener en cuenta que se emplean catalizadores en 
muchos procesos importantes, como los vinculados al petró- 
leo por ejemplo, de modo que la comprensión detallada de 
la dinámica de las reacciones químicas que son significativas 
para la industria puede tener grandes consecuencias eco- 
nómicas. 

Por otra parte, estos relojes químicos recuerdan a los ritmos 
naturales que aparecen en los organismos vivos. Si se tiene en 
cuenta que en las reacciones bioquímicas también aparecen 
las realimentaciones, producidas por el efecto catalítico de las 
enzimas, y que los seres vivos son sistemas alejados del equili- 
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brio, es de esperar que este nuevo enfoque ayude a compren- 
der el comportamiento de muchos mecanismos biológicos. 


VII ¿Cómo definiremos la complejidad? 


HasTA ahora hemos usado el término complejidad para un 
estado en el que muchos factores diferentes interactúan entre 
sí. Pero hay que dar mayor precisión a este concepto, ya que la 
complejidad de un sistema no debe confundirse con que un 
sistema sea meramente complicado. En realidad debería ha- 
blarse de comportamiento complejo en un sistema, ya que, 
como hemos visto, un sistema dinámico puede ser muy simple 
pero exhibir en ciertas condiciones un comportamiento ines- 
perado de características muy complejas que llamamos caótico. 


ENTRE EL ORDEN Y EL CAOS: LA COMPLEJIDAD 


Un cristal de cuarzo es un sistema ordenado, con sus átomos 
vibrando alrededor de posiciones netamente fijadas en la red 
cristalina; un virus tiene características de orden en su estruc- 
tura, similar a un cristal orgánico, pero cuando infecta una 
célula, súbitamente comienza a replicar genes como un orga- 
nismo vivo: es complejo; el movimiento de las moléculas de un 
gas en equilibrio térmico es verdaderamente caótico. La com- 
plejidad cubre así un vasto territorio que está entre el orden y 
el caos. 

No hay actualmente acuerdo acerca del significado de “com- 
plejidad”. Todos sabemos que los seres humanos tienen mucha 
complejidad, y que los peces tienen bastante; los hongos son 
algo complejos, y los virus, menos. Las galaxias son complejas, 
pero en un sentido diferente. 

En todo esto coincidimos, pero ¿podríamos ponernos de 
acuerdo en usar una definición más científica, que especificara 
cuantitativamente cuánta más complejidad tiene una planta 
que una bacteria? 
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En las ciencias físicas, una ley no es ley si no se le puede 
expresar en lenguaje matemático; ningún científico se entu- 
siasmaría por medidas de complejidad que sonaran así: “un 
pez es bastante más complejo que un virus, pero un poco me- 
nos que un mamífero”. 

Podría parecer que esta dificultad en hallar una buena de- 
finición de complejidad desaparece si nos 'olvidamos de los 
organismos vivos y la aplicamos sólo a los sistemas dinámicos 
formados por materia inerte, en los que es más sencillo saber 
si son o no complicados. 

Pero tampoco es exactamente así. Según las condiciones del 
sistema, se puede tener orden o complejidad y hasta caos, ya se 
trate de uno con sólo tres componentes o con una enorme 
cantidad de ellos. 

Es el territorio de la complejidad, que está situado entre el or- 
den y el caos, el que constituye el nuevo desafío para la ciencia. 

De todos modos, resulta difícil poder predecir que esta 
búsqueda de una definición precisa de la complejidad termine 
con el hallazgo de una única magnitud que nos dé un número 
para ingresar a las ecuaciones de la física, como ocurre con la 
velocidad, la presión o la masa. 

Al fin y al cabo, la complejidad de una galaxia implica una 
relación entre sus estrellas componentes cualitativamente 
diferente de la que existe entre las células que forman a un 
cordero. No parece que esta diferencia se pueda reducir a una 
única magnitud similar a la temperatura. 

Posiblemente la mejor solución al problema sea la propuesta 
por los especialistas en computación. En la actualidad la ma- 
yoría de los científicos que se dedican a tratar de definir la 
complejidad pertenecen a dicha área, cuya tradición es ver 
prácticamente todo como reducible a información y, por lo 
tanto, cuantificable en términos de bits y bytes. 


¿CÓMO MEDIR LA COMPLEJIDAD? 


Si se adopta este punto de vista, entonces se podría medir la 
complejidad en función de, por ejemplo, el tiempo que re- 
quiere una computadora para ejecutar un programa que si- 
mule un proceso físico complejo. 
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En principio no hay límite para los tipos de procesos que se 
pueden simular con un programa, de modo que ésta sería una 
buena base de partida para la definición que buscamos. Así, 
los fenómenos complejos incluyen no solamente enjambres de 
símbolos (números, programas de Se palabras), 
sino procesos físicos y organismos vivos. 

Dicho enfoque comenzó con la publicación en 1948 de un 
trabajo de Claude Shannon, de Bell Telephones, sobre la teo- 
ría matemática de la comunicación, que condujo a buscar un 
equivalente de la segunda ley de la termodinámica para la 
información. 

Como consecuencia de esa presentación, los científicos se 
habituaron a mirar los entes físicos —bicicletas o flujos de pe- 
tróleo— y preguntarse: ¿cuánta información requiere descri- 
bir este sistema? 

Por otra parte el fenomenal avance en la capacidad de las com- 
putadoras ha impulsado el enfoque computacional de los pro- 
cesos en sistemas físicos, combinando observaciones de dichos 
sistemas con la construcción de modelos para computadora. 
Se simula así un proceso, para obtener resultados que tra- 
dicionalmente sólo se lograban modificando las condiciones 
físicas del sistema y midiendo su efecto. Esto resulta especial- 
mente útil para sistemas como los que se estudian en cosmo- 
logía, donde no se puede alterar por ejemplo la estructura de 
un planeta, o en las ciencias sociales, donde no se pueden 
cambiar las condiciones económicas de una nación para ver 
qué variables son las importantes. 

Puesto que muchos procesos de la naturaleza son expresa- 
bles mediante modelos que permiten simular su evolución con 
una computadora, se podría intentar medir la complejidad de 
un sistema por la dificultad para representarlo con un algo- 
ritmo, es decir con un programa de computadora. ¿Cómo se 
mide esa dificultad? 

Se han propuesto diferentes fórmulas, como la de medir el 
tiempo mínimo necesario para que una máquina ejecute el pro- 
grama, o también medir la capacidad mínima de memoria que 
debe tener la computadora para correr ese programa. Pero 
como estas magnitudes dependen del tipo de máquina utiliza- 
da, es necesario referirlas a alguna computadora ideal que 
actúe como un patrón normalizador. 
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LA MÁQUINA DE TURING 


Esta máquina abstracta fue ideada por el matemático inglés 
Alan Turing en 1936. Se le puede considerar como un artefac- 
to mecánico con una impresora, a través del cual pasa la me- 
moria, consistente en una cinta de papel tan larga como sea 
necesario y sobre la que está marcada una secuencia de espa- 
cios; cada uno de los espacios puede estar en blanco o marca- 
do con una raya, lo que equivale a la notación binaria 0 y 1 
respectivamente (véase la figura VINL.1). La máquina puede 
realizar una de cuatro operaciones cada vez que utiliza la me- 
moria, o sea cuando va recorriendo los espacios de la cinta y 
los lee: moverse un espacio hacia adelante, o retroceder un es- 
pacio, o borrar una raya, o imprimir una raya. Como la cinta es 
tan larga como se necesite, se entiende que la máquina tiene 
capacidad ilimitada de almacenamiento de datos, además de 
que no se le fija limite de tiempo para completar sus ope- 
raciones. 


Memoria 


FIGURA VIN. 1. 


La máquina comienza a operar en cuanto se le introduce el 
programa, y continúa hasta que termina de imprimir el resul- 
tado. De este modo es capaz de efectuar todo programa expre- 
sado en código binario, o sea, en un lenguaje matemático for- 
mado por los dos únicos signos O y 1. 

Lo que caracteriza a la máquina universal de Turing es que, 
dado el programa de entrada adecuado, puede simular el 
comportamiento de cualquier otra computadora digital, aun 
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la más compleja. Obviamente es muchísimo más lenta, de 
modo que nadie ha intentado construir una, a pesar de su 


simplicidad. 

Tenemos ahora el instrumento para medir con precisión la 
cantidad de información que enviamos en un mensaje: la uni- 
dad fundamental de información es el bit, que se define como 
la información más pequeña capaz de indicar una elección 
entre dos cosas igualmente posibles. En la notación binaria, 
un bit equivale a un dígito, y puede ser 0 o 1. Así, por ejemplo, 
el número 29 se escribe 11101 en el código binario y, por lo 
tanto, contiene cinco bits de información. 

He aquí, entonces, la normalización buscada, que permiti- 
ría medir la complejidad independientemente del tipo de 
computadora empleado. 


HAY NÚMEROS QUE SON COMPUTABLES 
Y OTROS QUE NO LO SON 


Tratemos ahora de definir con más precisión la complejidad. 
Parece sensato intentar hacerlo empezando con objetos abstrac- 
tos como lo son los números, ya que las relaciones entre los 
mismos ya son, obviamente, cuantitativas. 

Turing distinguió dos clases de números: los computables y 
los no computables. Computables son aquéllos para los que 
existe un algoritmo, o programa de computadora, que corrido 
en la máquina nos proporciona el número, no importa cuán 
grande sea éste, que inclusive puede ser infinito. 

Para aclarar este concepto vamos a suponer que tenemos 
una computadora conectada vía satélite con otra que utiliza un 
amigo que vive en Japón, al que necesitamos transmitirle cier- 
tos números. Como sabemos que cuanto más segundos de 
transmisión requiera nuestro mensaje tanto mayor será la cuen- 
ta que nos presentará la compañía, nos conviene comprimirlo 
al máximo. Para ello conviene ver si tiene alguna característica 
que ayude a tal fin. : 

Supongamos que deseamos enviar un número como, por 
ejemplo: 


1234567891011121314151617181920212223242526272829303132 
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Al examinarlo hallamos que sus dígitos se forman escribien- 
do los primeros 32 números enteros en orden. Su estructura 
sigue una ley perfectamente determinada, por la que sé que el 
dígito 7 está seguido por el 8, luego viene el 9, etcétera. Por lo 
tanto, bastaría con transmitir un mensaje con las instrucciones 
para construir ese número. 

El correspondiente algoritmo a transmitir es un programa 
muy corto, que básicamente ejecute la siguiente instrucción: 


“Imprima los 32 primeros enteros en orden” 


La computadora de nuestro amigo le proporcionará el nú- 
mero, que puede ser tan grande como queramos cambiando 
muy poco la longitud del programa (por ejemplo, “Imprima 
los 1000 primeros enteros en orden”). 

En cambio, un número no computable es aquel cuyo único 
algoritmo posible resulta ser ese mismo número escrito en el 
programa. 

Esto lo entenderemos fácilmente si queremos transmitirle a 
nuestro amigo de Japón un número que hemos construido 
arrojando un dado muchas veces y anotando las cifras del 1 al 
6 que se van obteniendo sucesivamente. Así después de tirar el 
dado 25 veces obtuvimos este número: 


3546221356431142652142663 


En este caso, como cada dígito fue generado al azar, no tie- 
ne relación alguna con el que lo sigue o precede, de modo que 
el único algoritmo que puede reproducir el número es el pro- 
grama que simplemente lo copia: 


“Imprima 3546221356431142652142663” 


Del mismo modo, si deseara transmitir el número que obtu- 
ve cuando seguí arrojando el dado hasta completar mil tiradas, 
no tengo otra posibilidad que recurrir al programa: 


“Imprima 3546221356431142652142663...” 
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Donde los puntos suspensivos significan los otros 975 dígitos 
del 1 al 6. Deberé resignarme al hecho de que el mensaje con 
el programa sigue siendo tan largo como el número. 

Resumiendo, para un número computable es posible escri- 
bir un programa de computadora relativamente corto que lo 
calcule, aunque el número sea infinitamente largo. En cam- 
bio, para un número generado al azar, al no ser computable, el 
programa que lo calcule ha de contener ese mismo número y 
será, por lo menos, tan largo como aquél. 

Así, la complejidad de un número se podría medir por la 
cantidad mínima de instrucciones de un programa de máqui- 
na de Turing —es decir, su longitud-mínima— capaz de repro- 
ducir el número. : 

En consecuencia, un número generado al azar corresponde 
a una alta complejidad —el caos—, mientras que un número 
computable es ordenado, tiene baja complejidad. En el medio 
quedarían los números complicados, generados por una com- 
binación de aleatoriedad y orden. 

Pero, ¿qué relación puede haber entre estas consideraciones 
válidas para entes matemáticos y los sistemas de objetos de la 
naturaleza cuya complejidad intentamos estimar? 

Una pista para ver la relación está dada por el hecho de que 
el código binario que usa la computadora para expresar cual- 
quier número puede expresar también cualquier otra informa- 
ción. Las 28 letras del alfabeto se pueden hacer corresponder 
con 28 diferentes secuencias de unos y ceros, lo mismo que los 
signos de puntuación, etcétera; esto nos permite, por ejemplo, 
considerar la complejidad en la construcción de palabras en 
un idioma. 

Algo ordenado sería una hilera de letras como aaaaa, por- 
que se le puede escribir en un programa concisamente co- 
mo “5 x a”. En cambio, una secuencia al azar de letras, como 
defiksivgahglkjthlakijgueernsedgmk, es caótica, no tiene un pro- : 
grama que la imprima que sea más conciso que la secuencia 
misma. Entre ambos extremos está la complejidad que aparece 
en la forma en que se agrupan las letras para formar palabras 
y frases. 

Este método también permite codificar la información que 
corresponde a los objetos físicos. Así, un objeto ordenado, 
como un trozo de cristal formado por átomos de carbono (un 
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diamante), se podría formular en un algoritmo que especi- 
ficara la cantidad de átomos de carbono que lo forman como, 
por ejemplo, “10 x carbono”. 

Del mismo modo, las plantas, personas y bacterias serían si- 
milares, en este enfoque, a palabras, frases y párrafos: mezclas 
de orden y aleatoriedad. 


La complejidad algorítmica se define como la longitud del pro- 


grama más corto que puede hacer una computación. 


Un planteo semejante se encuentra en la teoría algorítmica 
de la complejidad formulada por Andrei Kolmogorov, Gregory 
Chaitin y Ray Solomonov. Esta teoría considera la compu- 
tación de la magnitud física Q especificada como una secuen- 
cia de dígitos S, y establece que $, y por lo tanto Q, es aleatorio 
si el programa de cómputo mínimo requerido para obtener $ 
es el programa que todo lo que hace es copiar S: 


“Imprima S” 


También se puede definir por este método la complejidad 
“computacional, que es la cantidad de tiempo que necesita 
una computadora para resolver determinado problema y que 
mide, por lo tanto, su dificultad. 

Se comprueba que hay dos tipos básicos de sistemas dinámi- 
cos cuyos movimientos se pueden computar: los sistemas orde- 
nados y los sistemas caóticos. Los primeros, como es el caso del 
oscilador armónico, tienen órbitas que requieren una canti- 
dad de información de entrada, o sea una secuencia de unos y 
ceros, cuya longitud es relativamente corta, pero el proceso de 
computación correspondiente da como resultado una infor- 
mación de salida mucho mayor. La información de entrada es 
de corta longitud debido a que las trayectorias inicialmente 
cercanas del sistema se apartan con lentitud a medida que 
pasa el tiempo, de modo que no es necesario conocer todos 
los dígitos de la información inicial para poder predecir me- 
diante la computación los estados futuros. 

Desde el punto de vista de la complejidad computacional, el 
tiempo requerido para computar una trayectoria es propor- 
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cional al logaritmo del tiempo del sistema. Así, si se calcula la 
trayectoria para el tiempo t= 15 segundos, se requerirá un 
tiempo de computación de 1.17, y para t= 150 000 segundos, 
un tiempo de 5.18, siempre mucho menor que el que requiere 
el sistema para recorrer esa trayectoria. 

En cambio, para los sistemas caóticos, la cantidad de infor- 
mación de entrada es la misma que la de salida, ya que las tra- 
yectorias inicialmente cercanas divergen exponencialmente 
con el tiempo, y esto trae como consecuencia que cada segun- 
do de tiempo del sistema se pierda un dígito decimal de pre- 
cisión en la información de salida. Para mantener el grado de 
exactitud en la información de salida es necesario, entonces, 
compensar esa pérdida agregando un dígito más de precisión 
en la información de la entrada y, por lo tanto, si se desea inte- 
grar con exactitud una trayectoria caótica, hay que introducir 
tanta información como la que se extrae. 

Si se examina esto desde el punto de vista de la complejidad 
computacional en lugar de la complejidad algorítmica, se ve 
que el tiempo que requiere la computadora para calcular una 
trayectoria caótica es proporcional al que requiere ese sistema 
caótico para recorrerla, de tal forma que la computadora ne- 
cesita tanto tiempo para calcular como el sistema caótico mis- 
mo para ejecutar el proceso, y por lo tanto, no es posible hacer 
una predicción. 

Todos estos métodos permiten cuantificar la complejidad de 
los procesos extremos de orden y caos, pero no son muy satis- 
factorios para la complejidad que está entre ambos. En efecto, 
la complejidad de un sistema dinámico como el de un organis- 
mo vivo es, siguiendo las definiciones anteriores, menor que la 
de un gas en equilibrio termodinámico, ya que en este último 
la entropía ha llegado a su máximo y el movimiento de sus áto- 
mos es totalmente aleatorio. Esto parecería implicar que com- 
plejidad y aleatoriedad son equivalentes, lo que no lleva a nin- 
guna parte si la intención es comprender la complejidad de 
los sistemas que exhiben una organización, como se encuen- 
tran abundantemente en la naturaleza. La idea de comple- 
Jidad no coincide con la de entropía, pues hay sistemas que 
pueden tener igual orden, o igual desorden, y diferir en com- 
plejidad. 

C. Bennett, de 18m Research, y otros investigadores han pro- 
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puesto un enfoque diferente, que dé como resultado una de- 
finición con la que la complejidad de un pez sea mayor que la 
de un cristal o que la de un gas. 

La evolución de la vida en la Tierra nos muestra organismos 
cada vez más organizados, es decir que algo que podemos lla- 
mar complejidad aparece como aumentando, a pesar de que 
la cantidad de orden en el Universo está disminuyendo. 

La propuesta de Bennett es medir ese proceso de organiza- 
ción, en especial para los sistemas autoorganizados que apare- 
cen en la naturaleza. Lo que.los caracteriza es que, además de 
ser complejos, parten inicialmente de sistemas simples (una lá- 
mina de fluido que se calienta, una solución de sustancias quí- 
micas, un huevo fecundado...). La idea de organización o com- 
plejidad de un sistema estaría así estrechamente vinculada al 
proceso por el que se puede ir de ese sistema inicial simple al sis- 
tema complejo totalmente desarrollado. 

La complejidad se mediría entonces por el tiempo que le 
llevaría a una computadora simular todo el desarrollo de ese 
sistema hasta llegar a su organización definitiva, y contando el 
número de “pasos lógicos” a lo largo del proceso. Se tienen 
como algoritmo de entrada las reglas básicas de su crecimien- 
to, que incluyan, por ejemplo, los datos genéticos, si se trata de 
los seres humanos. Así, considerando como teoría del origen 
de la especie humana la selección natural, el número de pasos 
lógicos sería el número estimado de veces que nuestro mate- 
rial genético ha sido modificado a partir del que estaba con- 
tenido en una bacteria inicial. La complejidad se mediría en- 
tonces por el tiempo que requeriría una computadora para 
simular dicha evolución pasando por todos esos pasos lógicos. 

Pero ninguna de las mencionadas definiciones de la comple- 
jidad que está entre el orden y el caos ha sido aceptada por la 
comunidad científica como totalmente satisfactoria, hay nue- 
vas propuestas y el debate continúa aun. Por otra parte, la falta 
de definición no es impedimento para que el tema de la com- 
plejidad en sí esté registrando progresos importantes, en un 
proceso que se asemeja al que acompañó a la comprensión de 
la naturaleza del calor y su relación con la energía. Los avances 
en la ciencia del calor se fueron realizando desde principios 
del siglo xvII1, pero a mediados del siglo xIx se le pudo definir 
con precisión basándose en la teoría cinética. Es de esperar 
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que los avances actuales permitan definir la complejidad en un 


tiempo mucho más breve. 


IX. Aplicaciones en biología y economía 


EXISTE una ecuación matemática especialmente adecuada para 
examinar estas propiedades comunes a varios sistemas y dada 
su importancia pasaremos a describirla. 


LA ECUACIÓN LOGÍSTICA 


Se trata de la ecuación logística, que es una ecuación simple, 
muy fructífera por la cantidad de aplicaciones en muchos cam- 
pos de estudio de sistemas complejos: ecológicos, biológicos, 
económicos, etcétera. 

Es una ecuación que opera sobre un número y lo transfor- 
ma en otro sobre el que vuelve a operar, y así reiteradamente, 
en un proceso iterativo. Sus características especiales sólo se 
evidencian cuando el número de iteraciones es grande, de 
modo que para aplicarla se requiere al menos una calculadora 
de bolsillo, aunque una computadora es la herramienta más 
adecuada. 

La ecuación logística produce sobre un número cualquiera 
dos efectos que se oponen: 1) lo incrementa, produciendo otro 
número mayor que, a su vez, es vuelto a incrementar por la 
ecuación y así reiteradamente; 2) va reduciendo esos números 
resultantes a medida que crecen, de modo que se tiene aquí 
un proceso con una realimentación controlada. 

¿Qué ocurrirá cuando la ecuación haya operado un buen 
número de veces? El sentido común nos indicaría que debe 
resultar finalmente algún número intermedio, ni demasiado 
grande ni demasiado pequeño. 

Pero aquí viene la gran sorpresa: esto puede ser totalmente 
equivocado, tan equivocado que nos podemos topar también 
aquí con comportamientos insospechados, de aspecto caótico. 
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APLICACIONES EN BIOLOGÍA 


La ecuación lógica fue propuesta en 1845 por el sociólogo y 
matemático Pierre Verhulst; sus sorprendentes propiedades 
fueron puestas de manifiesto por el físico y biólogo Robert 
May en la década de 1970, cuando la aplicó al estudio de la di- 
námica de las poblaciones de plantas o animales. En dichas 
poblaciones hay una realimentación en cada ciclo vital por la 
reproducción controlada por el efecto negativo de los depre- 
dadores o de la creciente escasez de alimento, que impide así 
que esas poblaciones crezcan explosivamente. La ecuación per- 
mite calcular, partiendo de las características de la población 
en un momento dado, cómo irán variando éstas con el tiempo. 

Deseamos saber cómo variará anualmente la cantidad de in- 
dividuos de una población, de la que sabemos que en el año 
inicial hay 1000 y que aumenta con un ritmo constante de 
10% por año. 

Si no existiera el control de los predadores o de la disponi- 
bilidad de alimentos, podríamos hacer la siguiente tabla: 


Año Número de individuos Total 
0) 1000 1000 
1 1000 + 100 1100 
2 1100+110 1210 
3 1210 + 121 1331 


Y así sucesivamente. 
Este procedimiento se puede expresar matemáticamente con 
la ecuación: 


X,,¡=KxX, 


+1 

Donde tindica la cantidad de años transcurridos a partir del 
año inicial que es t= 0, y Xes la variable que simboliza la canti- 
dad de individuos. Así, X, es esa cantidad en el año 1, y X,,, en 
el año siguiente. : 

El parámetro K nos indica la tasa anual de aumento de X; en 
el ejemplo anterior es K= 1.1. 

Por lo tanto, lo que nos dice esta expresión es que si sabe- 
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mos cuál es la población en el año t, basta con multiplicar el 
número correspondiente por la tasa de aumento K para deter- 
minar la población que habrá en el año siguiente. 

Para facilitar los cálculos se trabaja con un X normalizado, 
es decir que sólo pueda variar entre 0 y 1. Así, X, = 1 corres- : 
ponde al máximo posible para esa población, o sea al 100%, y 
X,= 0.5 al 50%, y no importa si estamos considérando 15950 
conejos o 12 millones de árboles; todo lo que nos interesa es 
calcular la variación anual de la población en relación con los 
valores anteriores o siguientes. 

Volviendo al ejemplo, evidentemente ese incremento de 10% 
anual nos llevará a una situación imposible: si comenzamos 
con 1000 conejos en el año cero, 200 años después habría 
unos 19 mil millones de conejos cubriendo la superficie del 
planeta; se requiere entonces agregar a la ecuación un térmi- 
no que refleje la situación real, en la que el alimento no va a 
alcanzar para alimentar a la creciente población, y además 
habrá cada vez más zorros y otros predadores que se alimentan 
de conejos. 

La ecuación logística es, en definitiva, la siguiente: 


AER 


| Ejecuta así dos acciones opuestas; cuanto más crece el factor 
| X,, tanto más disminuye el factor (1 — X,) reduciendo el resul- 
| tado final. 
Si, por ejemplo, Kes 1.1, para un X, pequeño de 0.1, resulta- 
ba con la primera ecuación: 


X,,1=1.1x0.1=0.11 
Pero como el factor reductor es aquí (1 — 0.1) = 0.9 la nueva 


ecuación es: 


X,,1=0.11x 0.9 = 0.099 
O sea que el factor reductor afectó muy poco el resultado. 
Cuando X, crece mucho, por ejemplo a 0.8 del máximo, este 
factor pasa a ser 0.2, y se tiene entonces en lugar de: 


X..1=1.1x0.8 = 0.88 


124 


El resultado corregido por el factor reductor de 0.2: 


Xx, =1.1x0.8x0.2=0.176 


t+1 

Éste disminuye la población a la quinta parte. 

Tenemos entonces una expresión matemática que permite 
calcular sin ambigúedad el valor de X,; en otras palabras esta- 
mos tratando con sistemas dinámicos donde el futuro depen- 
de de manera determinista del pasado, sin incertidumbres. 

Toda la información sobre el sistema se halla contenida en 
la ecuación logística, y aplicando ésta podemos conocer cómo 
variará en el tiempo. Así, si X, es el valor inicial de la variable 
cuya evolución en el tiempo queremos conocer, será X; al lle- 
gar al tiempo t= 1: 


XxX, =KxX,x (1 - Xo) (a) 
Y para X, tendremos: 
X, =KxX, x (1-X,) (b) 
Donde puedo sustituir X, por la ecuación (a), resultando así: 
X,= K(Kx Xo(1 — Xp) (1-Kx Xo(1 - X0))) 
Y similarmente para t= 3: 


X, = K(K(KXo(1 - Xp) (1 — KXo(1 — Xp))) x 
(1 - K(KX, (1 - Xo) (1 - KX(1 — X0))))) 


Como se ve, la ecuación se convierte progresiva y rápidamen- 
te en una fórmula cada vez más complicada; si pretendiésemos 
conocer el futuro para t= 20 (20 años, o 20 generaciones, o lo 
que sea que usemos como medida del tiempo), necesitaríamos 
unas 300 páginas nada más que para plantear Xop, y para llegar 
a un tiempo ¿= 50, que no es algo exageradamente lejos en el 
futuro, no alcanzaría con el tamaño de una biblioteca. 

En la era de las computadoras, no tiene sentido seguir este 
camino; ya que si hay algo que esos artefactos saben hacer, es 
reiterar con gran velocidad la misma operación todas las ve- 


125 


ces que se desee, y además lo hacen sin equivocarse ni abu- 
rrirse. 

De modo que, en lugar de seguir el método clásico de es- 
cribir una ecuación que sea válida para todos los 1, XA, y Kpo- 
sibles, para luego introducir en ella los números correspon- 
dientes al caso que interese, lo que haremos es darle a la 
computadora los datos iniciales X, y K, la cantidad de itera- 
ciones que necesitamos, y las instrucciones para que ejecute el 
cálculo básico que como habíamos visto es muy elemental, 
consistiendo en hacer una sustracción y dos multiplicaciones. 

Con el programa así preparado cualquier computadora per- 
sonal nos permite obtener valores tales como Xy 009 O cual- 
quier otro en muy poco tiempo, y además hacer experiencias 
tales como ver paso a paso cómo va variando X, al incrementar 
t, y qué pasa cuando cambia K. 

Quien desee experimentar el poder del cálculo iterativo 
puede hacerlo sin necesidad de utilizar una computadora: bas- 
ta con una simple calculadora de bolsillo y una buena dosis de 
paciencia. 

Si dispone de ambas, propongo hacer un cálculo iterativo 
muy simple, que consiste en elevar un número al cuadrado, el 
resultado volver a elevarlo al cuadrado, y así sucesivamente 
hasta diez veces en total. 

Introduzcamos el número “0.9999” en la calculadora, y pul- 
semos la tecla “x?”, o si no se dispone de ésta, la de multiplicar, 
x”, seguida de “=”. Repitamos la misma operación diez veces. 
Una calculadora de ocho dígitos da como resultado: 


X, = 0.9999 
X;o = 0.9026637 


Veamos ahora qué ocurre cuando se varía ligeramente la: 
condición inicial. ¿Qué pasará si utilizamos como X, un valor 
que difiere en sólo 0.1% del anterior? Para: 


Xy = 0.9999 
Es ahora: 


Xio = 0.3589714 
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¡De modo que el resultado final es 40% del anterior, habien- 


do cambiado el valor inicial en sólo 0.1%! 

Tenemos aquí un ejemplo de la sensibilidad a las condicio- 
nes iniciales, que, como habíamos visto, es una condición ne- 
cesaria para que aparezcan los fenómenos caóticos y de com- 
plejidad en los sistemas dinámicos. 


HISTORIAS DE PECES Y CRUSTÁCEOS 


En estas operaciones matemáticas iterativas aparecen otras ca- 
racterísticas inesperadas, que se ponen de manifiesto a partir 
de algunas decenas de ciclos de cálculo, y que pasaremos a des- 
cribir. Si alguien desea hacerlo por sí mismo, es preferible que 
utilice una computadora, a menos que le divierta hacer largas 
cuentas con una calculadora. 

El programa con instrucciones para la computadora tiene 
esta forma general en lenguaje Basic: 


INPUT K 

X=0.6 

FOR n= 1 to 100 
X=Kx*Xx*(1-X) 
PRINT X 

NEXT n 

STOP 


Donde Kes la tasa de incremento, y hemos fijado el número 
de iteraciones n en 100 y el Xinicial en 0.6. 

Para demostrar el asombroso comportamiento de la ecuación 
logística, comenzaremos comparándola con los resultados de 
los estudios matemáticos que hizo Vito Volterra en la década 
de 1920 para explicar las fluctuaciones periódicas de la po- 
blación de peces del Mediterráneo. 

Consideremos una población de crustáceos, y de peces que 
se alimentan de aquéllos, y vamos a suponer que los crustáceos 
tienen una tasa de reproducción pequeña, K = 1.01, y que su 
población inicial es X, = 0.6. 

Haciendo el cálculo resulta que la población disminuye cada 
año de modo que, al cabo de un tiempo, la colonia desaparecerá: 
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t Xx t Xx 


0 0.6000 5 0.1147 
1 0.2424 10 0.0732 
€x _——_—_ __===————+ + 


¿Pero qué pasa cuando Kes mayor, por ejemplo K = 2? 


t X 1 X 

0 0.6000 3 0.5000 
1 0.4800 4 0.5000 
2 0.4992 5 0.5000 


La población queda estabilizada en 0.5. 

Si la tasa de reproducción aumenta a 2.7, la ecuación mues- 
tra una fluctuación anual, variando entre 0.61 y 0.64 por la 
oposición entre el crecimiento y la acción de los peces, pero 
finalmente, al cabo de unos 15 ciclos, se estabiliza en aproxi- 
madamente X= 0.63, valor que hace, por lo tanto, de atractor 
para este comportamiento. 

Pero es a partir de una tasa de reproducción K mayor que 
3.0 que sucede algo nuevo: el sistema fluctúa fuertemente al 
principio, y finalmente la población de la colonia oscila entre 
dos valores estables, de modo que el atractor se ha bifurcado 
en dos. 

Así, para K=3.3: 


t XxX t X 

0 0.60000 14 0.47941 
1 0.79200 15 0.82360 
2 0.54363 97 0.47943 
3 0.81872 98 0.82360 
4 0.48978 99 0.47943 
5 0.82466 100 0.82360 . 


x€ÉEO___ A 


Los valores estables son aquí 0.47943 y 0.82360, y cada uno 
de ellos se repite cada dos años. 
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Al representar en un gráfico la variación anual de la pobla- 
ción de crustáceos (véase la figura IX.1), vemos que una vez 
que las fluctuaciones se han estabilizado, si en un año aumen- 
tan a 0.83360, esto resulta ser un verdadero festín para los 
peces, que, a su vez, se incrementan tanto que bajan la canti- 
dad de crustáceos al nivel inferior de 0.47943 en la temporada 
siguiente. Esto disminuirá la población de peces por la escasez 
de alimento, permitiendo que al año siguiente vuelvan a abun- 
dar los crustáceos, y así cíclicamente. 


Población 


0.5 


1 2 3 4 


Figura IX.1. 


Población Años 1, 3, 5... 


crustáceos 


Años 2, 4, 6... 
A 


Población peces 


Figura IX.2. 


El proceso descrito se 
puede representar también 
como fue planteado por 
Volterra, en un espacio de 
las fases, lo que permite vi- 
sualizar el comportamiento 
cíclico del sistema peces- 
crustáceos en la figura 1X.2, 
donde el eje vertical nos 
indica la población de crus- 
táceos y el horizontal la de 
peces. 

Si iniciamos la colonia 
en A con X, = 0.6 crustá- 
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ceos irán incrementándose a lo largo de la curva 1, que marca 
simultáneamente el aumento de los predadores hasta llegar al 
nivel B. A partir de allí la población de peces predomina, lle- 
vando al sistema por la curva 2 con disminución de crustáceos, 
que arrastra hacia un mínimo a la población de peces, y esto se 
repite cíclicamente con un periodo, en este caso de dos tem- 
poradas. 

La ecuación logística se ajusta muy bien a otros ciclos natu- 
rales como los de insectos y bacterias, que son los estudiados 
por R. May en la década de 1970. 

Continuemos explorando la ecuación mediante la computa- 
dora, que cumple en esta época de fines del siglo xx un papel 
fundamental como herramienta de la ciencia, análogo al del 
microscopio y del telescopio en siglos anteriores. 

Si aumentamos K por encima de 3.45, los dos valores que se 
repetían cada dos años vuelven a hacerse inestables, y cada 
uno se bifurca para producir una población que oscila alrede- 
dor de cuatro valores diferentes que se repiten cada cuatro 
años, de modo que el periodo se ha duplicado. 


Para K= 3.45: 
AAA A A A A A a 
t XxX t Xx 
== A A e o e 
0 0.6 76 0.42713 
72 0.42688 97 0.84470 
73 0.84405 98 0.45258 
74 0.45412 99 0.85474 
75 0.85524 100 0.42835 


q>2>--PIM>¿MMM»« 

Los niveles estables se hallan alrededor de 0.42820, 0.45290, 
0.84410 y 0.85320. 

En K=3.56 se produce una nueva inestabilidad, con bifurca- 
ciones que producen ocho valores fijos, con periodos del do- 
ble de los anteriores, y esto vuelve a ocurrir en K= 3.596, dando 
16 valores; luego siguen apareciendo mas y más bifurcaciones, 
hasta que finalmente en K= 3.56999 se llega a un estado caóti- 
co, con infinitos valores de X para la colonia que oscilan en 
forma imprevisible entre 1 y 0. 
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LAS RUTAS AL CAOS 


Vemos en este ejemplo que un sistema dinámico puede partir 
de un estado ordenado y bajo ciertas condiciones evolucionar 
hacia un estado caótico. ] 

En realidad existen varias rutas entre el comportamiento 

regular, ordenado, y el estado caótico, impredecible, por las 
que pueden transitar los sistemas dinámicos; además en algu- 
nas etapas de las rutas aparecen comportamientos complejos 
que pueden darles propiedades sorprendentes. 

Dichas rutas se pueden transitar en ambos sentidos: sistemas 
dinámicos ordenados pueden pasar a un comportamiento alta- 
mente complejo o llegar hasta uno caótico, como ocurre con 
los péndulos y otros osciladores, o pasar del caos a una com- 
plejidad organizada, como en los relojes químicos o las células 
de Benard. 

Hemos visto que para que aparezcan estos comportamientos 
en los sistemas con una cantidad de componentes que puede 
ser pequeña o muy grande, es requisito indispensable que el 
número de variables que intervienen en la dinámica del sis- 
tema sea limitado. Si el número de variables independientes, o 
grados de libertad, es menor de tres, los comportamientos ana- 
lizados no aparecen; por otra parte, si dicho número es muy 
grande, será imposible distinguir con estos métodos si su 
evolución se debe a los fenómenos que estamos estudiando o 
si obedece a factores fundamentalmente aleatorios. 

Las rutas entre orden y caos se pueden clasificar en tres 
tipos principales de acuerdo con las diferentes modalidades 
con las que se producen dichas transiciones: 

1) Casi periodicidad. El sistema se representa en el espacio de 
las fases con un atractor casi periódico inscrito en un toro y la 
transición lo transforma en un atractor extraño. 

2) Cascadas subarmónicas. El sistema presenta oscilaciones de 
un cierto periodo T, y a partir de la transición se produce una 
bifurcación, apareciendo otras de periodo doble 27 en la tran- 
sición siguiente se agregan oscilaciones de periodo 47, luego 
8T y así sucesivamente en cascada. Esto se observa entre otros 
en los fenómenos de convección térmica (capítulo VI) y en la 
reacción BZ (capítulo VII). 

3) Intermitencias. El sistema produce esporádicamente fluc- 
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tuaciones de gran amplitud. Estas transiciones suelen aparecer 
en los procesos hidrodinámicos y también en las oscilaciones 
de circuitos electrónicos, donde se manifiestan como un rui- 
do de baja frecuencia que aparece ocasionalmente. 

Con base en la clasificación precedente, en el caso de la ecua- 
ción logística, la variación del parámetro Kimplica una ruta al 
caos por duplicación de periodos, es decir, la cascada subar- 
mónica. : 4 

Como vemos, el papel de Kes definir la complejidad del com- 
portamiento, y resulta conveniente visualizar su efecto sobre la 
evolución del sistema mediante un gráfico como el de la figu- 
ra IX.3 de la página siguiente, donde representaremos en el 
eje vertical la población y en el horizontal el valor de K. 


DONDE APARECE EL CAOS 


La figura muestra el panorama que encontró R. May con su 
modelo basado en la ecuación logística, donde iba variando el 
K por variación del suministro de alimento. La figura repre- 
senta la población en el eje vertical y el parámetro K en el eje 
horizontal. 

Halló así las bifurcaciones sucesivas que indican el aumento 
de las oscilaciones en la población de insectos, y que luego, 


2 3 3334 3.7 3.84 4 
K > 


FIGURA IX.3. 


132 


después de varias bifurcaciones, se entraba en la zona caótica 
que indica la figura, donde la población en el modelo fluctúa 
erráticamente, tal como ocurre en la realidad con la de los 
insectos. ; 

Hemos representado en los recuadros de la figura IX.3 algu- 
nas de estas situaciones, que responden a un valor específico 
de K. Así para el de K=1.01, la población P que inicialmente 
era de 0.6, disminuye con el tiempo T, representado en el eje 
horizontal, hasta llegar a cero: se extingue. 

El recuadro de K = 2 nos indica que para este valor, después 
de algunos ciclos con fluctuaciones la población se hace es- 
table. 

Para K = 3.3, la población oscila entre los dos valores 0.47943 
y 0.82360, tomando uno de estos valores en un ciclo, y el otro 
en el ciclo siguiente. Estos dos valores están además represen- 
tados en el gráfico principal: si en K =3.3 trazamos una línea 
vertical, ésta cortará a las dos ramas de la curva bifurcada, en 
los valores mencionados de P. 

Del mismo modo, para K= 3.4 se tienen cuatro valores que 
se repiten en los ciclos sucesivos, valores que corresponden en 
la curva principal a las intersecciones de una línea vertical con 
las cuatro ramas. En K= 3.7, hay tantos valores distintos de B 
que no es posible hallar relación alguna entre un ciclo y el 
siguiente; en la curva principal, una línea vertical pasará por 
miles de intersecciones, marcadas como puntos en las franjas 
verticales. Se está en la zona de caos determinista. Finalmente, 
para K= 4, se tiene el caos con infinitos valores posibles para 
cada ciclo. 

Hoy en día, una de las preguntas que se plantean los ecólo- 
gos es la de averiguar si el comportamiento que predice este 
modelo se da en las poblaciones reales. 

Las opiniones entre los biólogos sobre estas teorías están 
divididas: desde los que piensan que la imprevisibilidad inhe- 
rente al caos determinista es un factor muy importante para 
explicar la evolución de las especies, hasta los que consideran 
que las poblaciones reales no pueden tener dinámicas caóticas 
porque se extinguirían, y que las investigaciones en laborato- 
rios y por simulación con computadora se alejan demasiado de 
los casos que aparecen en la naturaleza. 

Ésta no es una tarea sencilla, porque en los sistemas ecológi- 
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cos resulta muy difícil separar los múltiples factores ambien- 
tales, y tampoco se pueden variar parámetros como la tasa de 
reproducción para ver cómo varía su efecto. 

De todos modos se han hecho experiencias sobre pobla- 
ciones animales o vegetales aisladas en el laboratorio, donde sí 
se pueden alterar algunas condiciones, variando, por ejemplo, 
la temperatura ambiente para acelerar el metabolismo de mos- 
cas de la carne o de protozoarios. : 

Los estudios sobre estas poblaciones experimentales revela- 
ron efectivamente las transiciones que corresponden a las 
primeras bifurcaciones, pero sin una duplicación de periodos 
tan clara y nítida como la que aparece en los sistemas físicos. 

Lo que aparece inequívocamente es la transición al caos. 

El otro método poderoso que se está utilizando para el estu- 
dio de la dinámica de poblaciones es la simulación por com- 
putadora, en la que los investigadores crean un modelo que 
genera series de datos que representan el tamaño de una po- 
blación, a lo largo de muchas generaciones. Este es un mundo 
imaginario, donde el investigador especifica a voluntad los fac- 
tores que gobiernan al sistema, y luego los analiza usando los 
métodos que se aplican a los sistemas vivientes reales. 

Este método matemático permite construir un espacio de las 
fases con tantas dimensiones como variables independientes 
están en juego, y buscar atractores extraños, que, si aparecen, 
son una evidencia de que se está ante un sistema dinámico no 
lineal, determinista. 

Una de las dificultades que se presentan es que cada especie 
usualmente interactúa con muchas otras, y para cada una debe 
agregarse una variable o dimensión, de modo que hay que tra- 
bajar con espacios de las fases multidimensionales, donde es 
muy fácil confundir fluctuaciones estadísticas con la presencia 
de atractores. 

Se ha aplicado el método a casos como los de las epidemias 
de sarampión en Nueva York aparecidas durante un periodo de 
40 años, y éste reveló la presencia de un atractor tridimension- 
al en un espacio de cuatro dimensiones. 
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RITMOS EN LOS ORGANISMOS VIVOS 


Los métodos matemáticos se han aplicado también a la bús- 
queda de evidencias de caos determinista en la dinámica de 
organismos individuales, es decir, en los procesos fisiológicos y 
neurobiológicos rítmicos, lo que permite que se les pueda 
estudiar como conjuntos de osciladores que sé influencian mu- 
tuamente y que tienen ciclos de realimentación con una diná- 
mica no lineal. 

Un caso muy estudiado es el de los ritmos cardiacos, que po- 
dría arrojar luz sobre las arritmias y sobre la interpretación de 
los electrocardiogramas antes y después de un ataque car- 
diaco. 

El biofísico R. Cohen hizo una simulación por computadora 
de los ritmos cardiacos y comprobó que en el preludio de un 
infarto aparece una bifurcación de periodos de latido. Esto se 
explicaría considerando que los pulsos eléctricos que fuerzan 
a las fibras musculares de los ventrículos a contraerse y bom- 
bear así la sangre, obran como un sistema oscilador con un 
periodo regular. Si por alguna patología aparece una altera- 
ción del periodo de oscilación de los pulsos eléctricos de un 
grupo de fibras, estas dos oscilaciones diferentes combinadas 
pueden entrar en la situación que examinamos para dos pén- 
dulos acoplados, con una cascada de bifurcaciones de periodo 
2T, 4T..., hasta paralizar el corazón. 

Es claro que aquí se trata de un efecto que aparece en una 
simulación por computadora y que no es fácil obtener datos 
experimentales de un infarto, pues médicos y pacientes están 
bastante ocupados con tratar de superar el trance. De modo 
que lo más cercano al fenómeno real es hacer mediciones en 
un laboratorio. 

En 1980, el fisiólogo L. Glass inició una serie de experien- 
cias con células del corazón de embriones de pollo. Estas, en 
un medio de cultivo, siguen pulsando espontáneamente a un 
ritmo de 60 a 120 latidos por minuto, y son así un oscilador 
natural. Si se introduce un microelectrodo en esa masa de cé- 
lulas se les puede aplicar choques eléctricos periódicos, de 
modo que ahora se tiene un sistema oscilatorio con dos ritmos 
acoplados, uno intrínseco y otro forzado. Este último se puede 
variar, de modo que se observa cómo varía el latido cardiaco 
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apareciendo periodos 27, luego 47, y también situaciones to- 
talmente irregulares, caóticas que sugieren una fibrilación. 

Otra aplicación aparece en la neurofisiología, donde se ha 
analizado la complejidad dinámica de los electroencefalogra- 
mas (EEG) que registran la actividad cerebral de seres huma- 
nos realizando diversas tareas como, por ejemplo, contar des- 
de 700 hacia atrás de 7 en 7. Se computó la dimensión fractal 
del atractor que aparece al representar en el espacio de las fa- 
ses el conjunto de las oscilaciones registradas por los EEG, ha- 
llando que pasa del valor básico de 2.3 a un valor de alrededor 
de 2.9 cuando el sujeto está haciendo el esfuerzo de contar 
hacia atrás. La conclusión a que se llegó es que las formas de 
EEG de mayor dimensión fractal, o sea más complejas, corres- 
ponden a un estado mental más alerta. 

Pero todos estos estudios fisiológicos y neurológicos presen- 
tan el mismo flanco débil que el de las poblaciones animales: 
lo difícil que resulta aplicar estas teorías a los casos reales que 
aparecen en la naturaleza, y que son muchísimo más comple- 
jos que los modelos que se pueden plantear con una computa- 
dora. 

Las dificultades para la aplicación de estos métodos de estu- 
dio son reconocidas por R. May, quien, sin embargo, lo pro- 
mueve desde hace veinte años basándose en el hecho de que 
tanto las poblaciones como los procesos biológicos están go- 
bernados por mecanismos no lineales, y que, por lo tanto, de- 
bieran presentar comportamientos caóticos además de los cí- 
clicos estables. 

Parece todavía muy prematuro pronunciarse sobre el resul- 
tado final. Evidentemente es necesario seguir perfeccionando 
las técnicas de aplicación de un método que, al fin y al cabo, es 
muy reciente, antes de poder determinar hasta dónde estas 
teorías sobre las transiciones entre orden y caos constituyen la 
explicación de la dinámica de los seres vivientes. 


CICLOS ECONÓMICOS 
También en el campo de las ciencias sociales se está traba- 
. . Pp . .», 
jando intensamente en la aplicación de estos novedosos con- 


ceptos. 
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En particular resulta atractivo para los economistas intentar 
aumentar su capacidad de predicción de los mercados finan- 
cieros. Al fin y al cabo, los mercados financieros y las econo- 
mías de las naciones son sistemas dinámicos con mecanismos 
de realimentación y de autorregulación. Sabemos que si se ele- 
va el precio de un producto de manera excesiva, la demanda 
disminuye y el precio deberá caer. De modo que también aquí 
se puede intentar hacer un modelo basado en la ecuación 
logística: 

P,= AXP(1-P) 


Donde llamaremos P, a un precio en el día, mes o momento Í, 
Aa la tasa de aumento, y P,, , al precio en el tiempo siguiente. 

Como hemos visto anteriormente, las propiedades de esta 
ecuación matemática nos indican que para A menor que 3, 
después de cierto tiempo el precio P será estable, y que si A es 
mayor que 3, se tendrán fluctuaciones periódicas en el precio, 
con ciclos dobles que oscilan entre dos valores al principio, 
para pasar luego a cuatro, y así sucesivamente hasta que para A 
aún mayores, el valor de P puede tener comportamientos in- 
usuales, tales como no repetirse nunca, es decir ser caóticos. En 
tal situación, un gráfico de los valores de P puede ser fácilmen- 
te confundido con una serie de números generados al azar. 

En los últimos años varios grupos de investigación encara- 
ron este problema de aplicar el enfoque que utiliza la física de 
los sistemas dinámicos no lineales al campo de la economía. 

El más conocido es el que se nucleó en el Santa Fe Institute, 
en Nuevo México, donde en 1987, tres ganadores del premio 
Nobel (el economista Kenneth Arrow y los físicos Murray Gell- 
Mann y Philip Anderson) reunieron a economistas, físicos, 
biólogos y especialistas en computación para empezar un pro- 
grama de investigación sobre la economía considerada como 
un sistema dinámico complejo. 

Ellos han logrado construir modelos de computadora de la 
economía en los que los inversores son, a su vez, programas de 
computadora capaces de reconocer pautas de variación en los 
precios y de actuar en consecuencia, y que además aprenden 
con la experiencia. El resultado de esta simulación por compu- 
tadora comienza a acercarse a lo que se observa en las econo- 
mías reales. 
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Pero al igual que lo que vimos para las aplicaciones de una 
ecuación del mismo tipo a los sistemas ecológicos y a los orga- 
nismos vivos, la realidad económica es muchísimo más com- 
pleja, afectada por los constantes cambios en la sociedad, y el 
precio de cada producto está vinculado al de muchos otros. 

Así, es posible que los efectos de caos determinista que apa- 
recen en estos modelos matemáticos no estén presentes en 
una economía real, y que sea igualmente válido seguir descri- 
biendo la economía mediante procesos lineales con datos que 
tienen un importante componente de “ruido”, o sea de fluc- 
tuaciones al azar. 

Para decidir qué enfoque es el más adecuado, hay que uti- 
lizar los casos reales, que son muy difíciles de aplicar porque se 
requieren larguísimos conjuntos de datos, y aquéllos con los 
que se cuenta en general son escasos, además de tener un fuer- 
te componente de ruido. 

Pero en los últimos años se han inventado nuevas técnicas 
de análisis estadístico que son capaces de distinguir entre las 
fluctuaciones debidas al azar y las que pueden exhibir una re- 
gularidad si se les examina adecuadamente. 

La idea básica es que los sistemas complejos pueden revelar 
su estructura si se traducen los datos a un espacio de las fases 
con un número apropiado de dimensiones. 

Estos métodos se han aplicado con resultados alentadores a 
los estudios sobre poblaciones y otros sistemas biológicos, so- 
bre fluidos turbulentos, y ahora también para estos modelos 
económicos. Entre los más poderosos figura el de los algorit- 
mos que utilizan una sola magnitud (población, temperatura, 
precio, o lo que sea que nos brinde una serie suficientemente 
larga de datos) tomada a intervalos de tiempo regulares. 

En su forma más elemental, esta serie temporal es una lista 
de cifras que representan los datos experimentales, por ejem- 
plo, una población de bacterias medida cada hora: 


0.453; 0.671; 0.632; 0.661; 0.702; 0.799; 0.530; 0.501, ... 


Deseamos comprobar si aparece un atractor en un espacio de 
las fases que tenga como mínimo tres dimensiones. Por todo 
lo que hemos visto hasta aquí, se requiere entonces represen- 
tar las mediciones de tres variables independientes, pero los 
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datos que disponemos sólo nos dan información sobre una 
variable. ¿Hay solución para esta situación? 

Los matemáticos D. Ruelle y N. Packard hallaron un truco 
que resuelve el problema, y F. Takens logró demostrar que este 
artificio es matemáticamente correcto. 

El método consiste en fabricar otras dos series con los mis- 
mos valores de la medición, pero desplazados en el tiempo. 
Llamando X,, X,, Xz,... a los sucesivos valores obtenidos, se 
construyen las seudoseries: 


Y, = Xy Yo=X,, Y; 
L,=X3, Lo= Xy Za 


O sea que, en lugar de una sola serie temporal, se tienen la 
original X y dos copias desplazadas uno y dos lugares en el 
tiempo, Y y Z: 


X=0.453, 0.671, 0.632, 0.661, 0.702, 0.799, 0.530, 0.501,... 
Y =0.671, 0.632, 0.661, 0.702, 0.799, 0.530, 0.501, ... 


Z=0.632, 0.661, 0.702, 0.799, 0.530, 0.501, ... 


De modo que para el tiempo 1= 1 se representa en un espacio 
de las fases tridimensional el punto X= 0.453, Y = 0.671, Z= 
0.632; para t= 2, el punto X= 0.671, Y = 0.632, Z= 0.661, y así 
sucesivamente. 

Si el resultado es que los puntos quedan agrupados en una 
región limitada formando un atractor, esto es indicio de que 
hay una regularidad, algún comportamiento periódico. Si, en 
cambio, quedan esparcidos más o menos uniformemente por 
todo el espacio de las fases, puede tratarse de un fenómeno 
aleatorio, o puede ser que se requiera un espacio de más di- 
mensiones donde sí aparezca un atractor. 

Es claro que cuanto más elevado sea el número de dimen- 
siones necesario, tanto más dudosa resulta la presencia de un 
atractor. 

Una de las dificultades en la aplicación de estos métodos 
radica en que, incluso para sistemas dinámicos como los flui- 
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dos turbulentos no es fácil discernir cuáles son las variables 
pertinentes. Esto vale todavía más para los sistemas que se estu- 
dian en campos alejados de la física, donde es aún más difícil 
determinar la dimensión del espacio de las fases. 


E O AN 


X. Orden, caos y complejidad 
en las matemáticas 


HEmoOs visto algunas de las sorprendentes cualidades de la ecua- 
ción logística descubiertas cuando R. May la aplicó al campo 
de la biología y que causaron que muchos matemáticos co- 
menzaran a estudiarla en detalle haciendo largos cálculos con 
sus computadoras. Se trata de un nuevo modo de investigar las 
leyes de las matemáticas que difiere de los planteos donde todo 
el conocimiento se alcanza mediante pasos lógicos dentro de 
un marco abstracto y no mediante experimentos. Gracias a las 
computadoras se hace hoy experimentación matemática y esto 
es particularmente así en la que se aplica a la dinámica por el 
papel que desempeñan en ella los procesos iterativos. 


LA ECUACIÓN LOGÍSTICA REVELA 
UN MUNDO MUY COMPLEJO 


Una ecuación simple como la logística revela aspectos insos- 
pechados hasta que se hacen varios millones de cálculos que 
indican una relación muy íntima entre la dinámica de sistemas 
complejos y la estructura del sistema de números. Desde el 
momento en que las iteraciones para este tipo de ecuaciones 
revelan una sensibilidad a los valores iniciales se comprueba 
que no sólo los sistemas del mundo físico pueden ser a la vez 
deterministas e impredecibles sino que también es así en estos 
sistemas del mundo matemático. 

Solamente se podría eliminar esa imposibilidad de predecir 
un resultado a largo plazo si se midiese con infinita precisión 
en el mundo físico o si en el mundo matemático se calculase 
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utilizando todos los infinitos dígitos de los que están formados 
la mayoría de los números. 

Examinemos ahora los asombrosos resultados que provienen 
de hacer cálculos con la ecuación logística para un K mayor 
que 3.5 donde aparece el comportamiento complejo (véase la 
figura X.1). Vemos que el periodo T'se bifurca en cascada cada 
vez más rápido al aumentar K, y que la distancia entre los va- 
lores correspondientes de X se va haciendo cada vez menor. 
En una zona a partir de K= 3.57 se tiene para cada valor cre- 
ciente de K sucesivamente 1 024, 2048, 4096,... periodos, y ha- 
ría falta un microscopio para distinguir la estructura formada 
por los puntos. 

Otra propiedad notable es la de “renormalización”, descu- 
bierta por M. Feigenbaum en la década de 1970: para bifurca- 
ciones suficientemente altas, por ejemplo con 2048 periodos, 
si éstos se vuelven a bifurcar a 4096 repiten la estructura de los 
2 048, siempre y cuando se representen con un aumento muy 
preciso en las escalas; el aumento en K debe ser 4.66920166... 
y el aumento en X debe ser 2.502908... Estos números de Fei- 
genbaum son universales, como TT, porque la misma estructura 
de bifurcaciones en cascada y los mismos números de Feigen- 
baum aparecen también en otras ecuaciones siempre y cuando 
sean funciones continuas de Xy con un solo máximo. 

Las bifurcaciones en cascada y los números de Feigenbaum 
aparecen no sólo en los cálculos que hacen los matemáticos con 
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sus computadoras, sino también cuando se representan mate- 
máticamente muchos comportamientos de la naturaleza. 

Representemos en una gráfica más detallada la variación de 
Xen función de K (véase la figura X.2). Veremos que más allá 
de las bifurcaciones, para K entre 3.55 y 4 aparecen franjas ver- 
ticales cubiertas de puntos que corresponden a los infinitos 
lugares donde podría estar el sistema para un valor de K, de 
modo que si la computadora arroja un valor correspondiente a 
un punto X = 0.57739, por ejemplo, el siguiente estará en 
algún lugar entre los extremos superior e inferior de la franja 
oscura, y esto es todo lo que podemos anticipar hasta que la 
computadora nos revele el resultado del nuevo cálculo. 

Más allá de K= 4 los resultados no muestran estructura algu- 
na entre X= 0 y X= 1. De modo que los puntos son, en esta 
zona, indistinguibles de los que se hubieran podido marcar por 
azar arrojando dados, a pesar de que estamos utilizando una 
ecuación matemática perfectamente determinista. 

No podemos atribuir esta imposibilidad de predecir a la exis- 
tencia de factores desconocidos, porque la ecuación no tiene 
ninguna ambigúedad. El motivo es que como hemos visto, no 
podemos medir ni representar un estado presente con pre- 
cisión infinita, y si bien esta limitación no es importante en 
muchas situaciones, ahora percibimos que abundan los com- 
portamientos no lineales, donde causas insignificantes que 
pensábamos no afectarían a un sistema dinámico, producen 
efectos tan desproporcionadamente grandes que cambian su 
comportamiento de manera imprevisible. 

Podemos considerar que la región desde K= 3.55 y K=4 
corresponde a la complejidad y, a partir de K= 4, al caos con 
puntos que llenan todo el espacio. 


VENTANAS QUE SON COMO ISLAS EN UN MAR CAÓTICO 


En la zona de complejidad se pueden observar curvas más os- 
curas, que evidencian una mayor concentración de puntos, o 
sea, mayor probabilidad de encontrar el sistema, y además 
algo especialmente significativo: obsérvese que hay bandas ver- 
ticales blancas intercaladas, cruzadas por unas pocas líneas. Es- 
tas reciben el nombre de ventanas, e indican el regreso a un 
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comportamiento cíclico regular, donde aparecen dos o tres 
periodos diferentes. Tales intervalos con pocas frecuencias de 
oscilación son llamados de intermitencia. 

Si se examina una ventana de tres periodos como la de 
K = 3.83 con mejor definición, es decir, con más cálculos 
de computadora y se representa la zona central que hemos 
encerrado en el recuadro A de la figura X.2, en una escala 
amplificada, nos encontramos con otra gran sorpresa: el perio- 
do vuelve a bifurcarse, y así inicia una cascada con duplica- 
ciones sucesivas de 6, 12, 24, 48,... repitiendo en miniatura el 
esquema original, con ventanas que a su vez vuelven a mostrar 
bifurcaciones en cascada (recuadro B), en las que hay otras 
ventanas con sus propias cascadas y así tantas veces como se 
desee continuar explorando este universo desconocido. 

Se tiene así el mismo fenómeno de autosemejanza que ha- 
bíamos encontrado en las figuras fractales. 


CONJUNTOS DE MANDELBROT Y PLANOS COMPLEJOS 


Volvemos entonces al concepto de la generación de los frac- 
tales mediante el proceso de iteración. Cuando el matemático 
B. Mandelbrot revisó los trabajos de Gaston Julia, un discípulo 
de H. Poincaré, sobre cálculos iterativos con números comple- 
Jos, decidió que éstos indicaban un camino para construir las 
figuras fractales a partir de ecuaciones matemáticas. 

Es necesario aquí pasar revista al concepto de número com- 
plejo. Como sabemos hay varias clases de números, llamándo- 
E números naturales a los mas básicos y elementales: 0, 1, 2, 

. Éstos se pueden sumar o multiplicar para producir otros 
números naturales. Pero si se quiere restar con un método sis- 
temático, es conveniente introducir los números negativos —1, 
2,3)... 

Hasta aquí disponemos sólo de números enteros, y esto pue- 
de crear problemas cuando tratamos de dividir un entero por 
otro, de modo que necesitaremos las fracciones o números ra- 
cionales 1/2, -1/2, 1/3, -1/3, 3/2, -3/2,... y los números 
irracionales, que requieren infinitos dígitos para ser expresados, 
como T, y v2. 

Todos estos diversos tipos de números forman el sistema de 
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números reales. Pero queda aún una limitación: si se desea 
obtener la raíz cuadrada de un número, sólo se puede hacer 
con los que son positivos; así, la raíz de 9 es 3, que multiplica- 
do por sí mismo regenera el 9, pero la raíz de -9 no existe, ya 
que -3 x -3 también da 9. Las raíces cuadradas de números ne- 
gativos se denominaron, en consecuencia, números “imagi- 
narios” para diferenciarlos de los “reales”, pero como resulta 
sumamente conveniente utilizarlos para los cálculos, se resol- 
vió la dificultad “inventando” una raíz cuadrada llamada 2 para 
el número negativo —1. Resulta entonces 1? =-1 y no nos preo- 
cuparemos por la naturaleza matemática de ¿. Como cualquier 
número negativo, —a se puede escribir como —1 x a, su raíz 
cuadrada, v=a se puede escribir como V=1 x V=a, o sea iva. 
Con este subterfugio se puede considerar así a estos números 
tan reales como los “reales”, lográndose una mayor versatili- 
dad para el sistema numérico. 

Finalmente, si combinamos números “reales” con “imagina- 
rios” resulta un número que llamaremos complejo, no porque 
sea complicado, sino porque tiene varios componentes. Un 
número complejo es de la forma: 


ZL= X+1Y 


Donde Xe Y son números ordinarios reales, e 2 es la raíz cua- 
drada de —1. 

Similarmente a los números reales, se puede visualizar los 
números complejos representándolos mediante un gráfico en 
un sistema de coordenadas, pero en un plano complejo, don- 
de el eje vertical es ¿Y y el horizontal es X. El punto Zse ubica 
en la intersección de las líneas paralelas a los dos ejes que par- 
ten del valor de X y del de ¿Y (figura X.3). 

Los números complejos tienen su propia aritmética, álgebra 
y análisis, y al graficar el resultado de hacer operaciones ma- 
temáticas con ellos, se descubren aspectos de gran importan- 
cia y extraña belleza. 

Mandelbrot partió de una operación iterativa muy simple: 
asignar un valor inicial a la variable compleja Z, elevarla al cua- 
drado y sumarle un número constante que llamaremos C: 


Z=Z*4+C 
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FIGURA X.3. 


El Zasí obtenido se vuelve a elevar al cuadrado y se le suma 
el mismo C, y así sucesivamente, operación que requiere de 
una computadora para descubrir sus propiedades, pues éstas 
se ponen de relieve cuando aparece en la pantalla el resultado 
graficado en el plano complejo (figura X.4). La figura en for- 


Parte imaginaria 


1 -0.5 


FIGURA X.4. 
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ma de corazón, llamada “conjunto de Mandelbrot” representa 
a todos los valores de Zque no llegarían a tener un valor infini- 
to aunque se hicieran infinitas iteraciones. Rodeando la figura 
están todos los puntos que sí tienden a un valor infinito a me- 
dida que se van haciendo cada vez más iteraciones. 

La frontera entre ambas zonas tiene propiedades de autose- 
mejanza y una complejidad tan grande que sólo se puede cap- 
tar cuando se instruye a la computadora para que las sombree 
en distintos tonos de gris, o mejor aún en distintos colores, se- 
gún las distintas velocidades con que va creciendo 2 al iterar. 

Hay un libro clásico para admirar estas figuras, The Beauty of 
Fractals, de H. Peitgen y P. Richter, y programas de computado- 
ra como el Fractint, que han divulgado la belleza de dichas for- 
mas, en las que cada conjunto de Mandelbrot tiene adosadas 
figuras similares más pequeñas que observadas con la amplifi- 
cación de la computadora como si fuera un microscopio se 
repiten y repiten en forma autosemejante. 

Si se sigue aumentando la amplificación, aparecen otras for- 
mas características, que parecen dendritas, torbellinos, colas 
de caballos de mar, que a su vez se repiten en detalles cada vez 
más microscópicos, proceso que no tiene fin, ya que se podría 
iterar infinitas veces y, además, variando la constante (se pue- 
de tener una variedad infinita de formas (véase la figura X.5). 

Aun el observador más experimentado no deja de admirarse 
ante este proceso, mediante el cual partiendo de una ecuación 
simple se llega por iteración a un resultado con estructuras tan 
extraordinariamente complejas, lo que recuerda al resultado 
de iterar la ecuación logística. En realidad hay algunos aspec- 
tos semejantes entre ambas expresiones matemáticas, que lle- 
van a considerar que la ecuación de Mandelbrot es la versión 
en variables complejas de la ecuación logística, donde ahora, 
en lugar de la variable X'se aplica la Z, y en lugar del factor de 
crecimiento Kla constante C. 

Así como representando gráficamente la variación de X al 
aumentar Kse obtuvieron las sorprendentes formas de la figu- 
ra X.2, la representación de Zal ir variando C produce estas fi- 
guras de Mandelbrot. En ambos casos se tiene una región “or- 
denada” con valores estables, que es la región para K menor 
que 3 en la ecuación logística, correspondiente al interior del 
conjunto de Mandelbrot. La zona de “caos”, para K mayor que 
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FIGURA X. 


4 en el primer caso se puede comparar a la región más alejada 
de la figura de corazón del segundo, donde los valores de Z se 
van velozmente a infinito, y una región intermedia de una 
complejísima estructura, con propiedades de autosemejanza. 

Como hay ventanas para la ecuación logística, por ejemplo 
para K= 3.84, con un ciclo estable en el que aparecen réplicas 
en miniatura de la figura original, también el conjunto de Man- 
delbrot muestra para la posición con parte imaginaria ¿Y = 0, y 
parte real X = -1.75, una pequeña isla en forma de corazón 
(figura X.4), que, vista en detalle resulta ser un nuevo conti- 
nente semejante al conjunto de Mandelbrot principal con “bro- 
tes” que, examinados con el microscopio resultan ser réplicas 
diminutas del conjunto inicial, y que a su vez exhiben otros 
brotes más minúsculos. 

La figura de Mandelbrot y cada disco circular adosado co- 
rresponden a una órbita periódica particular: la de forma de 
corazón al periodo 1, el disco más grande al periodo 2, segui- 
do por discos de periodos 8, 16,... 

Tal complejidad nos muestra que, al igual que lo que obser- 
vamos en muchos fenómenos de la naturaleza, el comporta- 
miento complejo puede aparecer aun con leyes simples. 

Uno de los más importantes resultados de estas investigacio- 
nes matemáticas es que, con base en ellas, está surgiendo un 
nuevo método para encarar los estudios de los sistemas diná- 
micos: la dinámica compleja, que no significa complicada, sino 
basada en los números complejos. 


XI. ¿Se puede domesticar el caos? 


Hace ya más de un siglo que los matemáticos sentaron las ba- 
ses para estudiar los sistemas no lineales. Desde hace tres déca- 
das se les aplica para investigar los fenómenos Hamados caóti- 
cos que aparecen en la naturaleza. 

La actitud ante el caos ha sido equivalente a la que se tiene 
con las enfermedades: se investigan sus causas para evitar que 
aparezcan, ya que lo que uno desea son procesos predecibles, 
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comportamientos sin sorpresas, sistemas domesticados perfec- 
tamente controlables. El caos implica catástrofes, y se trata de 
mantenerlo a una distancia prudente. 

Pero también esta visión del caos está cambiando en los últi- 
mos años, ante la evidencia de que se le puede controlar y has- 
ta lograr que haga cosas útiles. 


EL CAOS PUEDE SER ÚTIL 


En la Primera Conferencia de Caos Experimental realizada en 
Estados Unidos en octubre de 1991 se presentaron informes 
que así lo demuestran, y que se refieren a investigaciones en- 
caminadas a estabilizar los ritmos cardiacos, controlar las os- 
cilaciones de los relojes químicos, aumentar la potencia de los 
rayos láser, sincronizar la salida de circuitos electrónicos. 

En todos los casos, los resultados de aplicar el caos y de con- 
trolarlo han sido muy alentadores. 

Muchas de estas aplicaciones han sido iniciadas por los físi- 
cos William L. Ditto, del Instituto de Tecnología de Georgia, y 
Louis M. Pecora, del U. S. Naval Research Laboratory, quienes 
ven dos razones fundamentales para que el caos sea de utili- 
dad. En primer lugar, el caos determinista de los sistemas di- 
námicos es, en realidad, una compleja estructura de muchos 
estados ordenados, ninguno de los cuales predomina sobre los 
demás, a diferencia de un sistema ordenado que tiene un úni- 
co comportamiento. 

Los investigadores han demostrado que si se perturba de la 
manera adecuada al sistema caótico, se le puede estimular para 
que adopte uno de todos esos comportamientos ordenados. 
La gran ventaja respecto de los sistemas ordenados clásicos es 
su flexibilidad superior, debido a que pueden saltar rápidamen- 
te de un comportamiento a otros de una amplia colección. 

La otra razón es que si bien no se puede predecir su com- 
portamiento, éste está determinado, por lo que si dos sistemas 
caóticos prácticamente idénticos del tipo adecuado son guia- 
dos o movidos por la misma señal caótica, ambos van a tener 
el mismo comportamiento aunque éste no se pueda predecir 
de antemano. La sincronización caótica es la mejor demos- 
tración de que se trata aquí de caos determinista, ya que si se 
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tratase de azar puro los comportamientos no tendrían por 
qué coincidir. 

Dos sistemas caóticos aislados no pueden estar sincroniza- 
dos, porque aunque sean prácticamente idénticos y comien- 
cen a funcionar al mismo tiempo, inmediatamente sus minús- 
culas diferencias serán amplificadas y los harán diverger cada 
vez más. Pero si son guiados por una señal caótica única ambos 
tendrán idéntico comportamiento caótico, con la condición 
de que sean estables, o sea que si se les perturba un poco, su 
trayectoria en el espacio de las fases cambie sólo un poco. Si se 
cumple esta condición, ambos sistemas suprimirán cualquier 
diferencia entre ellos y actuarán de manera caótica, es decir 
impredecible, pero sincronizadamente. Esto se está aplicando 
en sistemas de comunicaciones y procesamiento de señales. 


ÉL CAOS ES CONTROLABLE 


Las aplicaciones del control del caos se basan en el método 
ocY (desarrollado por Ott, Grebogi y Yorke, de la Universidad 
de Maryland), quienes lograron así que un sistema que presen- 
ta toda una colección de oscilaciones periódicas pase a tener 
sólo una. 

El método consiste en obtener en el espacio de las fases la 
información de una trayectoria cuando pasa por la sección de 
Poincaré del atractor extraño, y esperar hasta que la trayecto- 
ria vuelva a pasar por la sección. Si lo hace en las cercanías de 
una órbita periódica deseada, en ese preciso momento se per- 
turba al sistema dinámico modificando al parámetro adecua- 
do, hasta lograr que quede en esa Órbita. 

Lo que se comprobó es que como los sistemas dinámicos 
caóticos son tan sensibles a las condiciones iniciales, reaccio- 
nan muy rápidamente y con gran versatilidad a este control, de 
modo que su uso puede resultar muy ventajoso. 

Esto se ha aplicado al control de fluctuaciones caóticas en la 
intensidad de sistemas láser, logrando mayor flexibilidad y 
estabilidad, con una potencia de salida aumentada en un fac- 
tor 15. 

Otra aplicación de importancia es el control del caos en un 
sistema biológico: utilizando una parte del corazón de un co- 
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nejo se le provocó una arritmia, y en esa situación se le estimu- 
ló con señales eléctricas producidas con el método oGY. Estas 
bastaron para restablecer un ritmo cardiaco normal, lo que ha 
alentado a los investigadores para que busquen el control de la 
arritmia en corazones humanos y proyecten desarrollar marca- 
pasos y defibriladores basados en este caos controlado. 


Conclusiones 


HEMOS pasado revista a los múltiples aspectos que presenta el 
comportamiento de los sistemas dinámicos no lineales y a los 
diversos métodos que se utilizan para su estudio. 

Las aplicaciones de conceptos como la autosemejanza, los 
fractales, los atractores extraños, han despertado un gran en- 
tusiasmo en investigadores de las más diversas disciplinas, a 
partir de sus notables resultados para los sistemas físicos no li- 
neales. 

Se está dando un gran impulso a la mecánica estadística, al : 
estudio de las transiciones de fase, de los vidrios de espín, de la 
turbulencia, temas que además abren nuevas perspectivas en 
otros campos científicos. 

Podría ser, en efecto, que estos conceptos sean igualmente 
válidos para otros fenómenos complejos tales como las fluctua- 
ciones en las poblaciones de insectos o animales, en la econo- 
mía o el comportamiento del cerebro, etcétera. Pero para 
averiguarlo hay que saber cómo caracterizar el caos del sistema 
estudiado, y eso resulta tanto más difícil cuanto mayor sea el 
número de variables independientes en juego. Aun para sis- 
temas físicos como los fluidos turbulentos con muchos grados 
de libertad no se comprende del todo cuál es el papel del caos 
determinista en las diversas transiciones entre orden y caos. 

Cuando hay pocos grados de libertad en juego, se buscará 
un atractor extraño, pero esto se hace más difícil cuando se 
aplica a la búsqueda de atractores para el cerebro o el merca- 
do financiero, lo que exige entonces prudencia, porque las ex- 
pectativas generadas en este campo han fomentado una espe- 
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cie de moda, por la que hay casos en los que los mismos traba- 
jos científicos que hace unas décadas se presentaban con resul. 
tados obtenidos aplicando las herramientas clásicas se hacen 
ahora derivándolos de dimensiones fractales, sin que se vea una 
ventaja conceptual para la interpretación de los fenómenos 
estudiados. 

De todos modos, ésta es una disciplina muy joven, en cons- 
tante evolución. Las perspectivas son muy buenas, pues se está 
comprendiendo que, en lugar de evitar la no linealidad y la com- 
plejidad, se les puede emplear para disponer de sistemas más 
flexibles, rápidos y de comportamientos inesperados que pre- 
sentan una amplia gama de posibilidades. 

Pero las implicaciones del caos determinista para el campo 
del conocimiento se extienden más allá de su utilidad para es- 
tudiar los distintos tipos de sistemas dinámicos, y de las nuevas 
técnicas que se van desplegando a medida que crecen sus apli- 
caciones. Otro aspecto significativo está referido a la forma en 
que los cambios de nuestro concepto del mundo están condi- 
cionados por nuestras creencias, tema estudiado por Thomas 
S. Kuhn en La estructura de las revoluciones científicas. 

Al producirse un cambio de paradigma, los científicos ven 
cosas nuevas y diferentes en el mundo de la investigación, in- 
cluso observando con instrumentos y en lugares ya conocidos. 


Un ejemplo notorio es el de la forma en que se percibió la . 


trayectoria de caída de los cuerpos antes y después de los estu- 
dios de Galileo. Como nos muestra el historiador de la ciencia 
Pierre Thuillier en su libro De Arquímedes a Einstein, los sabios 
anteriores a Galileo estaban educados en la física de Aristóte- 
les, según la cual un cuerpo sólo se mueve mientras esté so- 
metido a una fuerza. Por lo tanto cuando se arroja una piedra, 
o se lanza una flecha, o cuando un cañón dispara un proyectil, 
el cuerpo se mueve en línea recta por el ímpetu recibido, has- 
ta que éste se agota por la resistencia del aire, momento a par- 
tir del cual cae verticalmente hacia el suelo, sitio en el que le 
corresponde estar según el orden natural. Es debido a eso que 
los grabados de manuales de artillería escritos por los expertos 
alrededor del año 1550 muestran trayectorias en forma de una 
L invertida. 

Podría suponerse que un experto en artillería debe haber 
observado alguna vez cómo caen una piedra'o una flecha dis- 
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paradas hacia un blanco, o la forma que tienen los chorros de 
agua o de vino vertidos desde un cántaro. Pero justo después 
de la difusión de las enseñanzas de Galileo, quien demostró 
que estas trayectorias tienen la forma de una curva parabólica, 
se comenzó a representar correctamente la caída de proyec- 
tiles, chorros de agua, etcétera. ¡Es como si los proyectiles, que 
antes de Galileo caían en línea vertical, pasaran ahora a mo- 
verse siguiendo una curva! 

De modo similar, un investigador de la talla de Leonardo da 
Vinci, apasionado observador de la naturaleza, hizo dibujos de 
lo que veía en las disecciones de cadáveres de animales y seres 
humanos. Pero cuando dibujó el corazón, lo representó con 
dos ventrículos y una aurícula, porque lo veía así, imbuido de 
lo que al respecto enseñaba Galeno, máxima autoridad médi- 
ca de esa época. 

En este siglo ocurre algo parecido, con un cambio de para- 
digma que comienza con la aparición de las físicas cuántica y 
relativista, y que continúa con la transformación que implica el 
interés actual por los fenómenos no lineales. 

Hasta el siglo pasado, sólo se estudiaban los fenómenos que 
obedecían a ecuaciones integrables, en particular las lineales. 
Si había interés por un fenómeno no lineal, se le transformaba 
mediante una aproximación lineal. Parecía que en la natura- 
leza lo realmente importante era la familia de los fenómenos 
lineales y que los otros eran una excepción, una rara especie 
indeseable por la dificultad de su tratamiento. Testigo de esto 
son los libros de física y de matemáticas escritos hasta hace 
Unas pocas décadas, donde los sistemas no lineales epenasss 
mencionan. 

Ahora, en cambio, se percibe que la inmensa mayoría de los 
fenómenos naturales son no lineales, y que los otros son la ex- 
cepción y no la regla. Si bien hay todavía mucha física de sig- 
nificación que puede hacerse en fenómenos lineales, cada vez 
más instituciones y científicos se dedican a la dinámica de pro- 
cesos no lineales. Antes no había fenómenos no lineales dig- 
nos de estudio, y ahora forman legión. Una vez más se com- 
prueba que la percepción de lo que cada uno tiene ante sí está 
condicionada por sus teorías y creencias acerca del mundo, 
hasta un nivel muchas veces insospechado. 

Como hemos visto, la incorporación del tema del caos impli- 
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ca un cambio de enfoque, una visión diferente a la que fue co- 
mún a los científicos hasta principios de este siglo: ya no se 
considera válido un determinismo “absoluto” como el que tan 
claramente formulara Laplace, y se acepta entonces que siste- 
mas que obedecen a leyes deterministas pueden tener un com- 
portamiento imprevisible, lo que hace inevitable su descrip- 
ción mediante las probabilidades, y además que esto es válido 
para la mayoría de ellos, incluyendo casos tan aparentemente 
simples como el del movimiento de todo objeto que esté so- 
metido a la acción de más de dos fuerzas. 

Por lo tanto, el estudio del caos ha revelado que la imprevisi- 
bilidad de este complejo mundo se puede conciliar con la exis- 
tencia de leyes naturales simples y ordenadas. 

Las matemáticas correspondientes a este campo también 
requieren un enfoque diferente, ya que la fundamental aplica- 
ción de la computadora permite visualizar, mediante imágenes 
que muchas veces son de una rara belleza, el comportamiento 
global, cualitativo, de las ecuaciones de los sistemas caóticos. 
También aquí se ve cómo ecuaciones muy simples, aplicadas 
por iteración, y con las características de la autosimilitud, dan 
origen a formas fractales asombrosamente semejantes a las que 
aparecen en la naturaleza: árboles, montañas, nubes, que 
jamás se soñó que pudieran representarse matemáticamente. 

Las matemáticas de los sistemas dinámicos caóticos son no 
lineales, y esto, debido a la consiguiente sensibilidad a las con- 
diciones iniciales, implica que, a diferencia de los que obede- 
cen a ecuaciones integrables, sólo. se puede conocer el com- 
portamiento global de estos sistemas. El ejemplo del estudio 
de Poincaré del movimiento de los planetas nos lo muestra: el 
único sistema integrable, lineal, es el estudiado por Kepler, o 
sea el que forman el Sol y la Tierra. Para lograr aproximarlo a 
un sistema lineal se excluyeron los otros ocho cuerpos que 
integran el sistema solar. De esta manera se puede obtener la 
trayectoria elíptica de nuestro planeta, pero si se quiere prede- 
cir a largo tiempo su posición con un mínimo de exactitud hay 
que agregar la Luna, que influye sobre el movimiento terres- 
tre, los demás planetas, y las interacciones gravitatorias entre 
éstos, lo que nos da un sistema no integrable, caótico, y donde 
aun los cometas que vienen desde los confines del sistema, o 
desaparecen por diversos motivos, pueden, merced a la extre- 
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ma sensibilidad a las perturbaciones, alterar a todo el conjun- 
to. De modo que no es posible analizarlo como se hace con el sis- 
tema lineal, por separación de los factores que afectan el com- 
portamiento de aquellos que prácticamente no lo hacen: todo 
factor, por ínfimo que parezca inicialmente, puede desenca- 
denar un cambio drástico del conjunto completo. 

En duraciones cortas, como son los miles de años de nues- 
tros registros históricos en este ejemplo, vale la aproximación 
lineal para los cálculos, pero cuanto más se desee explorar el 
futuro tanto más se deberá tener en cuenta el funcionamiento 
global del sistema. 
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GLOSARIO 


Abscisa: la coordenada horizontal de un punto en un sistema de 
coordenadas rectangular plano. 

Amplitud: en una onda, el valor absoluto del máximo desplazamien- 
to a partir del valor cero durante un periodo de oscilación. 

Análisis de Fourier: análisis de una curva periódica, a la que se des- 
compone en curvas periódicas elementales (senos y cosenos). 

Atractor: región del espacio de las fases de los sistemas disipativos 
hacia la cual convergen las trayectorias que parten de una deter- 
minada región, llamada cuencá del atractor. Los atractores “pre- 
decibles”, de-estructura simple,.son el punto y el ciclo límite, que 
corresponden a comportamientos periódicos y, por lo tanto, se 
representan con curvas cerrradas. 

Bifurcación: todo valor de un parámetro para el que se produce el 
cambio del número o de la estabilidad de estados y ciclos estables. 
Se dice entonces que el sistema dinámico pasa por una bifurcación. 

Bit: un dígito en la representación binaria de un número. Por ejem- 
plo 0110 es un número de cuatro bits. Es la unidad fundamental 
de información. 

Byte: ocho dígitos binarios o bits. 

Caos determinista: comportamiento irregular de un sistema dinámi- 
co cuando exhibe una gran sensibilidad a las condiciones inicia- 
les, por lo que es un sistema no integrable. Las trayectorias en el 
espacio de las fases que describen su comportamiento son en ge- 
neral muy complicadas, y se mantienen en una región del mismo, 
sin crecer ilimitadamente. Para que haya caos determinista, el es- 
pacio de las fases debe tener un mínimo de tres dimensiones. 

Ciclo: el ciclo describe movimientos periódicos que se pueden pre- 
decir, como las órbitas circulares. En el espacio de las fases, se re- 
presenta mediante curvas cerradas. 
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Ecuación lineal: para estas ecuaciones vale el principio de superposi- 
ción de las soluciones. Si 4, y u¿ son soluciones de una ecuación li- 
neal, toda combinación lineal de éstas, tal como u= au; + bus, donde 
los coeficientes a y bson constantes, es también una solución. 

Ecuación logística: la ecuación no lineal X, = KX,(1 — X,), que se 
aplica, entre otras ciencias, en la biología como modelo de predic- 
ción del crecimiento de una población. 

Ecuación no lineal: ecuaciones en las que no vale el principio de 
superposición de las soluciones. La respuesta de un sistema no li- 
neal puede ser sumamente sensible a las condiciones iniciales. 

Ecuaciones diferenciales: ecuación que, en general, tiene la forma: 


dx,/ dt= f(x)" 


donde x,(t) representa la variable : y la función f(x) da la evolu- 
. ción de x(t) con el tiempo. Los modelos matemáticos.en las cien- 
cias físicas se formulan a menudo como ecuaciones diferenciales. 

Espacio de las fases: espacio matemático constituido por las variables 
que describen al sistema dinámico. Cada punto del espacio de las 
fases representa un posible estado del sistema. La evolución en el 
tiempo del sistema se representa con una trayectoria en el espacio 
de las fases. 

Exponente de Lyapunov: una cantidad utilizada para caracterizar la 
divergencia de las trayectorias con el transcurso del tiempo en un 
proceso caótico. 

Fractales: conjuntos de puntos, o curvas, que exhiben una estructura 
cada vez más detallada a medida que se les observa con mayor am- 
plificación. Muchos fractales son autosímiles, es decir, que la figu- 
ra tiene un aspecto similar en cada escala de amplificación. Los 
fractales pueden tener una “dimensión fractal” no entera. 

Función analítica: una función que es diferenciable en la vecindad 
de cada uno de sus puntos. 

Grado de libertad: variable independiente en un sistema dinámico. 

Iteración: repetición de una operación o conjunto de operaciones 
matemáticas. 

Notación binaria: la notación binaria es la representación de los 
números enteros en términos de potencias de 2. 

Número complejo: un número conteniendo una parte real y otra 
imaginaria, y de la forma A+ Bi, donde ¿= SL. 

Número irracional: un número real que es no es racional, o sea que 


160 


no es el cociente de dos números enteros. La mayoría de los nú- 
mero reales son irracionales, lo que fue precisado por G. Cantor, 

Oscilación armónica: una oscilación cuya frecuencia es un múltiplo 
entero de la frecuencia fundamental. 

Periodo: duración del ciclo completo en una oscilación. 

Realimentación: cuando una parte de la salida de un sistema retorna 
a la entrada del mismo. ; 

Resonancia: una oscilación de gran amplitud causada por un estímu- 
lo periódico relativamente pequeño. 

Sección de Poincaré: una sección de Poincaré es la intersección de 
las trayectorias de un sistema dinámico en el espacio de las fases 
con una figura de una dimensión menor que la de este espacio. 
Así, para un espacio de tres dimensiones la sección es un plano. 

Sistemas dinámicos: sistemas en los que determinados parámetros 
evolucionan con el transcurso del tiempo. Se les describe mate- 
máticamente mediante la representación en el espacio de las fases 
de los estados posibles del sistema, y mediante las reglas de su 
evolución en el tiempo. Si la evolución de sus variables es conti- 
nua, la descripción se hace mediante ecuaciones diferenciales. 

Sistemas dinámicos conservativos: en la mecánica, los sistemas diná- 
micos conservativos, también llamados sistemas hamiltonianos, 
mantienen constante su energía. El volumen que los representa 
en el espacio de las fases se mantiene constante con el tiempo. 

Sistemas dinámicos disipativos: sistemas en los que la energía no se 
conserva, debido a pérdidas por fricción u otras causas. El volumen 
que los representa en el espacio de las fases disminuye con el tiempo. 

Sistemas dinámicos integrables: sistemas dinámicos con compor- 
tamiento regular, y para los cuales, por lo tanto, la representación 
del comportamiento en el espacio de las fases puede hacerse me- 
diante trayectorias sobre la superficie de un toro n-dimensional. 

Sistemas dinámicos no integrables: son los que se describen me- 
diante ecuaciones diferenciales no integrables. En general, su 
comportamiento en el espacio de las fases se representa mediante 
trayectorias caóticas. 

Trayectoria: una secuencia de puntos en el que cada punto genera al 

- siguiente de acuerdo con una función matemática, 

Trayectoria caótica: una trayectoria caótica en el espacio de las fases 
tiene tres características: 1) se mantiene en una zona limitada de 
dicho espacio; 2) no se estabiliza pues no pasa a un régimen per- 
iódico; 3) tiene sensibilidad a las condiciones iniciales. 
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La Ciencia para Todos 
es una colección coordinada editorialmente 
por MARCO ANTONIO PULIDO 
y MARÍA DEL CARMEN FARÍAS. 


Entre el orden y el caos: la complejidad. 
El Universo, la Tierra, el ser humano... sur- 
gieron de un estado inicial de desorden o 
caos, así lo imaginaron las antiguas Cos- 
mologías, y también la Biblia: “la Tierra es- 
taba informe y vacía, y había oscuridad 
sobre la haz del abismo, y el Espíritu de 
Dios se cernía sobre la haz de las aguas”, 
para establecer el orden. 

A lo largo de milenios el hombre ha ido 
descubriendo, a través del razonamiento 
científico, que detrás del aparente desor- 
den o complejidad de la naturaleza exis- 
ten comportamientos iterativos —que se 
repiten— ordenados y regulares, que se pue- 
den registrar, predecir, cambiar y explotar, 
esto es, son dinámicos. Como señala Gali- 
leo: “el libro de la naturaleza está escrito 
La en lenguaje matemático.” 

En esta obra el autor descñbé el lenguaje característico —ecuaciones no li- 
neales— de estos sistemas, llamados caóticos —que están entre el orden en el 
que nada cambia y el estado de total desorden—, así como los diversos méto- 
dos que se utilizan para estudiar su comportamiento. Esta nueva línea de investi- 
gación en el campo de las matemáticas, que se inició hace apenas tres décadas, 
“el estudio del caos”, emplea un enfoque diferente. Utilizando las leyes conocidas 
de la física y con la ayuda fundamental de la computadora analiza, por ejemplo, 
cómo aplicando la iteración y autosimilitud a ecuaciones simples, es decir 
a través de una representación matemática, se originan bellas formas fractales, 
como árboles, montañas, nubes, etc., muy semejantes a las que vemos en la na- 
turaleza. 

La complejidad constituye, pues, el nuevo desafío para la ciencia. En este 
campo, investigaciones matemáticas actuales se están aplicando también para 
explicar fenómenos en los procesos fisiológicos y neurobiológicos rítmicos, 
como las arritmias del corazón, y en las ciencias sociales, como los mercados fi- 
nancieros. 
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